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PREFACE . ' 

, Problems from the sciences that are inje/ited^into mathematics texts«(to s6tve as 

p * A ' y . ' > \ '.-^ • -' 

vehicles for the^eVelppnlent of theory(or as examples fpr the use of theory) often^ex- 
hibit the mathematical sciences as. short -oMer menus and applied mathjfematic^ as the 
corresponding cook-books. v)k sought to cowiter'this in the S&ISG 12th grade text 
C^lcuflf^ , ,In adcfttion to exploiting applications to motivate the development of the sub- 
ject and to illustrate some of its ^pecjal consequences, E.G. Begle (Director of SMSG), 
A. A, Blank (Chairman«df the' Calculus team), and I decided that it would be of interes.t * 
to attempjt ^ systematic development of applications. As -part of this program, I 
*prepared one chdpter (9; Growth, Decay and Competition) to^show how essentially one . 
general method of the ca^culjis is used in maily different sciences, and another (15: 
Geoijjietrical^ptics and WavesJ to show how different methods of the calculus have . 
Jurthere(f the development of essentially one field of science. 

In the present volume, I have modified ;these two chapters to make them independ- 
ept of the original text, so that they may serve as a math-science supplement to other 
programs. As such, the present material forms part of the series started by Volumes 
X and XI of the SMSG Studies* in Mathematics, volume? base'd on the lectures for high 
school mathematics teachers given by Max M. Schiffer ( Applied MatK^matics in the 
High School , Volume X) and by George Polya (Mathematical*MetHods in^cienc^ / Voi^ 
ume XI). m art attempt to provide, simila't. motivation /or the audience ajidressed by • 
Schiffer and by Polya, this volume begins with an introduction suggesting the inclusion 
of math-science material" at the secondary level; much of tbis appeared fti the Calculus * 
'(l-3j The Scope of Calculus, or in the Teache;:s Commentary), and in both books it 
h elps, prepare the' reader ,fir sub se(luent material. ■ ^ ^ . 

^ Although this volume consists largely of personal contributions to thl^ Calculus , * - 
I am very much awa,re of the stimulus an^ help that I received'from itiy colleagues, 
particularly from 'A. A. Blank (the "Al" of Chapter 2), and fron^ F. L, Elder, M. S. 
Klamkin, C. Leeds in, M. A. Linton, Jr., L Marx, R. Pollack and H. Weitzner. Lam- 
also very ple^^^ed to acknowledge that the work would not have be^^i attempted in ^Mjr^^^ 
other environment than that created by Ed Begle and SMSG. \ . 

Victor TWersky * - • 



UniverSfly of Illinois at ;Chicago Circle , ' 
Department of Mathematics, December, d9j36 
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^ » thapter 1 . 

' { '" INTRODUCTION 

How did you he^r about* calculus? 
How did you ^"ear about Hejen of iTroy ? 
\ And have you yjpt heard the story about Al?^ 



^1.^ Tl^e Scope .oi^Calculus. • . ' ^ ^ . 

C^culus^s the sAudy of the derivative and the integral, the relationship between 
these conc^Vts and their applications. The great advance which takes the calculus 
beyond algebra and geometry is based on the concept of limit. The basic limit proce* 
dure of tfte^differential calculus is typified by the problem of 'finding the^ slope of a, 
curve; the basic limit procec^rVof the integral calculus is typified by tlie problem of 
finding the, area enclose^i'by a curve. The slppe is fouild as a.iderivative, the are^a as 

* an integral, and superficially these appea^r to be unrelated. But there is only one 
cal^lus^dejrivativ^ and integral' are complementary ideas. K we take the slope of^ 

graph-of the area function, we are brought back to the purve itself. If we take the 
area und^r the grapl) Of the slope fuiiction, we find.the original curve again. The limit 
concept, in its. guises of , derivative' and'integral, together with this inverse relation • 

* between the .two, provides the fundamental framework for the calculus. ^ 

• ^ The darfvative and the integral may be mte^reted geometrically as jslope s/tid as 
'♦area, but these are only^ two among a wide range of interpretations. We usually begin 
with slope and area^'ln order to introduce parts'of the subj^t in an intuitive, geomet- 
ricalway. But ^though an intuitive geometrical introduction is useful and suggestive, 
it mtTst necessarily be based on very familiar steps/ The steps are so familiar that 
it may not seem that they could tead ultimately to entirely new methods for solving . 
completely* different problems th'^ those, encotintered in earl^r courses. We should 
therefore stress that tj^e concepts of derivative and ir^e^i^ral are universal, and their 
incorporation int^ a calculus, a system of reckoning,' enable? us to solve Significant 
problems in all branches of science. We can set the stage for a systemati<rdeveiop- ^ 
ment of the subject?, and emphasize the universality 9f the concepts by mentioning netv 
kinds ol prpblenis at tl;e very' start of the course. Then in addition to solving problems 
that are primarily vehicles for the»»development of theory or for th^ illustration of 
tedhniques, problems which do not begin to suggest the full scope of the subject, the 
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course itself should also ii^clude significant applications of the calculus. Problents " 
for which 410 methods of the earlier mathematical courses are particularly helpful can 
'be fo.und in very familiar contexts -^Droblems that range from the Spreading of rumors 
to the colors of the rainbow - and their very familiarity serves to emphasizp the 
' novelty of the metho'ds. ♦ ' "'^^ 

No methods of t\ie earlier courses help answer questions such as: How diSTyou 
first hear about calculus? How did you first hear, about "Helen of Troy?' To frame 
such questions mathematically, we must first isolate some essential features: Some- 
stories spread like fires; others die out. If the story is too dull, nobody bothers to / 
repeat it. But if the^story is good, some of the people who hear it (and remember it) 
pass it along. Starting from these Ideas, how far could one get by pre-calculus * 
methods? - , • _ 

*tk The same concepts that a.re basic to the spreading of stories are also basic to the 
processes of forgetting and learning*, So many of the facts and procedures we stuffed 
.into (fur heads and n6ver used afterwards seem to have vanished. Others, that we met 
repeatedly and actually worked with have become so much a part of ourselves "that we 
feel we have always known them. Our first^exposures to these facts may not Jiave 
taken, but repeated encounters in different contexts finally left their-mark. (We take 
tliis as the fheme of Chapter 2 J ^ ' , 

Can we jbonstruct a mathematical description for the way t1iat stories spread — x 
the way thatlwe learn, and forget? Starting from, approprriate assumptions (the math- 
^ematical model) we can discuss some limited aspects of sUch processes with the aid ' 
V of elementary calculus. Suqh processes illustrafe a broad class of phenomena wJiO^ 
unifying features are the basic math'ematical models.for growth, decay, and competi- 
tion. Besides helping to describe the.spreading of rumors, and learning or forgetting, 

these- saiiie mathematical models serve to clarify our observations of radioactive 

decAy, the attenuation of sunlight by a cloudy sky, the progress of chemical reactions, 
the growth of bacterrial colonies, or the spread of disease through a city. In each of 
''these situations, the essential feature is that the amount of some, quantity is changing 
(with respect to time, or distance, or whatever) at a ra^^ proportional to the amount 
already present. A process of this sort can be mathematically described by a cer- 
taih type of equation (a differ^tial equation) whose solutions, at Feast in the simplest 
cases, are combinations of exponential functions. » . . . 



Other processes of jiature change in a cyclic or periodic way; they repeat in 
identical form each year, each second, or every inch. The planetary motions, the 
'tides, the harmonious chords of music, the propagation of x-ray's through crystals — 
even the colors of oil films on water a*ll depend on periodic phenomena; -^For such 
processes, the rate of change of the rate of change of some quantity is (negatively) 
proportioijal to the quantity itself, and the mathematical modef leads to a different 
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class of differential equati3ns whoseosolutions, in the'simplest gases, are tombina- 
tions of trigo'nometrtc functions, , , • 

With^the calculus we may also* investigate more complicated natural processes 
that involve^a combination of grojvth or decay with some sort' of cyclic behavior. We 
may also solve much simpler problems. How much time \vill ft take to drive 300 
miles if you start at a-speed of 20 miles per hour, but incres^e your speed by 2 mlles^ 
pe^liour for each l^our of driving? At what angle sho^ld you throw a ball for it- to 
traVel as far as po&sible? In what directions with* respect to the sun are the rainbow { - 
colors strongest? * * * , ♦ 

These problems and many others which the calculus solves involve rates o\ 6hange, 
this is the province of the differential calculus, A secgnd broad vari6fy of questions 
is concerned with totality — the summing of small effelJts, this is the province the 
integral calculus. -By Recording your speed during a long trip^ Nvlth many accelerations^ 
and decelerations, can* you calculate how far tlie trip has taken you? If we know how a 
single drop'ot ink spreads on a blotter, can we predict what happens if we spill the 
whole bottle? Starting with a simple source V)f radiation^ can we predict the total , 
radiation from 4n extended "distribution of such sources? Knowing the Nvay a single' 
'droplet of wat^r perturbs a ray of §ul|ftght, c*in we determine how*much I^^gTit reache^ 
us on a cloudy day from the entire overcast sky? , 

s 

Such summation or integration problems are closely related to the rate -of -charge 
or differentiation problems: -the total effects result from an addition'of small ♦ 
• variations. Therefore we do not 'study separately an integral calculus and a'differen- 
tial calctilus. We study a calculus cbmprising' both differentiation and integration, and 
each aspect helpj; Us to undejrstand and apply the other, 

'Most.of the applications of calculus empllasize the) effects of variation or sumiiia- ' 
tion. "Calculus" was tailor-made to treat such problems. Except for the simplest 
f)roljlems of thisJype, the methods of arithmetic, geometry, Und algebra ^re inadequate^ 
and even for the simpler problems the methods of calculus are the more efficient or 
the more suggestive. 

The calculus w-as invented to treat probleins of physics. In tlirn, to the groNvth/of 
* * * « . 

the calcii^us, and to, the cieyelopments- it led to in the larger branch of mathematic^ 

knowi^ as analysis, we owe much of the progress in the physipal sciences and modem 

engineering', and^4nore recently iri'the biologltal and social sc^enceS, Jhe concepts 



and operations introduced by the calculus provide the righ! language, the right/ tool^, 
even the right reflexes 'for tjie major part of the applications of mathemajticafto the 
sciences'. '* * . ' • • . 



We Should therefore emphasize applications, in a general course, not*ust to show 
that calculus provides useful metho^lfe^concepts for the sciences^but because so 
much o{the calculus was developed ta solvfe specific problems. We should show the 
student how the effort to solve physical problems led to,mg^hods of the calculus, how 
the attempt to make the best use of these methods an^'to understand their full scope 
and limitations led to the 'development of the calculus as an independei)^ study, ^nd how 
the products of this, study in turn led'to deeper insight into the original problems. We ' 
. should show that science enriches mathematics' by providing significant problems and 
. concrete difficulties,' and thit mathematics enriches science by providing, solutions of 
problems, and system and organizatien. We^should balance accounts of physical moti- 
vations of mathematic?.! pvo/edures with mathematical formulations of physical phe- 
nomena, and match concept for concept that each has acquired from theK)ther. 

An ideal general course in calculus would maintain unbalance of topics and of 
viewpoints that would meet the requirements of students who will become ma'themati- 
cians, others who will become scientists primarily interested in applications, and 
still others for whom 'mathematics wHl become simply one of maAy cj^ep iluellegtual 
experiehq^s 'during ihe't^ education. For the stWent of science, a fluent intuitive grasp 
of mathbmatical methoas is a primary need, for the student of mathematics, a ca»eful 
deductive dei^elopmept is essential. These drffjerent viewpoints conflict on occasion, 
.but tRey also suppleiyent each other, and both thcj scientist and the mathematician gain 
by a deeper ppreciation of botli these views and of their interrelations. > 

" Histprically;41ie replacement of an intuitive basics for the calculus (the method of 
infinitesimals) by a careful log^ical structure (the method of limits) marked a vital 
phase^in the development of mathematics. This phase is far from complete. We are 
, still attempting Xp learn how to combine intuition in approaching tjew jj^oblems with, 
.the effective u5e of logic, not only to^^jnper and to verify our intuition, but to permit 
generalizations of broader applicability. Today /most mathematicians appreciate the 
'essential ^oles of l^oth intuitive W deductive procedures, not only for creating mathe- 
matics, but for learning it ahd for teaching it, and it is particulj^rly important to stress 
their interplay in^introductqry courses. ~ . ' r 

• 2. Systematic Approaches to Calculus-Science Programs . 

The primary purpose of a modern introductory course in Calculus is a systematic 
development of fundamental theory that stresses the novel featurdfe of calculus ai^ • 
indicates its relations' tp the earlier arithmetic, algebra; and gdtoietry, as^welfas to 
the later analysis. The initial level of rigor should be appropriate for the students,-^^ 

^ and the problem's that motivate the development of the subject, and the exercises that 
illustrate some of its special consequences or develop manipulative skill, should be' 

-^/meaningful to the students. After some grounding in theory and with some facility in 



techniques, we bhould attempt some systematic development of applications of the 
calculus to, the sciences, . , V * 

To serve up applications only as isolated fSrpblems and exercises provid(?s the 
student with an indigestible stew of unrelated ingi^edients,. To exhibit the, mathemati- 
cal scienc«s as. short-order menus and applied mathematics as the corresponding cook 
book recipes is criminal. To shuffle together menus and recipes and rules, of cookilig 
without mnen>onic guides yields' 9nl> an unteachJble and unleaftnable mess. When we 
teach calculus, let us teach calculus. And when we consider applications within a sys- 
tematic development of the basic concepts of thi^alculus, let us attempt to do so in a 
systematic way as concrete, steps in the ascent, . - . 

There are two different systematic approaches could follow, We could pick 
one gerfer^ method of the calculus and show how it is used in many different sciences, 
or we could piick one general topic of science and show how differeat methods of the 
calculus have furthered its^evt^opment. Because these two procedures are essen- 
tially different-, and because both are indicative of the.wa^ applied mathematicians and 
mathematical' scientists work in practice, we follow both procedures. 

Chapter 2 '(Growth, Decay, and Competition) illustrates the first approach. It 
develops mathematical mod^s and shows how they are used in va^ous sciencQ,s. To 
stress the generality uf the mathematics and that the equations are completely inde- 
pendent of any scientific discipline, /he models are introduced to describe the spread- 
* mgof stories. The equations are then applied to radioactive decay, electric circuits, 
^ attenuation of sunlight and the color of the sky, propagation of nerve impulses, for- 
getting and learning, chemical reactions, and to sociology. This' chapter sti'esses that 
despite the marked differences in the classes of phenomena studied in the various 
sciences, there ard a number of processes th^t are fundamental, to all. The same . 
mathem'atical structures, and consequently the same unc^erlying^concfepts, are ^n- 
countered again and again in different contexts of nature j The equations are the same, 
bufthe functions and vari^les represent different meastftreable quantities withdiffer- 
ent name's th^t change fro^ sj^ience.to science, ^ 

Chapter 3 (qjeometrical Optics and Waves)' illustpites the^ second systematic 
approaph to ap'plicatipns. It follows a sequence of physTtal concepts (axioms called / 
"laws' of nature") and shows how various methods of th^ calculus supplement each 
qther in revealing,the implicit consequences. It starts with yery restrictive laws, 
weakens them, and thereby generalizes the development. The 'procedure is "quasi- 
.axiomatic," The laws contain, undefined ternis (as ^9 the axioms of a ;nathematical 
system), ai)d implijeit restrictions on the associated observationajl and computational 
techniques (essentially as in mathematics wherje axioms are usually stated in an 
implicit dpntext)/ The developjnexrt suggests the hearistic search for first principles, 
ariH shows that physics, like mathematics, i^. cumulative : basic concepts thousands of 



years old are first accepted as laws and then exhibited as .special consequences of a . 

• more general, set of axioitfs. The chapter considers geometrica-i propagation and re- 

^ flection (EucUd, Hero, caustics, shadows, edge diffraction, eikonals^, refradtion (Snell, 
Fenftatr rainbows," mirages), energy conservation (Kepler -Lambert), waVes and super- 

• position (Huygens, Newton, Young, Fraunhofer, RayleiglT-Born), Rayleigh's theory for 
the color of the sky, Ffesnel diffraction and the method' of stationary phase, and some 
features of the complete mathematical model for scattering. 

•In the*course of Chapter 2 we introduce many special labels from different 
sciences. However, these are not essential to the mathematics, a^id we rarelytlabo- 
rate on them: iiigtedd we eiriphasiz^, that as far ^is the present development is con- ' 
cemed, they are merely different iem^v^xy names for the same mathematical corii- 
ponetitsr The primary aim is to illustrate how one mathematical model based on ' * 
calculus covers a variety of different phenomena. On the other hand in Chapter 3, ^e. 
seek to present a sequence of physical co«eepts in a mathematical way (the structured 
development of n(ajor topics is' ah essential mathematical feature of this chapter)'and 
to provida the physical motivation for several mathematical techniques that are basic 
to wave^hysics. We introduce these^techniques.in easy contexts where the resjiits 
can be obtained directly,.and then apply them (s6me in considerable depth) in more 
com"plicated situations . ■ ' •■ .: , 

, , ■ ' , • ' ' ' <f ' 

' Both Chapters can be fitted strufcturally in a calculus program, n6t only to illus- 
trate many elementary concept* and methods, but'to help provide motivation for sub- 
sequent developments. The first, which uses only eleinentary procedures to solve *he ■> 
equations that ari%e, provides a. lead-in to a systematic disci{gsion^ elementary dif- 
ferential-equations arid techniques of integration (and it serves this* purpose for Chap- 
ter 10 of the SMSQ-Calculus) ; similarly for the second, an^ultiple integrajg' line 
and surface integrals,' partial differential equations, etc/"^^' 

The two chaptersvpn mathematics;apd science a,re VeryLdifferent., .Both- attempt " 
systematic approaches, but the first uses matheiftatics as the guidijig thread,.^hd tire 
second usee science' as the thread. .Tjje^.chapters supplement each other in indicating 
the ways tha^ mathejijaEtics ^nd science interact, and may'sr%st crossing threads of 
a fabric. T% two,spedific threads we.follow intersect at lilyleigh's theory for th'e , 
color of the sky: in Shapt^r 2, it isia special case-of a general attenu?tion-process; 
in Chapter 3, k speciar case of a ^eiieral scattering process.' The cha'^rs together' 
may suggest that mathematics and Vcience provide the threads tfcat glVe structure. to ' 
oijr concept of nature. ■ ' - . , ' ' ■ " . V 

^ ,. / 

3. Style, Rigor, Notation, and All'TKaf. ^ , ' ; i 

' • , . , ^ « 

Parts of Chapter 2 (whichis a self-contained revision of Chapter 9 of the.SjilSG. I 
Calculus) 'are written ih a' relatively lig.bf Vein in ^ Attempt to break'th^ monotony of . 



the constant iteration of the basic material. It is easy to r^ad (at least by 12th grade 
students), and the model fox stoty spreading picks up the threads of Sectioi^ l.l^on 
**Hw did you hear about calculus? How did you hear aboufHelen of Troy?. . ei6. 
WTien we convert, the equations to a model; i^P^ forgetting anfl memorizing isolated facts) 
•we pick up other threads. About half of the appliqatioAs m6ntioned in Secti(Vi 1.1 are 
covered in the next chapter, anS many of the others are covered iri the'tinal chapter. 
The sizfe of the space oV the^dmplexity of the equations allotted'to the^iffer^^t dis- . 
"ciplines touched on in Chapter 2 is not meant to correspond to their irapoiianoe or to 
the writer^s preferences (except foT story spreading). - 

If Chapter 2 represents a set of hors d^oeuyres with ^ compon flavor, then Chap^ 
ter 3 (which is much harder)' represtots a sub'stantial dinner w ith a structured se- ^, 
: quence of dishes, '^^i.s material (in somewhat jjJifferent form) was initfally assembled 
to provide the introductor> Jectures f of the writer^s graduate courses in Scattering 
Theory in the Mathematics Departnaents of Stanford University and of the Technion- 
IsraeHifstitute of Technology. The different backgrounds of the students^ (mathemat-* 
ics,^gh^ics, and eiectrical,engineering)'call^d for a general elementary survey to 
introduce fundamental concepts'and terms,: and this a]so served to' providejmotivation 
foi^the mathematical model of scattering, for considering certain classes of problems, 
and for th^ development of both analytic and heuristic procedures. The present 'version 
Chapter 3, is in generat more elaborate; it$was amplified initially to serve as a larger 
vehicle |or calculus teclAiques, (as^Chapter 15 of the SMSG Calculus) , and ,then revised 
to make it more self-contauied., Ho^^^SverTthe final portion of Section 3.9 is largely a 
^top-gap for a developmWtljjeft should l^ad into Greenes theorem. Maxwell's e^tions, 
the acoustic ^equations, and Schrodinger^s equation. , * ^ . ^ . 

Chapter 3, which attempts to give an individual esthetically complete picture of ar ( 
portion of wave physics, is hindered because it Ijas only a quii^ liniited number of ' 
physical concepts and mathematical methods^available for the development. We pro-^ 
'ceed heuri^tically in some cages (sometimes the ^rigorous procedure is mentioned 
parenthetically), and in other cases we visualize experiiRents and anticipate refine- 
ments in attempts to pr^ide some substitute fo'r the intuition -that comes with first , 
Imjid experience. 'It is hard to be both- simple and honest, and to avoid-using implicitly 
tfye more complete results that are available. Attempts have been made to wee^ out 
blunders and oversights, but tSj^orms'of various sections are still only provisional 
versions of what the writer would like them Vo "be. In any case,%iuch of the mathe- 
.matical deyelopn^ent^is informal, and some discussions^ (probably many moi:©-than 
Intended) are n6!^more than suggestive,* ^ '^^ * : 

As for notation, the following wamifag for the writer's material was included (by 
another member of the team) in the Calculus; Teachers t^ommentary (p. 6^5): 



»an -scientific usage there are often notational ambiguitjie^ which we have not 
eliminated. When \^ say N is a function of t we mean that there Sxists a function 
V t —N, not that N' itselt is a function. Careless parap^irase of this expression" 
leads to the use of N(t) for v(t), a particularly bad but comm(Vi.notation in science., 
'While at first encounter, such ambiguiti^es may distress the, student, he^ay take' 
comfort in the realization tfiat comprehensioli of this and future scientific material 
he rftay r^ad can be enhanced by less-cumbersome — if lesV precise -notation " 

A les's specifi(v^v<^in'gon notatv)n that is also appropriate was written by the 
late I, Marx and terturlfed by the ^vVit^r in 1964- ' ' / / ' * * 

^ ''Whei;! a function f is defined by a formula such as fix) = x^(^ - 4x)/ we, distin- 
guish between the function f, an assignment pf val^ues tp^all x.in^e dbn^ain Of f, and 
the values f/x), the numberls assigned. Si^i^rly, in Geometry >Ve distinguish between, 
a point P' and the coordinates (3,-2) of the point, an4bet\veen a cur^ C and the 
equation (x + 1)2 = 2(y 2) of the curve, noweyevl for brevlty''and?^^heLre confusioji 
is unlikely/ive shall omit certain woydi^ (which shbuld be included^ for mathei^aUoal^ 
accuracy) with the Understanding that they are in/plicit. We^shall shorten ^tthej^- 
tion f with Values f(x) = y?(72 - 4x)t^^o 'the funbtion f(x) = x2(72 - 4x)U 'thepoiht a 
P with coordinates (3,-2yv to »the popif (3,-2)y and 'the cu^ve_C"wLth equation • 
'(x+ 1)2 = 2(y'-'2)» toVhe curve (x Vl)^~= W- 2)», In redding such' phrases, we\eep 
in mind that the mathematical pUject ari^ its Dtesentation by a formula or set of num- 
bers are quite distinct: the representation id used to replace the object in order to 
make the material easier to read- Of t^oursi^ such al5'^reviations are usually used 
conversationally by professionaymathemat/cid^ but the^^ avoid them when they write 
fbr mathematics journals:^ they/avc^id the 



(wJiicK would make the article Absolutely, 
^ords j^hich makes the articfe only Coi 

4^ Collateral Reading sLndr^Mehrhnoei 
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[qX by putting in all the extra words 
eadable) , but^by leaving out almq^al^ ^ ^ . 
lonally unreadable).''^ 



Those who may not yetliave r^d^"" preceding books of the SMSG math-sciende 



series should read the excdlent an^rMvily volume^'by 



lax M. Schiffer, Applied Mathemati 



■nfhe High SchooK Vol. X of SMSG Studies in 
' • » * ' - ' 



Mathematics, 1963; ;;. 

George Pojya, Mathematica.1 Methjc^s in Science, VoL XI of SMSG Studies in Mathe- 
matics, 1963. ll ; ^ ' 

Part of the following discdssio|is.of\i^of)Ulation. growth and Euclidean optics ia based 
on Volume X (which assumes no caldiiluL and uses difference equations instead of dif- 
ferential equations). Both volumes ai^.e rVmarkalJle for their content, style, and for the 
insight they convey^ 
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' " ' • 

Practicalljxe^erything'in the following chapters is included in Calculus , Sttident 
Text plus Teachers Commentary (SMSG Calculus Team, Chairman A, A, Biank), 
SMSG, 4966,' ' • * . ' ' , > . 

^ |bme 6f this text^s alternative versions of the present material may well prove 
more palatable, and some, of the results that are simply stated in the present volume 
- were assigned as'problems in the Calculus oftith worked solutions included in the 
Teachers Commentary. This text also includes H, Weitzner^s Chapter on mechanics 
,and oscillatory phenoxnena, many additional problems that extend the utility of the 
present chapters,, and much acjtiitional material of interest to math-science programs, 

Many of .the physics topics touched on in Chapter 2 are covered in a good general 

college text on 'physic's that juses calculus^ 
' ; ^ ; / '4 

. F,\V, Sears and Mark W. ^mansky, Umversity Physics , Addison-Wes'ley Pregs, Inc., 

' j^ass., 1952. 

..^ "More results onj ral^ioactil^e disintegration are given by 

Jienry Sematr, Atomic Phvsics ," Rinehart and Co,, Inc?, New York, 1946. " ^ 
Chemical jreaC.tions are discussed in detail by 

S. Glasstone, Physical c|igmistry , Ch. XI, D, Van Nostrand Co., Int,, New York, 1940, 
The model for learnin^l is* based on ^ • ^ ' ' ...^^^ 

H. Von Foerster, Ql ^ntum Theory of Memory , Transactions* of Sixth Conference 

Cybernetics, pp. 112-134/ Josiah Macy, Ji^. Foundation, New York, 1956. ^ ^ 

N ' . - 

\ The Section on propagation of nerve impulses is based on the work of Blair and 

Rashevsky as discussed by , ^. * '* ^ *, 

' N. Rashevsky, Mathematical'.Biophysics , Ch. XXIII, University of Chip^o Press, , 
\ • / 'hlindis, 1948, f ^" ^""^ 

and the discussion M. Socio lofflTwa^cfeived from 

N. Rashevsky, Mathematical Theory of HUman Relations ^ Principia Press, Indiana, 

As for Chapter 3, an exciting, nom-mathematical book on natural visual phenom- 
ena' is that by 
Jo. ' / 

M. Minnaert, Light and Color in the Open Air , Dover Publications, Inc., New York^ 
1954, ^ . ••^ / ' * * ^ 

• — ;and^>a*ctetailed histcTrical dewlopment of optics is given by ^ ] j 

EmsTl^Tach, The Principles of Physical Optics , Dover 'Publications, Inc., New*Xork,* 



A good general cfollege tekt 6i^.geometrical and waye optics is . 

F. A, Jenj<ins a4id H, E. White, Fundamentals of Optics . McGraw-Hill Bodk Co,, Inc, 
; New York, 1957,. . ! ^ ' * ' ' ' 

A detailed discussion of tl^^ wave fronts ^or reflection from a cpnckve hemisphgVe is . 
given by / , • « . * 

R, W, mod, Phygical Optfc^ , p, 54ff, The MacMillan Co., New York; 1*934, 

*/ *• ' ' " . ' ^ »> * 

Detailed discussions of y^r/ous aspects of wave physics are giveliby" 
/ • 
^ R. b/ Lindsay, MBc6anfc|X''Radiation , McGraw-Hill,^ ^w Yt>rk, I960; ' 

Photographs of the caus|i|s oredge waves (or edge rays) are ^veri by' 

3/. CoulSon and G,'G. Becknell, Reciprocal Diffraction Relations Between Circular and 
Elliptical Plafites,^ Physical Review pp. 594-600 (19^); An Ext^^on of the 
Principle of the Diffraction Evolute and Some of Its Structural Detail , Phys/cal 
Review XX, pp,. 607^^12 (1922). . * ' 
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The discussion of Rayleigh scattering is based dn 'papery number ^.ind 247 of 

U>td Rayleigh, Scientific Papers , Dover Publications, Inc, New Yjofk, 1^4, 

The extension of geometrical optics to diffraction and to other "physical optics" 
phenomena is given by 

k * A 

- J.i. Keller, A Geometrical Theory of Diffraction, pp. 27-52, ['raj 
' Vol, 8,.Amer. Math, ^oc, Providence,%.I,, 1958, 

Intrpductory discussions to mathematical aspect^ of wave pheAomeji^ involving more ^T* 
than one scatterer arfe included in - 

^ . Twersky, MultiRle Scattering of Waves and Optical j>henoir i 



Symp, Appl, Math, 

I 

enawJpumal of the 



Optical Society of America 52, pi>, 145-171 (1962)^ 

in particular. Section 1,^ on scattering^ by two small obstacle 
discussion of one obstacle in Section 9 qf Chapter 3,' 



i.,could be read after the 



/ 



Chapter 2 

GRONA/TH, DECAY, AND COMPETITION 



Introduction, • ' . " ^ 
1 • ^ 

We f^yognizie^ttiat^^^ areas'of ioday's sciences have linder^one lo^g 

historicaLtievelopment to r^acli their predsnt systematic deductive stage, T 
j havebeeifi earlier stages of observation and classificat^ of superficially related 
phenomena, of measurement and the collection of data, of seeking detailed interrela- 

tions, of isolating essential facto of creating new concepts to link the more obscuryS 

L/ / > / . • ' ^' ' h 

interralatioris, of proposing ''laws of nature*" to summarize and generalize observa-j 

tions, of tefeting such laws to determine, the class of phenomena to whion they applj 

' of u(v6nti|^ more general more abstract laws, and of predicting less /3|^v;^ous 

f 'phenomena, Cycles'|Of observation, speculation, and Verification m^p^'our attei|;(pts 

to/systematize oi^r feiowledge of nature. We ^jq)loit yesterday's hciejiice td invesnf 

tjc/day»s(tools to dis(>over the science of tomorrow, ' v 

'/ / ' ■ ■ ' ' • ' -J. 

Mathematics plays its role oil every stage of the develppme^t of a science info a^ 
'deductive system — from the initial classification of related phenomena, to the se2p«oh^ 

/ ' * / r " ^ ■ - - - ^ . t ^ . 

for*tlie least number of fundamental principles on which they depend — from the estab- 
/ * * ' % ** 

lishnient of laws of nature which serve es;5entially 2§ axioms fqt a deductive systei^. . 

to the unfolding of the consequences implicit in such mathematical models of natural 

phenomena, Science deals ^vl^h phenomena — with observations, predictions, and, 

experiments on natare. Special areas of science requite special,equipment, and ' 



$p"ecial i]9.e^rement techniques. AJl areas.require mathematical thinking, 'mathem- 
atical 1 

d'^ 



matical tools, ^and mathematical Aodels. 

ll,the]sciences, frpm astroiiomy to zoology, ui^ mathematical models for com- 
plijfred'pheftotiiena observed in nature. To construct a model, ^ isqlate the'effects 
^that appear to l)e fundamental, ^d we define relevant variables, pRrameters, and 
functions. As suggested by oi^r observations and measurements, we seek appropriate 
equations for the dependence of, the functions on the essential jvariables; f^r complete- 
ness, we may have to introduce auxiliary^pnditions that specify, for exam^ile, the 
initiaa values pf the l^mctions androf the variables atjhe start pf a process. The 
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- -solutions of the equations subject to the auxiliary restrictions may then be compared ^ 
\viih additional measurements to determine their domain of applicability in nature. 

We make observations, \ve create models,' we make predictions;, we make more 
observations, more models, more,iJredictions, we day-dream and jump io -conclusions;^ 
we seek to verify our guesses, and keep the Very few that f>as^ the tests. By such 
means, by a mixture of measurements, mathe,matics, and mysticism we*seek to 
"understand" what is goings on around us. If we can predict and describe a process • 
and relate it to analogous processes that we know about, we are content — for a while. 

Despite the marked differences ih the classes of phenomena ^stud^ed in the various 
sciences, tliere are d number of processes that are fundamental: to all. The same 
mathematical structures, and consequently the same underlying concepts, are en-- 
countered again again in different contexts of nature. The equations are the san^e, 
mt the functions and variables represent different measureable quantities with dif- . 
ferent names that change from science to science" The stages and settings are very 
different, and the overall plots vary; but the subplots are routine, the actors go 
through the same motions, and* only the' names of the characters are ohangecT^ 

r ' ' " " ^/ ^ 

Let us consider a class of processes that is common to all the sciences^ proc- 
esses tl^at.involve such notions as '^growth," ''decay," and "competition." We can 
construol; mathematical models for such processes by using either vei5y general - 
abstract mathematical language,'^or by^j^ing the special language pf one of the special 
sciences, Jo show complete impartiality, we do neither. Instead we use the language 
for story telling. ^ 

We begin with a story -a story about stories. Then we show how the same story 
jcan be told again arfd again and again .... ' . • ' 

' • . . . - 

2. The Spread of a 'Story:. Model for.GroNvth . . , ^ h 

^ Once upon a time (to) I told a number (Nq) of friends a story about my good 
friend Al. Months later (time t) someone came up and asked, "Did you hear the one 
about Al?" Since I had started the whole business, I didnH have to listen.' Ini^tead I , 
asked myself, "H\dw many [N(t)] people have now heard the'story about Al^V 

How many people know the story about Al? Good stories spread, and this was a 
good one; the number of people that know it grbws with time. The number N(t) of • 
pebple w^ know it at timcyt should be proportional to the origSial ^uiiiber Nt) that 
were told the story at time tQ — to the Nq storytellers that couldn't keep a good thing 
to themselves. The older the story, thrfmore people know it. Therefore N(t) in- 
creases with the length of time t ^ to that the story has been circulating, as well as 
with the number of people available tp spread'Ht. If N(t) know the story at time t, 
how many N(t) know it at a slightly later time t? Since N(t) - N(t) must vanish jf ^ 



eilhfer 9(t)t.or t - t vanisKes, it is plausible to expect that the number, N(tX - N<t) of 
people who iearn«the story in the small time interval r - t is directly proportional, to 
both Nrt) a\id to the interval r - t. accept these ideas a's the initial aj^sjam^tions^^ 
and exp.ress them mathematically in the form ^ ^ ^- 

< V N(t) = N(t)o vf AK(t)[T - t] ,* ' Niftd) No ^ . , 



(1) 



where A is a positive constscnt — the gro\\^h coefficient . (Have we left out anything? 
Sure. We ^Halk about that later.) 

^ Accepting (1) a^s an adequate mo'deUfoj the'(Jhange in^N over a small time inter- 
^val still does not tell us how N(t)' is related to the initial value Nq, To determine 
this, we let tKe time interval approach zero and thus rei5laded(l) by a-differential 
equation, and then integrate over timQ to obtain N(t) «in terms of Nq. 

From (1) we have * • » ^ • • 

• i . ^ 



- (2) 



• , . ^^(^) y = AN(t) •■ 



limit as t approaches ^,'we obtain , 



(3) 




If we dis'count the fact that frrends come in integral packages (usually) an^o to the 



Equation (3) states that the instantaneoudi^Se of' change of N is proportional to N: 
this is th^ basic equation for growth. Later on, we will also consider the cas^ where 
A isneg?itive;^with A negative we havey|^>asic equation for decay. (If A is zero, 
then'N is a constant, and there is nothingHp talk about — neither for potential story- 
tellers nor for us.) 

In order to^express N(t) explici^y in termsT^^.N, we integrate 

.f-^-.'^^y.^ . ...... ^ , /. ti t 



fog mo - log N(to) p log =J^^) ;.. 



tj N(t)_ dt ' 



• and obtaS 
(4) 



Expresslng both sides of the* equation in terms of exponentials, we have for"l:i = t,-,;- 
(5) \< 



N(t) =N(t„)e^(**'-'V=^N„eA(*--5.) 



As a check, we differenti^ (5), and Verify directly th^ 



(6) 

in accord with (3), 



dN(t) 
dt 



ANoe"^' ''o|/ak^). 
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20 



For.cofJvenience m all that follows we take to = 0 as the initial time. Thus (5) 
becomes' , , , ^. 

(7) N(t) ='Noe'^' , No = N(0) , ' . . 

where t is the time that has elapsed since the start of the process. . / 

• From Equation (5),. we se^ that N «<> as t oo (i.e.^ n ^creases beyond any 
, bound a^t approaches infinity) which is not realistic for^ whert we know about story- 

tellihg (and other growth processes^ *Later on we consider a more realistic model. 

The present model is incomplete and should be restricted to moderately small time 
, interval t. " A • i 

•We have told a story about stories to get to (7). Now that we have (7), recog^r 
nize pat the result has other iiiterpretations and that the analysis has 6ther applica- 
tions.' Equation (7) provides an.elementary mddel for the growth of timber and 
vegetation, the growth'^f populations (people i bacteria), y{e growth of moflev^ttj banks 
(generous banks T^here they credit tlie interest fo the. capital instantaneously), the 
growth of a substance^in t^ie course of a ch^ical reacti^;i, and so on. ^ 

^ We can now ans\vex such question^ a^; ^' ' * . 

C ' . ; . ' 

If I tell 2 people the story at t = 0, and if the number' N(t) that i^iows it grows 
at a rate proportional to5v =^1 per day, then.how muc^^*N(t) at time yt ^ 7 days? 
The answer from (7) is 2e^ or approxiniately 2193; thus more than 2000 fteople know 
the story a week after I started to spread^it, 




If I deposit ^10 at 5% hiterest per year and theibank ad(^|^e^ttteresFtot^^^ 
inal amount continuously, tlien when wil^ it rekch $20? FoV^ . ^ r 

A = ^ it follows from (4) that log (f|) = (^)- t . ^Thus 

t = 20 1og2 = 20(0.693.-.) 13.9 years. * * ' . ' 

* - -» * 

In the following we refer to an equation, say (3), of this section by writing (2:3), 
etc. " ^ * * 

3. Radioactive Disintegration; Decay witl\'Time. ' ' • - 

The same considerations that led us to our simple niodel for growth apply^equally 
to the ^analogous model for dec&y (negatTve "growth)T We~take a'negatiye constant of 
proportionality in (2:3) to correspond id N(t) decreasing in time, and apply 

(1) . ' ^ .AN(t) C N(0) = No N(t) = Noe'"^*'. 

to' the problems r^ieactfve decay. Different radioactive substances disintegrate at 
different rates corresponding to different valiies of the decay .doefficient A. It is con- 
venient to exfiress thfe coefficient in terms of the half -life of fhe substance, the time 



it takea half of the initial amount of substance to disappear. (Why not thfe wl^ole"-life?) 
If T is the half-life,«the5n frorii (1) we have ^ ^ ' , 



' so that , * • '-^ 



(2) 



..i No ^ - 2 



-log£ ^ 10^2 ^^- 0.693 



Htilf of the material Nq will be leff'^t time r/one-quarter ^11 be left at time 2t, etc. 
When will it all be gone? We see from (1) that in order for N to approach 0, we 
require that t approach infinityj-we h^xe N for all A (all substances) as^^'- «, 
and thi^ is why the wbole^cUfe is a useless measure, 

LeFus consider a specific example. The half-li|p of radium is about 1600 years, 
and the corresponding decay coeffic.ient A is . ' . . ' * ^ ^ 

- ^ ^ 0.000433 per yei^r . ' . ' 

If we start with some given amount (Nq> and wait a hundred years, we get 
log^ s^^-0,^^, and consequently N « O.DSSNq is the amount left. "Thus only 4.2% 
disappears in one hundred years. 

The basis for applyii^ (1) to radioactivity is statistical, i.e., it holds in the sense 
of an average. Although the physical process governed by probability,. and we can- 
not tell when any one atom will disintegrate* it is quite useful to determine the mean' 

life -time per atom. We start with Nq atoms at t = 0 and end up with 0 atoms as. 

- 

t approaches infinity, and we are interested in the average length Of time that an 
atom exists. . / / ' 

If nj^ atoms disappear at tixpe tj , and n2 atoms at time t2 , etdL, where^ 

tj < t2 < • • • tk 

jjid \ m = Nq, then the mean life -time of an atom is the average value _ 



"i=i 



If the total nunlber of atoms present in tfee interval (ti-^ ,ti) is Ni , shd Nj+j isjhe 
numberVre sent in the interval (t^ ,ti+j), then. 



and the mean becomes 



m ^ Ni - Ni+i , 
4 



k 

^0 



^ ^ t, (N, - N,,,). ' ' 



i=l 



Treating the relation between N and t as though it \('ere continuous (essentially as 
for (-2:3)), we define the mean life -time as «. 



(3) . T = -j^- / tdN= ^ /^"tdN. 

We now seek the relation^b'^veen the mean" life -time T per atom and tHe^S'cay 
'coefficient A. ' . . ^ 

From (1), we have * ' % 

7 

, (4) dN(t) =^ -AN(t)dt = -NoAe-^*dt . - > , 

^Using (4); we change the integration variable in (3) to t, and the limits Tip and J) to 
the corresponding values t = 0 and ooi ^ 

(5) ^ . * T ='A r te-\dt. ' ^ ^^ . 

To evaluate \hevintegral in (5), we regard it as the form /udv (witl\^ = t aad 
dv = e"^*dt) and integrate by parts:, 



"f-'te-^'dt ^ I'^.^T = i- 
•'O L -A ^ -A^ Jo A2 



. Substituting in (5), we obtain 



"ilius the mean life -rtime.is the reciprocal of the decay Coefficient. ' 

^We have a model for simple radioactive .decay. What is left af^br an atom _ 
disiiltegrates? Many things, including Vdaughter*' atoms which can also disintegrate. 
Later on ^ talk about some mother-daughter relations^ and discuss the decay of . 
both pogulatlons. \ 

The simple decay model we have been considering also describes essential 
features of many other phenomen^^As another example within the same mathemati- 
cal structure, we need chaj^ge only the names of the characters in ordelFfor the^^ 
results to" apply to the molecules of air in the room. Suppose that % is the total 
number of mplecules present, N(t) is the number that have not had collisions by time 

t, and that the^ma^iJlSe between collisions i^T. »Theh ^ , the probability tHfttj 



any one molecule goes foretime t without a collision, follows directly from (1) and 
(6): 



(7)- 



m - 

■No. . 



-t/T 



/ ■ 



1 



Same elation. as previously, but the characters are now playing different physical. 

roles, and, of course, the overall plot is qult6 different. Were \ve to continue the 

present story we would require mucli additioYial structure. ,lf the mean velocity of the 

molecules is V, ttien L =^Tv is call^ tl>e mean free\path— the average distance a 

^ - . • t 
molecule travels between collision^ — but this path leads to statistical mechanics 

(which bases the physical properties of^ matter on the motions of molecules) and would 

jtake us too far afield. / • 

' " ' ' / ^ . ^ ' ' 

Electrical circuits offer several, examples of decay processes. The fundamental 

roles are played by charge q , curr'^t I = i voltage or electromotive force V. 

Simple circuit components (capacitors. er condensers, resistors, tnductors»oi' coils) . 

-''V 

^jare c^aracterized^by constants: capacity C = ^ , Tesistancfe R - y , and inductance 

/dl ■ / ' ^ ' * ^ 

L = V/gf / If a condenser at voljtage V discharges Across a resistor it initiates a 

process described by 



(^) 

and consequently 
(9) 



dV -V. 
dt 5r ' 



VOt) = Voe 



-t/CR _ 



qo 



"f/CR 



where 1(J is the initial voltstge on the condenser and qo the initial charge. A current 
I flowing through a circuit coni^3ting of^a^oil and a resistor satisfies 



dl _ 

dt " "L^ » 



(10) 

• and consequently, 

- . aS: 

T t 

(11) . '7 I = loe 

where Ij is the initial value of the current. 



(, 



Rt/L 



'4. v The Colors ofjhelSky: Attenuation with*l3istance» . * . 

The model of Secjtion 3 for a process thdct'decayg with time can a^so be'Sppuld 
fon processes that attenuate (weaken) with distance. Ifet us relabel tl;ie yariablecS 
(3:2) and call itj^distance," and write it as "x*' to keep the neW^le in view. We ' 
thus have / ' ' -* * ' x ^ ^ ' ; 
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(1). ,^11^ = -AN(x) , N(0)'=No;' N(x) = Noe"^* . 

which we call the attenuation equation . 

The- attenuation of the earth's atmosphere nvith altitude is described appi-oyimately 
% (1) ^th- x-.as>th^ height above the earth's. Surface, N(x)' the number of molecules 
per unit volume (the density) of atmosphere at height x.'and Nj the density at.gromd 
level. There are different kinds of molecules in our atmosphere Mth different 
• masses m, so that we sTiould introduce m as aphj^ical parameter and" write .'.'"^ 



(2) 



N(m,x) = No(m)e-^H<^ 



.for each constituent! If one introduces i^ore structure intp the model, it turns out 
that A = ma, where a is independent of m (but depends on temperatui'e, and the 
acceleration of gravity). . . " ' ' • 

As an alternative setting for (1), visualize a narrow beam of particles incident pn 
the face of a medium of more massive particle's as in Figiire i.. 
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FIGURE 1 . 



There ar^j particles p^er unl^t volume oj the incident beam anjl^othing.happens to- 
thpm until they encounter the medium that starts, at x =' 0 . Then as' the beam p^^e- 
trates, itg. lighter partifeles hit the heavier ones of the medium and Jo off in other 
divections: particles of the incident beam^are lost to other directions by scattering . 
(Visualize a column'of/jhildren frying to charge through a crowd of mijling adults.) 
Thus N(x) , the density of particles in the incident beam at a distance, x within the 
medium is less lhan No"; this attenuation is governed by the scattering coefficient per 
unit length A . (For other processeBTaenfifytloi| We beam may attenuate because 
its' particles combine with the heavier -ones bf-the medium.) 



The principal characters in the above are particles —lopse characters that cdn 
stand for electrons, proton^, children, etc. If we now relabel N as energy density 
per unit volume or intensity , then the same plot also holds-for light-rays, x-rays, 
Y-rays, and all" other kinds of waves meeting appropriate ob,stacles. For tbe 
particles, we mentioned one physical parameter, tWe mass m , and spoke^of particles 
being scattered by more massive ones. 'FofwavWs, the appropriate physical {Pa- 
rameter is the wavelength \ — e,g'y the distance between* adjacent crests of equally 
spaced ripples on'a lake. The lowest wavelengths associated*with visible light give 
the sensory impression callgd'*'red,'* and they are about twice as long as the shortest 
of the Wave^ngths associated with visible ''blue.'' From blue light to ultra-violet to , 
x-rays to -y-^rays we go to shorter and shorter electromagnetic wavelengths, from r^d 
li^ht to infrared_to microwaves to radio waves we ^o.to longe,i»-electromagnetic wave- 
lengths. We cau^also talk about sound waves, water waves/and even the waves of 
"probability amplitude*' associated with electrons, neutronL and other fundamental 
particles. ^ , - 

With N for intensity, (1) in terms of an appropriate A describes^ the attenuation 
■ of a beam of sunlight penetrating a cloud or a layer of fog, etc. We can use (1) to 
^determine the thic'Bness lead shiefe's to be used<£with dental x-ray equipment or 

with a nuclear reactor to reduce stray ra^tion to a tolerable value. We could talk . 

in greater detail on aqy of the above, but iitst^d let us talk about something more 

colorful. ' * • y \^ 

Let us consider Rayleigh^s theory for the color of the sky. ,The jBSsential feature 
of sunlight is that it is^ade up of light of different colors from red tP blue (the 
visible spectrum) \^^h.associaj;ed wavelengths Xr to \^ such that (approximately) 

(3) I * ■ Ar = 2\i, . : ^ 

The wai^length \ of an intermediate colc/t (orange, yellow, green) satisfies 
Xj. >\ > Xb-*' Rayleigh showed that when a beam of Jight of wavelength \ Is 
. scattered by the molecules of the earth^s atmosphere (mainly nitrogen and oxygen),. . 
the intensity N(\,x) along the beam is governed approximately by (1) with 

(4) " ^ ^ •A(x)=^ ; ^ 

where C is independent of \. (In the chapter on optics and waves we discuss this 
more fully.) 

. From (4)'knd (3) we have- » 

A(Xb) ^ ' ' 



(5) • = ^ = = 16 



and consequently 



^ ' / " • / / . ' 

Thus the blue componefif of wljite light is 16 times more strongly attenuated than'thQ 

^6^1 A beam of \V^ite sunlight, reddens with penetration into the earth's atmosphere 
. be^SQ it loses i1;s blue cbmponent mor^ rapij^y.than its r^d. The blue that is lost 
Jrom thQ sunbeams by scattering gives the sky its blue color in directions away from 
the sun. The direct beamsrfj*pm th§ overhead sun are still relatively white because . 
they have not lost that much blue. The reddening of the diregt^ beams is best seen 
when the sun^s'low on the horizon and its rays traverse maximum distance through 
the scattering atmosphere; the cloUds in the path of these rays^^are bathed in red. " 
Such colored effects and other scattering phenomenaj^arising ftom water drops, dust 
p^icles, and other impurities in the atmospherV ar^&or.e fully discussed in poetry 
courses. * , v ^ / 

5. Mother-Daughter Relations; Birth and Decar . 

^ As mentioned previously, when a radioactive ^tom (the mother atom) disintegrates 
it may give rise to a daughter atom which can also disjjitegrate. Let us now consider 
such mother -daughter relations. ^ ^ - ^ 

Suppose we have Nj mother atoriis which decay ^t the rate 

(1) ' . ■dr--AiN,, Ni(0) = N,o. 

The rate at \^hich the mothers decay equals the ,rate at which the daughters are 
created. But the dau^hteil^ also decay on their own; If Nj mothers with birth coeffi- 
cient A^^ve rise 'to d&ughter^ with decay coefficient Ag , the rate 6| change of 
the number of daughters is given by \ 

' . . dNo ^ ® ) ' • - \ 

(2) dr^= ^i^ir ^2N^' .,.N2(o)-p. I " . 

Equations (1) and (2) form a pair of simuTtanlo^liMquations for determining Nj and 

us first consider a limiting case sueh'ihat the mothers decay very slowly 
frompared to the daughters,^:e., the decay coefficient ofthe first is much smaller 

-^than that of the second, Aj << Ag . This corresponds, fo^^example, to the behavior 

- 'f 1 

for the pair radium-radon. For radium mothers, the half -life is approximately J.600 

^ # - -m- 

1 

ye^rs: ^Aj = per y^ar. The radon daughters have a half -life of about 4 days: 

A « 1 360 days ^ 1 ^ 90 - - ' 

^ 4 days""^ year ^ 4 days ~ year " ^ 

Thus A2 « 144,000Aj ,.and we may take the number Nj of mothers constant for the 
purpose of obtaining a /irst approximation of Ngjft) , • " ' 



Thus we tegroup the terms in (2), replace t by tj * and integrate from tj Kfi. 
to ti'= t, ^ ^ [ ^ . 

subject to the^ approxiigi|ition that Nj is constant. We obtain 
^\ * logK(AiNi ^- A2N2) =.M2t, ^ 

and; - ' . ^ . f • 



K(Aikj - A2N2) = e^. 2' . 



where the integration constant K is to be determined from the initial conditions at 
t = a. There are no daughters at t =*0: ^ ' • v ^ 



CoiiLsectuehtly, the number^ of daughters at time t^is given approximately by 

N2 = (^)Ni(l-e-^2t); 

- ' > A / A.N 

If is vdl^ large, then N2 approximates ^ ^ , i.e., the number of daughter atoms 

approaches a fixed fraction of the relatively inert mother substance* (This is called 
long term or secular ectuilibrium.)^Wh'^does this mean? It corresponds, for . 
example, to the case Where N2 is a gas (such as radon) in a closed container, and a 
situation where just as much^ N^is created (from Nj) as is destroyed by radioactive 

f. * ' ' dNo . 

ly, The/birth rate of N2 equals its death rate, so that -g^ is zero; our result as 



t approaches infinity in (3) is^us the same as that obtained by equat^g (2) to zero. 
Equilibrium corresponds- to i ♦ ^ 



(4) ' - - §- = 0;, N,=|||)^I^," 



Betore continuing the development let us feedback the above to our earlier ^dis- 
cussion of electrical circuits (3:8). Equation (2) with i^j held constant is also the 
mathematical relation for tlieiburrent I in an electrical circuit having resistance R, 
* inductance L7and^extemal Yoltage V: r 



dl.^ V R 



V 
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If V is constant, and if I = 0 at t = 0 , then the solution o! (5) for the growth, of 
current in time is obtained by. inspection of (2) 'an(U3): 

(6) . I - e-^^A). 

'If V = 0 and I - Iq at t = 0, then (5) reduces to J[3;iq). 

N6w back to ihothers arid daugljters, and let the* mothers also decay in (2), To 

take this int^account, we substitute, the solution of (1), i,e.^Ni(i) = NjqC ^ into (2) 

and obtain , ' ^ -.^^^0^ ^ T * " 

' , \^ i 

We solve for^Ng by using e"^^* as the "integrating factor." We note that 

* • . f ^ . 

d/„ V\ dNg At ^ At /^j^ V A t 

To exploit this, we multiply both sides of (7) by e'^^* obtain 

d / ' A t\ A t- A t ■ ' . 

and then replace t by tj and integrate from t, = 0 to 1, = t : 

" V'-t = A,...j;,v.-v.,.,, W^;.--. 

Since N2(0) = 0 , weHhus have ^ 



A. -A. 



which reduces to (3) if. Aj « Ag , and A^t 0 , In distinction to the approximation 
(3), the present complete form Ng of (8) vanishes both for t - 0 and t - «; conse^ 
quently N2 must have a maximum at a specific value of t , 

, If. we differentiate (8) with respect to t we^ obtain ^ 

. I^-'AT^A '^ +A2e 
^> This vanishes, and N2 has a maximum, when - * ^. 

From the logarithmic form, we obtain 

■ /im -■*■ - . - logAf - logA; 

^'"^ . . : ' = A, -.A5 ' . * 

as-the time when the number of daughters" is largest. The maximum value of ^ 

^augWerk is N2 - Njoe""^"* = IT^lNioe "^'^ = /p-jNi as in (4), 



6. Neutat Processes: Stimulation and Decay , 

" One methodpf studying the excitation and propagaticm of nerve impulses is based 

_,on stimulating a ner\"e fiber electrically. Characteristic measurable effects ai:e pro- 
--ddced in .response to an electrical stimulus V (the voltage associated with a direct 

current, the discharge of a ijondenser, o,r an alternating current), provided that V is* 
_greater than-J^jtbreshold yalue Ve —the minimum value of V that is just sufficient to 

cause the effects, A simple mpd^l (introduced"by Blair on the basis of experimental 
~;data) describes thp bnset of the effects in terms of a local latency (ali^o called the 

"local excitatory functiog?^^^(t) such that ' , ^ 

- - - =KV(t) -.AN(t), 

Where K is the growth of the latency per second per unit stimulus,^nd A is its decay 
, coefficient. Thus the growth of N increases with the magnitude of the s^mulus and 
decreases with N, (The function N may represent the difference between the con- 
centration of an exciting ion at an electrode while V is applied and its concentration 
for V = 0,) If N(t) reaches (or exceeds) ajhreshold valu^ Ne, then the nerve 
becomes excited (and a characteristic wave of physical -chemical changes \vith an 
associated electric potential propagates'' along the fiber), 

TJje simplest/application of (1) is to the situation 

(2) N(0) = cr, ^ V(t) = V = conslant , 



which corresponds to the application of a constant stimulus at time t = 0, By 
inspection of (5:2) and (5:3)'^the solution of (l)^and (2) is , 



(3> 



KV 

Thus as t we see that N approaches its largest value Nmax = . Conse- 
quently "excitation will occur if , ^ * 

KV 

\ ' 

' or equivalently if the stimulating voltage satisfies " ' ^ ~ . • 

', (5) •■ . • " -T^^='V., ' . 

where V© is the threshold stlniilus mentioned previously. (The value Is known as 
the rheobase , the threshold or liminal value of the constant voltage necessary for 
excitation,) ^ . - ' 

Assuming 1;hat V > Ve (so that excitation must occur), therfthe nerve becomes 

excited at the time t© corresponding to N of (3) reaching the th'reshold value 
^ , ' / , 

(6) N, =-^^l-e 
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^^^or, equivalently, 

• .IX 

(7) . . . , 



which is the latent period that dlapseis between the (istablishment of the constant 
stimulus and the release of excitationV 'If V < \i, iio le^^l v^ue of tg exist^^ If 
V.= Vc*, then te ^ however, tijis is\an inconveniisnt length of time for measure- 



ment pui;^ses^ A more convenient measure is the 
y = 2Ve: 



(8) 



_ \log 2 \ 0.^93 



value of te corresponding to 



This is known as the c hron axie f— tne latent tiir e l:efor0 excitation for the case of a. 

\ \ I • I 

stimulating voltage equal to twice the threshold v^hie. \ ^ 

What hav^ we bgen doing in the above? 4-essen ially changing the names of the 
characters introduced for radioactive decay and shewing how a few simple ideas link 
parts of biology, physics, and chemistry. LetVis no^ generalize the mathematical 
^ development to nonconstant values of V in (1).^ 

If V is a function of time,, we solve (1) in teVms of V = V(t) by proceeding ^ 
essentially as for (5:7): we transpose AN to the yeft side, multiply by e^* to obtain 

.(^ + ANje^^rp ^(Ne^M =.KV(t)e^* . integrate th\ le^t side explicitly, and regroup 

to obtain • " 



(9) 

If N(0) 



0 = 0„ and V i^s -a constant, then (9) reduces'to (3). 

If we stimulate the process by discharging a condenser of initial charge q, 
capacity C ,^and resistance R, then as in (3:9), 



(10) 



V(t) 



Substituting in (9) and integrating, we obtain for Nq = 0,* 



(11) 



KqR 



(CRA - 1) 



which is simply, (5: 8f) witjjgifferent labels. Thus the excitation function N has a 
maxlmiun when * - o - . 



5^ . f 



(12) 



•"t = 



CR 



(1 - CRA)""^ \CRA. 



log 



)■ 



I! 



If the maximum value of N is precisely the threshold v^lue, thefi t of (12) is the 
corresponding latent time from onset of stimulus V to release of a wave of activity 
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injjife nerve. The corresponding initial voltage V(0) = ^Ts the threshQld, initial 

voltage of the condenser for excitation to occur. If the maximumr' does not ^qual the- 

threshold, then we relate the condenser^s chalracteristics to the threshold by - . 

AN ' " 
equating <11) to and using Ve - to eliminate K, ^ — 

If we stimulate the process by a sinusoidal alternating current, then the applied 
voltage is, ^ '* . 



(13) - , 1 v(t) = Vo sin-tjjtv — 

- '^^^^ 

where Vq is the constant amplitude. Substituting in t9)^fop>J^0) = 0 we have 

(14) , N(t) = e-^*KVo ^sincotidti. ' / 



Ti'ansposing th^* integral fj'om'^he right-hand side to the left and.dividing through by 



1 + —5 we obtain 



' k^^' sin U dt = -^eAtcosa;t -f Ae^^ sin cot 
f > o)^ + A^ V 

'T^nsequently tlie solution'of (14) is \ 



To evaluate the integral in (14), we regard it las ±he-^tor-m JiidV (with u =^^* and 
dv = sin cotdt) and integrate by parts: -\]\ 



This.idoesn't look as if it will get us anywhere, but let us handle fhe new Integra by 
again integi^atii^^ (vyith dv* = cos aitdt)': ' ^ " . .. 

i-' fe* sin ^X. dt = -^^*cos.a;t .^AeAt _ A? J^At a;t dt. ' ' " - 



1 



KV ' \ 

(15) N(t) = ^2 ^ (^ sm cot - CO cos cot + coe^^*.) , ^ - , ' ; 

The eKpohential term of (15) is significant only for small values of t< As t 
increases, e"^* becomes negligible and (15) reduces to 

(16) • N(t) a ^2^J\2 sin cot - CO cos cot), 

V This periodic approximation has equally spaced extrema in time, which occur when*' 

(17) ' ^ ^ ^|Voco ^^^^ cot 4- CO sin cot) = 0 , t = i tari^^fl^ , ' ^ * 

dt CO + A^ - CO . \ CO/ ' 

Substituting these values oM int<>Tl6y, w6 find that ^he^maxima of N equal 

(18) - N„ 



KVo 



max - (0,2 + A2)'/2.* 

O ' - 1 ' 25 DO 

ERIC / ^2 



If we equate Nm»x ta tb^ threshold value Ne , then Vq of (18) corresponds to ^he 
threshold valu^df the=^mplitude*of the sinusoidal stimulus, say ; however, the 
♦ result oc + has only very limited validity in nature, 

7, A More Realistic Story About gtories: Bounded Growth. 

' ^ Let us return to our model for the spreading of stories (or of diseased, or of ink 
blots), and introduce mor^ structure. Previously we assumed that the rate of change 
of the number that knewlhe story at time t was proportional only to the number 
itself: ; ^ • . . ^ ^ . 



(1) • - = '^N(t) , -^(0)' = No 



This is all right as far as it goes„but it ignores the fact that tKete is an upper bound 
(say Nm) on the number available^ to hear it: there is^ a fiivite number of people on 
earth, same don't talk yourjanguage, some don't talk at all, and most never listen. 
Furthermore, although we may tell tl^e same person the same story a dozen times, 
, each listener should be' counted only once. \ 

In view of these con side ration^, we replace (1) by 
(2) f = AN^^^1q^^ N(0)^No,> 

where N = N{t) know the story at time t and are available to spread it, and Nm - N'- 
do-notknow the story, have good hearing, and are enthusiastic' listeners and potential 

(N - N> ' ^ 
gossips. The factor ^ is the ignorant fraction of the available population. 

Dividing both sides of (2) by N^, we introduce f = as the fraction that know the , 

Story, and work with' . . - , ^ 

(3) 



|=Af(l-f), f(0)=fo=|-^ 
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Our original model (1) gave N as t '^cc." What, does the present model give? 
We exi)ect that f 1 as t i.e., -that eventually everyone knows the stoiy. (Even 
this model is far from complete, but at least this kind of resutt is acceptable,) t'rom 

df * • * y - ' . 

(3) we see that ^ 0 as f X, i.e,*, that f stops changing wh^V eVeryott<5 knows the 

stoiT; from the discussion for (5:4) we may surmise that ^ - 0 as- 1 - <^o, but let 
us solve (3) and see the details, * ^ * '5 

From (3), we write ij^^^Zi^ = J Adtj, where fj and ti are dummy variables. 
• • • ' ' 

We decompose the integrand into "partial fractions," Since 
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r f(i -f) " f 1 

we obtain * 

r[A.-^\, 

Thus solving for f , w§ have- 
foe'^' ' •• 



4. 



f f 

= log TT? - log TTTT = 

^0 



(4) * . ^ f = 



1 + foCe'^* r 1) \ 



If t is small, then the denominator approximates unity and f = fge^* , in accord 

f 

with the simplified model (1), On the other hand if t is large, we re\yrite (4) ^s 

fp ' • " ' \ 

' ^ ~ fo + (l-fo)e-^^ . . * " > 

from which we see that f as t «, 

The above model indicates some of the essentials but it is still incomplete. 
However, it is good enough to show that although you may ^ill not have heard the 
stery about Al (see Appendix), you should by now have heard about the Calculus — or ^ 
at least about Helen of Troy^ 

8, Population Problems: Growth and. Competition . 

A more general equation which inclines (7:1) and (7,: 3) as special cases is 

(1) . ^* ^ = AN J BN^, N(0) = N^, ^ ^ V\V ^' ^ 

This is called the logi^^ics equation . We still call A the growth coefficient , and we 

may caU B fiie braking coefficient because the term -BN^ slows the growth. The 

dN * * 

equation of unchecked growth, "gf = , permits N «> as t « , but Equation (1) ^ 

does not. What boimd does (1) impose bn N? As in Section 7, we see that = 0 

when^ = BN: the correspondijig value N = must.be the equilibrium value . _j 

which /n approaches as t <«. ,x * . » 

At one time, essentially ^ ~ AN ^feuged as a model for the growth of popula- 
tions of different countries, and this led tQ dire predictions as to the fate of mankin(i 
" (lasVof Malthas). Then, essentially (1) was introduced (by Verhulst) and -appropriate 
\ ^ ^ ' * ' ' ' • • ' ^'^^ 
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A^s and B's for various countries were obtained from their earlier census records; 
the projected growth curves were remarkably accurate (at least for adl countries 
except Verhulst^s^^s^elgium). The bu^Jdup qf population gro^vth arising from A has 
been interpreted as due to cooperation between people, and the slow-dowa associated 
\vlth B ^as due to competition between pec^ple Jpr limited resoui:ces. The ccAnpetition 
interpretation is qiiite plausible: if p is the probability that a person, wants a partic- 
ular thtng, then p^ is the probability that two persons want it simultaneously; if ' 

there are N persons, then there are ^^g" competing pairs, and the total proba- 
bility of competition 'is P ^(^^1) . incorporating the linear term into the g^rpwth ^ 

term of the differential equation, we may take BN^ as a plausible measure of the 

Jmultaneous desire or of the competitivejurge. However, the reason for regarding 
^ as, a measure of codperation is not clear. X probability interpretation similar 
to that for BN^. indicates that AN corresponds to N persons acting quite independ-. 
ently of each other, this miay well be ag* close to cooperation as one pan expect from^ 
a group, and SMSG authors have therefore taken this as the guidii^ principle for 
preparing their textbooks. \ * • 

Let us solve (1) by the, same procedure \ye used for (7:3), The steps are essen- 
tially tlie same,, and we get 

and consequently ■ 

Noe"^* 



N-tB 



If t j[s Ismail^ then tl)e denominator is approxipiately 1 + NgtB^ « e , and (2) 
reduces to 



(3) 



t(A-BN^^) t(A-D) 
^ Noe ^ =Noe 



where D = BNq is introduced as an abbreviation. Thus for small t , the result has 
the same form as for the simple model in terms of the grov^rth^coefficient A - D . ' 
On the other hand, if"t «^ then the limit of (2) is A : ' - 

/ » , 

i- A * 

in accord with the remark after (l),that N - g represents the long term equilibrium 
value/ A 



9. Quizzes and Nonsense: Forgetting an^i Leaminff . 

' The previous sections also provide a simple model for forgetting and learning, 
.at lea^t of unconnected chains of nonsense syllables invented. by psytfhologists for 

4est purposes. Thus (as proposed by Yon FoersterJ we consider 

0- \ , ' J 

(1) ' ' ' ^ = ^N(No - N) , " N(0) = No, ^ ' 

where Nq ,is the initial number of items memorized (dates, telephone numbers, un- 
<;onnected theorems, etc.), A is a forgetting coefficient , and BNq is^a memorization 
coefficient . The notion behind (1) is that your head is originally filled with Nq 
''carriers^' of information;' some (AN) carriers just lose their information forever; 
3ome (BN) lose informatlpn in the sense that they pass information on to the empty 
N0 - N carriers. 

Integrating (1) l?y partial fractions, or'by comparison with (7:1) and (7:2) (the 
prespnt (1) is (7:2) with a new growth coefficient BNq - A), we write the .solution of 
*(1) as ' . ' ^ 



D - Ae-(^-A)t ; ^-^oB. ^ . ^ ^ 

.The limit of for t ~ «, the re membrance R (as defined by Von Foerster) de- 
spends critically on the magnitude of ^ , If D > A , then ^ , 

On the other hand, if D < A , then 

Thus the remembrance of things past is zero for ^ :^ 1 , and then increases towards 
unity as ^ increases from unity, 



^10, Chemical Reactions: Multicomponent Processes , 

Suppose we have a chemical substance with initial concentration C"'(gram- 
molecules per unit volume) which is reapting in time with something unimportant and 

./.plentiful to fom another substance with concentration N(t) , The rate of cbange of 
,N is proportional to the concentration of thp original substance 

^ C - N(t) : • ^ ' 

, ,(1) ^ = A(C -N), N(0) = 0, 



at time t , i.e., to v 



where N is the concentration of the new substance and A is called the Reaction rate. 
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Equation (1), which iff known as the law of mags actioa , is essentially , the special 
case of (5fi) for Aj mucir^ smaller than ; by inspection of (5:3) the solution (i.e., 
the concent rat ioir ofAhe §o lotion) ! 

(2) . : ^ ' N = c(i-e"^*). ; '\ .:^,>r'* 

Equivalently, Ecj^tiort (l), is a shifted versi(?n of the simplest decay equation (3: 1); 
sub,stituting M = C - N in (1), we obtain , ^ ' 

^"^'-^^'^ M(0) - C - N(oft'cV ' ^ / 

which leads directly to (2) for ,N = C - M . 

, From (2), we see that if^t = 0, then N = 0 ; if t «, then N , so that. all of 
the original substance eventually reacts. We may isolate A in the form 

which used in chemistry courses to determine A by measuring C, N, and t . ^ 



• -In a bjmolecular reaction , we have two different substances with initial concen'- 
trations Cj and C2 which react at a rate determined by A ^to prbduoe a third sub- 
stance^whose conceatration is N(t): ^ 

^) ^ ^ ^ A(Ci - N)(C2 ^*N), N(0) = 0. * 

This is just another variation of the logistics equation (8:1). The solution can h^. 
obtained from the previous ones, or directly: ^ . / 

r dN 1 r[_l 1 1 ... 1 , c/n ^ 

J (C, - N)(C2 - N) ;,C7Tc^ J [pPTn ^CT^nJ^N log ^ At + K, 

t ' ■ " • ^ 

where < , , . ^ ] ! ' 

( 

• . ' _ log C1-/C2 

follows from, N = 0 at t = 0. Thus , ^ ' ' • ' 

' * A ^ . l' • 1 <^2(Ci - N) 

■ /-•t(C;-C2)l°gc,(C2 -N).' . . 

_ . _ , ^ - ~ . f 

and ' ' " ' ' 



• -,r . /' 1 - e 

16) • . - . N = - 



TCi-C2)At . 



1 , (C,-r )At 
1 - .(C\/C2)e' ^ ^V. 



The. case C, = C2 = C may be obtained from the limit of (6) as C^ approaches', 
Cj. Equivalently, we start with 



130 ' ^''^ 
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The constant equals K = -yq , >and tliereifore;* - V, ' ' • -f 

^ ■'- * i- V* ■ ' - • ° ' J 

' • \ ^ ^ C(C> N) ' ^" ^. - 1 + CA • ' ^ ^ ' 

The equation for opposing iflSmdlecular and bimolecular reactions ha$ the forni 

^(9) , A(C.- N) - BN^, N(0) = 0. ^ * ' ; 

We do not discus^ this case but merely reduce it to a previous form. Thus we 
-iuti*oduce 



■ ■ ,■ V . .. ■ •. 

/ ( V in.order-to rewrite (9) as . , _ 

.V. :)• , •1' ' ^ 

^ (11) ^ ^ = -B(D, - N)(02 - N). . • 

*^ \ ^ . ' 

We now have the form (4) with the previous A, Cj, C2 replaced by -B, Dj, P2, and 
the CO rrespondin^e suits may be v^itten do\9h by inspection. 

We could go'pri to higheix-order reactions of the fcfrin — j;>; . ^ ^ 
(12)- • . , ,^ = (C, - N)(C2 - N)(C3 - N) .... 

T>ut we must finish the story. ; 



/ ' ' 1L> Sociology; The-E?nd. 

Now we could rehash everything. We could change the names of the characters 
in the previous-equations and tajk about profound sociological problems. Instead we 
introduce, a more general model for the growth ot populations, one which inpludes , 
practically all of our previous equations as special cases, and^'carcely talk at all. 

In Section 8, the growth of a population of N individuals was described by 

rJ (1) ^ ^ AN BN2,Z, W).=.N(>, . 

^ . -w , ^ . r ' ^ ^ ' , — 

whSre A is the growth coefficient^ and' B is the biding coefficient. Let u^pow 

^ . ' t \ • \ 

Introduce more structure. We may \vrite k ^ ot - ^ ^ where a. is the birth\c5*effi- 

cient (the birth rate per individual) and.wherq, ^ is one of two death coefficients. 

' o •*■ '.' . . 31 38 ^ 
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^"^TTwe assume that the population is confihecTf^tn' area S , then we may write the 
other death coefficient as k.. i.e., the death rate ^p^jUn^vidual increases as S 

^ dex?reases or as N increases (no room to Jive). Thus the total death rate is 

(CN\ C - ' 

+ -g-j N . Using y = g instead of B (merely for esthetic reasons) we rewrite 

_ (l)as ' ' ' ^ 

.-.■(2) ' . _ ^ = aN (^,+ yN)I^. 

* * • > . / ^ 

A more general mo^el (considered by Rashevsky) is that for the .growth of a 
population consisting of t\vo types of individuals with different birth and death ' ^ 
. characteristics-. The total population is j ' ^ \ "V^ , ^ 

and N, and N2 are specifiedTby the simultaneous equations 

dN, ' r T ' • 

-gt" = OiiN, + aaNi - \p^ +-y,(N, + Nj) N, , . . 3 > 

-g^- = a2-iN, + a22N2- +-'>u(Ni + N2)JN2, " 

Avhere the a's, p's and y»s are all constants. The tSlrms proportional to q repre- 
sent the contribut^ons-of the two groups to the birth rates; the death rate6 that depend 
on yi (with i = 1 or 2). depend not only on Ni , but also- on the total population 
Nj + N2 = N . The system of Equations (4) generalizes practically all the other 
equations considered previously in this chapter, , * f- 

We do nothing Vith (4), but as a set of exercises you could obtain all the pre- 
vious equations that we considered and that follow from (4) under suitable restrictions. 
• Talk about a:ctiVe individuals and passive individuals; talk about active and passive 
disobedience; talk about social ag|pegates, freedom^ crime, war, propaganda, etc.^ ^ 
Write a book about it and caH it War an^Peace. 

What have tried to illustrate with this chapter? As you apply mathematics to 
*^e various sciences, you soon discover that at a fundamental level ther6 appear to be 
" only a few different kinds of processes going on: the same mathematical structure . 
arises again^ and ^gain in different contexts of nature. We have seen that the equa- 
tions are the same, and only the names of the functions and variables change from 
science to science. The stages and settings are very differentTand the overall plots 
vary; but the^subpllots are routine, the actors go througji thfe same motions, and only^ 
the namQs of tbe characters are changed. ... . -r 
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The Director of SMSG has rjiled that the story about Al may be disseminated only 
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\ ' * Chapter ^3' 

" GEOMETRICAL OPTICS AND VVAVES 

1, Introduction . 

There can be wonder and excitement in following a thread of mathematics through 
several sciences and recognizing'their kinship in the concepts they share. There can 
be deeper wonder and greater excitement in following a single scientific thread through 
the cunlulative concepts that trace its evolution. ' 

l£X US start from' "laws of nature" that were isolated only 'after long years of 
observation, speculation, and verification, and show how various methods of calculus 
supplement each other in revealing the consequeuces implicit in these laws. Thle 
particular laws are limited in their^omains of .applicability in nature, and correspond 
to suitably restricted classes of observations. ^ we progress along the scientific 
thread, we' trace part of the development from the early very special laws of geomet - 
rical optics to the modem very general laws that constitute the mathematical model 
for wave theory. We use optitial terminology (light rays, mirrors, etc.) but the |con- 
cepts^e consider are^basic to light, sound, and all wave physics. 

• - We start with^^ery restrictive laws, weaken them, and thereby enlarge the domain 
of the subject. The^ procedure we follow is "quasi-axioma(ic," The "laws" or "axioms^* 
we list cont^ undefined terms well as implicit restrictions on both the observa- 
tional' and computational procedures that are associated with^them, (The'aSamns of a . 
matheiiiaUcal system fHso contain undefined terms, and are also usimliy Stated ' 
implicit context, but the situation is ixiOre obvious in J3ci6nce.) ^ ^ 

We follow a thj^ead that suggests the heuristic search for first principles on the ' 
, baisis of limited initial data, and the testing of principled by the new data they predict. 
JVe show that science, like mathematics, is cumulative. JBasic concepts never die, but 
are exhibited differently at different evolutionary stages.. A two thousand year old 
concept that was once accepted as a law — an axiom for a deductive development, — is 
, . now a speciaFrestricted consequence of today's set of axioms. ThQ concept lives not . 
only "within the new laws and as an intuitive guide for their exploration, but it still has 
a ^ife of its oWn within an appropriately restricted 4pn^ain — a domain* whose boundary 
isfetoday determined analytically by the new Taws, instead of empirically. . 
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We do not attempt to define "light," anymore than we \^buld ^tempt tp'deftne "L:" 
L is for light, L is for life. L is for laugh, L is part of the ABC of communicating by 
symbols that represents perhaps man's most profound abstraction. Light is part of , 
the- ABC of physic s» 

We define neither light, nor various subjective attributes:^ T5tarf4 .with "let 
there be light," an'tl introduce various matheipatical C9nstructs that delineate certain 
measurable properties we associate with light. /We are aware that you are reading 
with the aid of light, and that although this page is at present practically fully illumi- 
nated, you can change th6 situation by closing the book,^ You have seen landscapes by 
sunlight with no trace of Btructure to t8e flood of illuminatignr but you have a\so seen 
dancing spots of light traced lly shafts of sunligh^ op a tree-shaded path, and beams of 
light entering darkenec} rooms through narrow cracks^of slightly open ^ors; or you 
n^ay have become aware of straight Ime^characteristics revealed- in floods of sunlight 
-by the shadows that they cast. You have handled light sources such a& electric lights, 
flashlights, candles, and have seen the^stars as distant sources of light. You have 
seen "the image of your face in a mirror, and the fractured appearance of things (fin- 
gers; silverware, etc) 'partially submerged in water. 

The most primitive constructs for sufth situations are geometrical, and they were 
introduced by the geometer Euclid (about 300 years before the current era) , Euclid 
represented light as something "propagating" (traveling)\long "rays" (straight lines) 
and as being "reflec^" (thrown back in a special way) when it encountered a mirror, 
A geometry of rays and ttd idea that light travels at different s^^eeds^n different trans- 
parent materials serve also to accouyit for the "refraction" (breaking or bending) of a 
ray passing from, say, air to wafer, ' ' . 

Why introduce the idea of traveling? The candle that is consumed as it gives forth 
light, and the monthly electric bill make clear that som§,thing is being used up to pro- 
vide the light. We transform energy from some other form to the form associated^with 
li^t, and the rays we donsider are guides for the flow of ^ergy, 

Jn the! next two Sections, we use the calculjus only to consider the geometry of 
' ' * , I ' ' i 

rays — mostly straight rays^ but also some curved oAes; mostly familiar effects, but 

not necessarily familiar interpretations: we introduce a signed ray (a "shadow form- 
ing ray") to account for shadpws, and some of the rays mdy split into many ('^diffraction") 
to describe some aspects of what occurs when light strikes a sharp ed^e. Sections 2 
and 3 ^pal with geometry, so that we need not mention enetr^flow until we associate a 
magnitude with a ray (as we do in Section 4), Howeyer, for t)ie sake pf the phy3ical 
content, it should -^ways be kept in mind that these rays, in some sense, are the direc- 
tions for ei^rgy flow. In what sense is energy flow associated with light? \Ve do not 
answer this, but merely provide an analogy that is appropriate for most aspects of 
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visible phenomena, " k> J>* 



\ 

There are two familiar ways in which energy travels: if you and a friend are in a 
' swimming pool an^ you splash him, you shower him with water drops, eacfh drop carry- 
ing some amount of energy from you to him, alternatively,^ if you plunge your hand into 
the water, much of the energy of the effort travel^ to him via a wave on the water's 
^^rface. As another pair of examples, you can attract his attention with Energy in a 
packaged form by throwing a ball at hiija, or you can reaph him with energy traveling 
as a sound wave by shouting at him. ( 

Simple observations on propagation, reflection, refraction, and scattering of beams 
of li^t can be interpreted.either as energy traveling as a stream of particles or as a 
coUimated wave; the rays (Sections 2,3,4,7) are either the trajectories for the parti- 
cles, or the normals to^the wave fronts (Section. 5). For the more complex phenomena 
that we may observe readily, only the wave interpretation is adequate: these involve 
/'interference" of light beams (Sections 6,8,9). Two beams of light may intersect and 
produce additive effects iri the region of intersection but then emerge from this region 
unaltered in form: streams of waves show this characteristic, but streams of particles 
do not in general (i.e., some particles end up going in other dk^ctions than the original 
ones)^ For stil^more complex phenomena, some aspects are described more simply 
by a wave analog and others by a particle analog: Which is light? Wa,ve or particle? 
Neither. Light is light. However, since we are more familiar with particles and waves 
(practically the A and the B of the AB^C of physics) we may exploit these better known 
entities in learning and teaching about light. 
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2. Geometrical Optics. ' 
*2.1 Euclid^ s Principles. ^ ^ 

Early observations of light sources (sun. star, lamps) and of the reflections of 
such sources and.objects in smooth surfaces (water, polished metal) suggested that ' 
many phenomena could be described in terms of two "laws" of nature. We call these 
Euclid s principles of propagation [El] and of reflection (E2]: 

[EIJ: light travels along straight lines (rays) ; . • 

IE2J:- when a ray is incident on a smooth plane surface, the incident 

ray. the reflected ray. and the normal to the^urface alltie in * 
the same plane, and the two rays make equal angles on the • 
opposite sides of the normal.""* - *" - ' 

Figure 2-l(a) illustrates [E2J ; it shows the plane of incident, I ' 

Z TT ^^"^ intersection -: 
point (I) on the mirror; the rays are at angles a with the surface normal (N) in" ' 
Figure 2-l(b). we see that the reflected ray (I to P) is the extension of the mirror ' 
ilMSe (S' to I) of the incident ray (S to I). . T ' 




(a) 



(b) 




FIGURE 2-1 



ERIC 



38 



44 • 



If a set of rays diverging from a source S is reflected from a plane mirror as in 
Fi^re>^-2(a), the corresponding set of reflected rays appear to originate from IJie ' 
image source S* as in Figure 2-2(b), T hits as far, as the reflected set of rays ' , 
ffgflected ray system) is concerned, we may replace the mirror in Figure 2-2(b) by * . 
the source and reduce a reflec/ion problem specified by [ E2.] to a propagation 
problem specified by [El]. (This image method was essentially introduced by Heron 
or Hero several hundred years after Euclid.) 

We regard' [El } as defining geometrical propagation in a uniform medium, and. * 
[E2) as defining geometrical reflection from smooth planes. These cover the situa- 
tipns /of Figurje 2-1 ^lid 2-2 as well as more complicated reflection problenTsI 




(a) 



J 




FIGURE 2-2 



A set of parallel rays incident on a planar reflector as in Figure 2-3(a) is re- 
fleeted as a set of parallel rays." If we regard this reflector as consisting of two hinged 
planes, and swing one'^way from the source as in Figure 2-3 (b), then the reflected rays 
are said to diverge ; if instead, we swing the plane towards the source, then the. re- 
fleeted rays converge as in Figufe 2-3(c)» In Figure 2-3(c), the Reflected rays inter- 
sect, v^hile in Figure 2-3(b), their extensions "behind" the mirrors intersect; the first 
(c) is called a real intersection , and the second (b) is called a virtual inters ectibn . In 

either case, the reflected rays appear to originate at the intersection, 
»■>■,' ^ • 

If we have a.spt of rays incident on a complex reflector consisting of many planar 

portions, then we can determine the reflected set by applying [E2], Equivaletitly, once 

we have determined the intersections (real or virtual) of the set of reflected rayS we 

have reduced the problem to dne^ specified by [ El ], It is of particular interest to 

determine the iptersecfions of rays reflected from curved mirrors . But before we 

consider curv.ed mirrors, we introduce a more general law of nature than fEl], 




• . (a) . W ' X (c) 

^ _ , FIGURE 2-3 

2,2 Hero's Principle, • ^ 

Euclid's principles describe the essentials of many directly observable ph^omena 
(ai^d also suggest applications not found freely in nature). However, the principles are 
very restrictive, and their description is very wo^dy. The restriction to plane miri:ors 
, * is removed and the description is cpnvpressed by the jnore general principle of Hero: 

^ [H]: a ray follows the shortest path between points/ ^ 

Before applying [H] to more general situation^ than covered by [E], we use the 
calc/iilus to show that [E] follows from [H]. Since, by definition, a straight line is the 
shortest path beUveen points, we see th*a.t'[H] cov;ers [El] directly. Similarly in con- 
sidering [E2] we need not discuss wriggly pathe. We consider Figure 2-4(a), and, seek \ 
the shortest path between S and P via a point J on the surface of the mirror, froxd' * ^ 
the start wKtake J in the plane-throUgh S and P that is perpendicular to the mirror: 
any displaceraent of J perpendicular to this plane will clearly lengthen the component 
paths Li and L2 . We introduce the lengths and angles of Figure 2-4(b) and seek the 
smallest value of L = Li + L2 (a§ required by [H ] ) and show tha^t this coVresponds to 
y -*a (as required l3y { E2 ] ). ^ ' / / » ' 

In order for » ; 

' < ' * , 

(1) L = Li + L2. = s/h\ + x2'+ Vh^o + (d - x)2 

to be a mjniraum for fixed hi and h^ , %JS^d J. on the surface of the mirror, we require' 

' ^"'^^ ;* dL ' 

(2) e-. ^ = 0 . 

^ * *' ' ' \ 
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' (a) •''(b) 

FIGURE 2-4 



Thus 



- , dL ^ X _ (d-x) ^ JL . l=x = sin'« sin 7 ='0." 

and consequently sin 7 = sin a as in [E2]. 

Equation (2) by itself states only that L is an extremum, or equivalently that the 
ray path'is stationary for first order variations of L. I^^^ver, aince 

d^L _ cos^ a . cos^y ^ n \ 
• dx2 "TT L2 ' 

we see-'thaf L is in fact a minimum. (Of course Hero did not use the calculus: he used 
,the image principle and geometry to show that y = a corfesponds to th$ shortest path.) 

' It is clear from the above that applying (H ] to reflection from a point on a curved , 
surfape is equivalent to using [^2 ] for reflection frjDm the tangent pl^e at the point, 
(and practical applications priojr and subsequent to J H]. have Been ba^ed on (E2] pljuus 
the "tangent plane approximatipi^^"). Figure 2-5(a) shows reflection of a Yay from a ^ 
point oh the concave side of a reflector, and/Pigure 2-5(b) shows tlie' corresponding 
rpfl^Qtion feoiirthe convex side; note the relations of the directions^. Similarly for a 
set of rays Incident on a curved mirrq^, we can construct the reflejjtedl^^t (or equiva- 
lently, tbeii: intersections) by geometry-, note that/both situations in Figures 2-§(c) and 
(d) give rise to the same intersection point ^i^eal for c and virtual for d)^ 

There ^re reflection situations that are not covered by [ H ] but are covered by 
f E21 and the tangent plane aiiproxiniation. For example, if we consider a source .j(S) ^ 
^(jjjgan the.-circimf^^ circle and a diametrically opposed observatioii point (P), 



(a) 






CO 



FIQURE 2-5 ' 



then the geonfetrically reflected ray from S to P via a point I on the circumiterence 
as in Figui^e 2-6(a), is the longest of all such paths (as is clear geometrically). For - 
an arbitral^ point on the circumference, we may write SI = SP cos ^ and PI ^ SP sin0, 
so that L = SP {cob 0 sin 0), and (2) in the form = SP (-sin 0 + cos ^) = 0 gives 
^ "^'J ^ - Oi Si-s i^the figure; but for this case 



d^L 
d02 



-SP (COS0 + sin0) = -2SP cos t*< 0. 

4 * • 



so'that'^L is a maximufii. There are also situations of 'int^^^est covered byn(2) for 
which the second derivatives vanish. For all sugh cases, independently of the second 
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FIGURE 2-6 (a) 
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derivative, the first derivative of the^path length is zero, Le\, the ray path is stationai^ 
for first order variations. To cover all such situations, we replace [HJ by the more 
general principle 

[HM: a ray foUov^s a stationary patfej^ * * • 

We may distinguish two clashes of curved reflectors and illustrate the essentials 
for the case of a concave' cylindrical mirror. The simpler class corresponds to a 
"small aperture" mirror as in'i'igui;e 2-6(b); for this case the semi-aperture of the 




'FIGURE 2 -6(b) ^ 



nurlx)r AB/2 is very small compared tpi:|s radius of curvature a , If I is the center 
point on the luirrbr, and S and P He m the mirror's normal at I , then from*the law 
of reflection, it follows that to a first approximation the rays from S reflected at all 
points of the mirror go through the point P such that 

sr PI sc ' ' 

In particular if SI ^ then the situation cprrespondS'to an iiicident set of parall 
rays ^s in Figure 2-6(c),,and we obtain simply 

■' ' . ' Vl ^ } . 
J' ' > » • ' ^ ^ 

<H , ' a " * 

'Thus for the small-aperture mirro;x all reflected rays intersect at ^ (the focus ) 





\ 

HGURE 2-6(c) 



^he more general problem of ^reflection from a miri:or with arbitrary sized aper- 
*ture is illustrated in Figure 2-6(d) for a set of rays from a source on the axis of a 




.•FIGURE 2-6td) ^ , ^ - ' 

semicircular mirror, Figiire 2-6(b) shows only to the situation in the, vicinity of the 

* axis. If we rotat^e these figiures around their symmetry axis, the situations^correspond 
to reflection from portions of spherical mirrors. If the distance of the source from 
the reflector in Figure 2-6(d) becomes infinite (parallel set of rays incident), then the 
envelope of the reflected rays (the locus of intersections of neighboring rays) is an * 
epicycloid" (to be discussed subsequently); the cusp of this curre*4^ at I • 

^ The envelope of the rays is called a.caustic ; we may de^al with a causti& surface . 

a caustic line; or a caustic point : the last is also called a focus. For specific sets of 
^ rays incident on specific curved reflectors we could determine the caustics geometri- 

cally; however a geometrical procedure is usually too tedious. Instead, we apply the 
^ calculus to a quite general situation, and determine the caustic that specifies tl*' re~^ ^ 

fleeted field for a parallel set of rays incident on a cylindrical mitror* We restrict . 

attention to the plane perpendicular to the cylinder's generator^^^^that|he problem is 
\ essentially two dimensions^, and derive the correspoftding^in^austic. Since the re- \ 

fleeted rays are t^gent to the caustic, once wCknow the caustic we specify the^. * A 

* refl€?cted field by means of [El], ■ o h ^ ' t'o 

^ ■ ' ; " ^. '-" ■'. ' - -^^^ 

2,3 '^Caustics, • / ' 



We consider a set of rays parallel to the x-axis incident on a fwo-dimeftsional 



I 



mirror. For each incident ray, we could deteimiine tfie corresponding reflected ray 
; geometrically. Instead,^ as the initial step for a subsequent |dQtvelopment, we consider 
an pljservation point P(:^,y) on the same 3ide of the reflected as ^e source, apd apply 
[H] to relate^e incident ray that strikes-the mirror at I(|, 7?) to the reflected ray # ^ 
through P i see Figure 2^-7, ^We could specify the point I ancl the incident ray that 
strikes it in term's of the parameter of arc length along the curve; however, it is 
mozTe; convenient to use t^e. angle that the incident ray make^itlvthe normal at I as 




FIGURE 2-7 



the par^eter, the an^e a such that tan a is tlve'slope of, I^. The length of the inci- _ 
dent ray measured itom the yr3pds is | ; the lengSi of the reflected ray from I to P is ^ 
R = y/ - x)2 + (tj - y)^ , and its inclinatioAyto the k-axis is '®. 

Tffe total length from the y-axis via I to.^ equals ' ^ 



(3) * * L = I R = I + y/ii-xf + (Tj -y)2 . 

Differentiating, we have^ 

(4) 



^where the prune indicates differentiation with respect to q . Using.[Hl, essentially as _ 
for (2), we equate U to zero to obtain ' ' ^ j 



(5) 



sin €> 



Iff l +^cos @ 



2 sm -g" cos -g- 
2 cos2 ^ 



tan. J 



Thus from the chain rule we have 



tan 



dTj 



(6)'^ .■. . 

^ Since ^ is the slo|)e of the tangfent of the reflector at the point of incidence, ^ 
. is the slope of the normal , anol&quals tan a . Consequently / 



from which- 
.(8) 



— ^ = tan 7 = tan Q! , 



.4- 



J: ' as could^ave been obtaine^irectly from ] E 1 on^ inspection of Figure *2^7. 
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'The equafton of the re£)ected ray arising from the ra> incident at an angle a with 
^ is ' ' / 



(^) * . ' ri y ^ ii' x)nan 0 = (| - 5^'tan 2ot , ^ / ' 

^which we rewrite as r ' 

'M^^) . * g(0L) = (I - X) tan 2a - (??' - y) = 0 . 

This specifies the set of reflected rays ^Ci3^esp1)nding to a set of incident paraUel 
rays. The parameter a describes not only the cur^of tlie reflector '[^(ol), t?(q:)J^ 
it al^o picks out the ray incident at |, n and the corresponding reflected ray (10) that 
reaches x,y. The point of intersection of two neighboring^rays g(a) = 0 and 
g(a + 4a)3f .0, corresponding to a and a + Aa , is determined by g(a) "= 0 and 

;~Aa ^ ~ 0 r In the limit Aa — 0 , the point of intersection of the,rays 

falls on the envelope and is specified by the simultaneous equations 

(^1). -^m^ or g» (a) = 0 . I 

Differentiating (10), we obtain 

(12) ^ g' = i^iBXiZa 4- - - , 

cos^ 2 a t>f * - ' 

* ' * 

and with (Tfwe 'eliminate n ' =, • 

tan- a ■ ' ^ 

(13) ' ■ g'.= ^' ('tan2a + T-i-) . 

^ ■ \ tana/ co.s22a . _ 

* ' 1 * ' 

Sinpb.tan"2Q; tan a + 1 = * , we reduce (13) to 

^Z'" ^ -^w COS" ctOt 

>^(14) . V g' = ^ /^'cos 2g~ _^ ' 

V , >rf<.-. /V_, cos2 2a \ tana ''t^ ''V ^ .-..I 

Thus, Since (ll)_requare-s g' (af'= 0 , the X-cp&^nate of the point on the envelope is 

, : * . '2 tan a ' 

* which we may rewrit^e in various equivalent forms, e.g., , " ) 

. (16) < ^ X = I' - -iL- cos'2a V • ^ \ 

^ We obtain the^correspondih^y-coordiijatejiy using (16) to pliminate 4 - x from (10): 
thus ,g = ' "I*; cos 2q^ tan ?a - (7? - y) == 0 , and consequently) * , ' 

(17) ^ i y = 7?: - sin 2a / i * ^ ' 

- ^ ■ ' , / ; 

Equations (15) and (17) .(from whicl^ we^cotdd eliminate d) specify the caustic > 
curve, the envelope oi the reflected rays (the locus of ^he intersections df neighboring \ 



rays). Given a specific reflector we can use its parametric representation to elimi- • 
nate i and r? , and thereby determine the caustic%xplicitly. We illustrate"this for two 
simple reflectors, the parabola and semicircle. 

Parabola . We consider a set of rays incident on the parabola 



(18) 



.2 = 



-4p| 



as in Figure 2-8. The parametric equations of the parabola in terms of a , where 




v^FIGURE 2-8 



tan a is the slope of the normal, are 



(19) 

arid consequently 
'720) . 



rj = 2p tan a , i -p tan^ ct , 



cos^ a 



, Using these expressions for tj and 17' in (17), we have 
■ (21) 



2'tan-a ' 
cos2 a 



y = 2ptana- — sin 26? = 2pYtan a - ^i^) = 0*. , 
^ > cos2 a > J ^ . . cos-^y > 



Similarly, we use^the corresponding^expressions foi; | an<J ,4' in (15) to obtain 

t|^^refl< 



^ ' ^ Thus the envelope of the reflected rays* is ^ » ' \ 



(23) 



X =^^P , ^ ,y 0 7' 
1 V 



i.e., the focus of the parabola as in I^igure 2-9.*The focusing.^jpjJerty of the 



\ 



FIGURE 2-9 

parabola accounts for its many applications (as telescope mirror^, microwave and 
sonic "dishes,''' etc.) for collecting practically parallel radiation (th*e rays from very 
^ distant sources )^by reflecting the incident rays to a snjall appropriate detector placed 
at its focus, Siriiilarly, parabolic reflectors are us^ed for the inverse problem of con- 
verting the Radiation from a source at the focus into a parallel beam of rays. ' 

The above exanipl^is practically trivial in that (23) cjould have been obtained by. 
much simpler procedures than the one we followed. In the next example we follow Mk 
essentially the same procedtoe to obtain a far less obvious resolt. 

Semicircle . The param^ic equations of a circle of radius a* for the problem of 
Figure 2-10 are , 

(24) r/ = a sin a , | = a cos o: , ^ 
and consequently . > ^ . * 

(25) ^ 17* = a cos a t ^* = -a sin a , 
From (17) we* then obtain for the caustic ' ^ 

(26) y , y = a sin q: - a cod a ^^^^^ = a sin^ a , 



r ' 
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and from (15), . . ' e> - 7 

a si n cos 2a _ ' a cos a A;? ^,^o 

(27) x = acosp^- . 2tana ~2 <^ +'2sin2o^). 

SqUaririg^and a&ding (26) £11^27), we obtain . ^ 

(28) r: (X? + y2) = 1 + 3 sin^ a •= l + 3(^)s , 



a 

the equation of the epicycloid traced by a point ^en a circle of radius j rolling on the , 
V 'a • * / 

outside of a fixed circle of radius ^. ' , ♦ 

• . 4 -. ' 

The cusp or focal point of the caustic is. at x ^ ^ , y = 0 ; this corresponds to 
g" - 0 , and occurs at a = 0 . The r^s incident near the center of the mirro^ (q ^ 0) 
aije,known as paraxial rays of '*smalV3!perture** mirror theo'^; only these give rise to 
reflected ra3^ thaf appear to originate at the cusp ^ . [ Fpr parallel rays incident on 
the Fjarabola, the entire caustic consists of the point focus; similarly for a source at 
one focus of aniellipse, all reflected rays go through the other fociis (hence the l5|)el).] 

Virtual Caustics . In the above we considered reflection from concave mirrors; for 
such cases the reflected rays intersect and the caustics are real in the sense de- 
fined in Section 2.2, Similarly for incidence on a convex reflector'the extensions of 
the reflected rays behind the reflector intersect on a virtual caustic. The identical 
caustic'curve specifies reflection from either side of thfe mitron. Figure 2=11 shows 




FIGURE 2-^11 



the i^ituation for incidence on a convesC parabolic reflector, and/ Figure 2-12 shows 
the an^ogous situation for a semicircle. Figure 2-13 sho^yS the geometrical method 
of constructing the *epicycloidal caustic of the semicircle. . > ' 

Sinqe tho caustic of Figure 2-13 is the envelope of tjie set ol exteiided reflected 
rays, it is tangent to all members of the family. From the figure w^ see that thesis- 
tance from the mirror along the ray extension to its point, of tangency with the catistic 
equals § cos a, (Without the geometrical construction, the result^foUows on sub.-r 
tractiiig from ^ ^ (the totaj ray extension from cylinder surface tq x-axis) the , 

' \ 
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siif2a ~ 2 cos a sin d length of extension between caustic and x^^s.] 
, \ THds neighboring reflected-rays of real length R (where > is th^istance frcim the 

,\\' inirro^) appear to diverge from a source (ifeir point of intersflption) at a distance 
• \ ' \V.'^ ■*■ f ( along thfeir extension." > ^ ^ 

^■V4' - ' ■ • . ' ' ^ 0- ^ ■ 



Since the reflected rays arle tangent to th^ caustic, we maj^^reat the caustic as the 
evolute (the l|)cus of the centers of curvature) of a s^^"^ curves which are orthog- 
cmal to the ijays. These curves, the involutes of the "caustic, are called the eikonals 
or eikonal curves in ray theory; the radius of dlirvature at ,a, point P on such a curve 
equals^R + cos a . The rays (the orthogonal traiecta^s of thfe^konal curves) ase 
tangent to the caustic and normal to the eikonals, and this provides a geometilical con- 
struction for the' eikonals: they are trace^ by the points of a taut string as itf unwinds 
from the caustic. ^ 

2.4 Shadows . . ' 

^ ^ in the discussion of (3)ff we restricted consideratijon to an observation point P , 
lyixig on the same side of the reflector as the source (the "ht side'' of the refl.ector). 
If we drop this restriction, we obtain <an additional s<$|||ion of = 0 with L' a^ gtven 
ill (4),i.e., % , . / % * . . 

• • si'**; 

^ ^29) i\ ^ . ' . . L» =^ 0 if © = TT , ' ^ • 

wherfe the^geometry is shown in !Figure 2-14. Thus in addition 5?he geometrically ' , 

■ •/ >* ' . " ^ ' " 

0 = TT . 
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feRected ray shown in Pig|Jire 2-7, we»^ee from L' .= 0 (i^., from [H']) thai the inci- 
dent ray also gives rise to another ray -ona traveling along the 'original direction of 
incidence. Were the reflector "abseht, we would iftterprgt this* ray as the incident ray 
itself (i.e. the situation of [El]). However,^we insist on then)resence*.of the reflector 
and seek a physically significant interpretation of ^th^ rays cdrrespondinglo (29)^. 

Whence interrupt a broa^d beam of light by a*mixror, we notice essentially tivo 
effects: because of the mirror, there is^not bnly sdn'e light-observKl in a region qf 
s^ac^ outsjge of the original beam, ^ there is a'fed some light^missing f^m a itegion* 
of space, originally filled by the beam before we ifi^erted tjie obstacle. Were \ye inter- 
ested solely in the original team,* then we might simply say that some dtihe light has 
beSn ^Tjent" from its original direction (reflected) and let it go at that. -Howeye^ in> . « 
^^order^to ultimately specify the full effect of .the obstacle analytically, we assign, it a 
more positive roTe^ We say that the incident rays'"excite** the. obstacle to prodiice not 
only the set of 'reflected rays but also ar^^ of jshadow forming rays parallel to the ' 
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"missing" incident rays in the dark region of space. It is theselk^dow forming rays 
that we read into ,(29).; these must cancel the incident ra_ys on the "dark side" of the^ 
mirror to "create" the geometrical shadow. (This idea of shadQw forming rays may 
be hard|to reconcile with mental images of the reflection of rays based on a ball ^ 
bovincing off a wall. However, were we interested in specifying the total eJEfect of the 
wall in the ball-^wall problem^ \y;e CQuld also do so in^terms of reflected balls and 
shadow forming balls,)' . * ' . . , 

To make the role of the obstruction more explicit (and to set the stage for ouj 
subsequent' discussion of scattering), we introduce a symbolic representation fQr the 
rays. We represent the effect of an incident ray by Ej , of the geometrically reflected 
ray by Eg , and of the shadow forming ray by E^* Wfe represent the total effect E^ 
corresponding to a ray E* incident on a i/eflector as 

( Eg • in lit space 

(30) Et = Ei + E: . E = < ^ ■ ^ 

( Eg in dark space, 

Thus in the lit space the total effe.ct is' E^ = E^ + Eg as shown by the t\vo rays onthe 
left hand side of Figure 2-14, On the other hand in the. dark space we have Et = 
Ej Es corresponding to the dashed ray on th^ right hand sid^^ Figure 2-14; in ^. 
order^hat Et represent the physical situation of the geofnetrical shadow, i.e.,^ , 
order tha^ E^ vanish, we^requ ire ^ , ? ' ^ 

(31) • , ^ JE, =-.E, , 

We take'(^l) a^O^s^pplementary assumption to [H'J : the first sblutioji (0 = 2a) of 
U = 0 accounts for geometrical reflection (and we subsequently determine a ma^i- 
tude to be assigned to such tays); th^ second solution (© = tt) plus (31)*' accounts for 
shadow formation. - ^ / 

%* . « 

The symbol E in (30) represents the scattered>'part of the total effect E^ = - 
Ei + This is the part of 'Et that we may regard as originating at the obstacle to 
Ei , or as outgoing from the obstacle^ * . ' < 

If we consider a system of parallel rays fricident on a convex semicircdlar. cy- ^ 
Under (or equivalently on a full circular cylinder), then thecorresponding scatt^refl ray 
system (reflected pjus fhadow forming rays) ii as pketchednn Figure 2-15, 

• .The family of curves perpendicular to these rays is th^ corresponding infinite 
set of eikonals. Figure 2-;L^ plus its reflection in the x-axis shows several of th^se ^ 
ctirves. These curves may be obtained geometricatly from the caust^ca (fhe caustic' 
for the shadow forming rays is the point at x = -«), or by constructing the normals 
of Figure 2-15 geometrically, or analytically, W^ consider^an analytical derivation 
in a followi^ section. At larger and larger distances from the scattereif the eikcwials- 
of Figure 2-16 become more and more circular, \ v ' 
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2.5 Edge DiffractedRays. 

There are adStioiikl sets of rays implicit in Hero's principle, an^^elr utility^ 
has been shown by the reclsnl investigations of J. B. Keller. In particular we consider 
^dge diiSfracted rays arising when, a ray is incident on a sharp edge (which '^breaks up*f 
or "diffracts" an incident rajj). ^ In order (o motivate introducing such rays,, let \is 
review the precedping materiaL ^ . 



We have discussed reflected rays and shadow forming ravs . and we saw in cpn- . 
nection with the semicircular cylinder that both kinds of rays were required to obtain , 
a complete coverage of space by scattered rays (or equi;;^alently to obtain closed scat- 
tered eikonals). However if the scattereras a strip as in Figure 2-17(a)>, such rays 
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FIGURE 2-17(a) : 



alone do not cover space, which implies that the scatterer^s influence is restricted to 
the two directions shown in the figure. To construct a scattered Vay system that covers 
all space, we Introduce the edge diffracted rays of Figure 2-17 (b); these rays are in- 
cluded in [H], i.e., an incident ray striking the edge is diffracte<|'to P via the shortest 
path. ' " ' 




FIGURE, 2-17 (b) 



From Figure 2-17(a) and 2H7{b), we see th'at there are- es'sentially three different 
cases that arise tor a fully illuminated strip; these correspond to the three different 
Qbservation points of Figure 2t17(c), An observation point at' Pj receives two 



• P2 



FIGURE 2-17(c) 



diff?kcted rays; P2 receives one reflected ray and two diffracted rays; P3 receives 
one shadow-forming ray and two diffracted rays. In a subsequent section we show that 
the magnitude (of energy flow) associated with a diffracted ray is in general much 
smaller tfian the magnitude of the other rays in Figure 0- 17(c). If we assume this 
result for present purposes, we neglect the diffracted rays in the regions correspond- 
ing to P2 and P3 and obtain the scattered ray system of Figure 2-1,7 (d); this figuj-e 



FIGURE 2-17(d) ' 



shows only the "strongest" scattered ray at each observation point. A cpr,r(Bsp9nding 
'eikonal curve nonrial to the rays of Figure 2-17 (d), is shpwn in Figure 2^I7(e), and it 
is clear that such surfaces becomejtnore circular witl\ iftcreasing distancelrom^the 
scatterer. * , - - ' 
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. FIGURE *2-I71(^);,^ ] 



The variQus rays of Ftgure 2-17 correspond only to the scattereci ray, system, 
i.e., to the effects in sgace arising" from something tjiat obstructs' the incident rays; 
this figure does not take into account that the observation point is also reached by an 
incident ray. In particular, as discussed for equations (30) and (31), the incident.rjays 
and shadow-iorming rays cSncel in the shadow region corresponding to P3 . Thus the 
net effect in the shadow 'region must, arise from the edge diffracted rays as in Fig- 
ure 2-18; such' effects have been discussed in detail by J.B. Keller. (Bright areas in* 
the. shadow region of obstacles with very regular edges were first commented on by 
Grimaldi, 1613-1663,) ' , . (' 
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FIGURE 2-18 



For present purposes," we consider only the caustic of the edge rays for the anal^ 
ogous problem of a disc. Thus if a parallel set of rays is normally incident on a 
circular disc as in Figure 2rl9, each point of the edge gives rise to a "full fan of rays 




FIGURE 2-19 
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nprmal to the e<^ge ilfhat point.^An off-axis^bservation point receives edge rays, 
only from two points of the citrcumfer^nce on the disc, i.e., from the ^ametric ally' 
opposite points cut by the plajie containing the observation point and. the discos axis. 
However, a point on the axis of the disc receives*edge rays from th^ entire circiun- 
ference: the axis.is.a caustic ^olthe .edge , rays «.Thus.the center of the shadow-bf^a - 
nopnally illuminated circular disc should show aJ)right spot, the Arago bright spot , 
or 'Poisson bright spot (asvpi^dicted originally about 1800 via a wave argument — a 



1)lg ktgument) 



For the circular c^c', the line caustic of the edge rays is the envelope of the 

planes normal *to the, edge of the disc. For a disc of general shape* (an arbitrary 

planar scatterer) normal to the parallel incident rays, the corresponding caustic of 

the edge rays is a cylindrical surface, .the envelope of the planes normal to the edge, 

) Sinde two such planes intersect^in.a line normal to the disp, the caustic cylhidrical 

^ • ' . r i- 

surface generated by the lines of intei^section is al^o normal to the plane of the disc, 

" The cross section of the caustic cylinder cut by the plane of the disc (or as viewed on 

a screen in the discos shadow), is the line envelope of normals to the edge in the plane 

of the 4is.c;^ it is the evolute of. the edge. . , . . 



; In particula* for an elliptic edge 
(32) 



a2 : 



the equation of the evolute, the four-cusped curve sketched in Figure 2-20, is 
(33) . ■/ 



1 • 1 • 1 
(ax)3 + (by)3 (a2 - b^)^ 



FIGURE 2-20 



To derive (33) we write a point on the ellipse parametrically as 



0 



(34) 



^ = a cos (p t ^T? = b sir^(p < 



The corresponding normal through x, y is specified by 

ax by 



(35) 



cos(p .sin<p 
and the derivative with resi>ect to (p gives 



3? + b2/ = 0 , 



(36) V 



ax 



n3, 



0 



COS^ ^ %)X^ ^ 

Substituting (36) in" (35) to eliminate either x or y, we see*that 



(37) 



ax 



T - 



Consequently the locus pf the normals is ' 

a2 - b2 o 



= -by t a2 



b2 



(38) 



X = 



a 



cos^^p , 



b2 . q 



Eliminating (p from (38), we obtain the required result (33). [ Note that for the geo- 
metrically reflected rays we started with a set of lines ^(the rays), determined their 
caustic, and then identified the caustic as the evolute for a set of involutes (the eikonal 
curves). For the present case, however, we started with an involute, (the edge of the 
. disc) and determined the corresponding evolute, the envelope of its normals (the loclis 
of the centers of the circles tangent to the involute, the^lociis o^the centers of * 
curvature).)/ ^ 
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If we visualize an experiment in which we start with a circula^r disc and gradually 
convert it to one of elliptical cro^s section (or equivalently if we rotate the circular 
disc so tjiat it is no long.er perpendicular to the incident rays), then on a screen nor- 
mal to the direction of^ incidence .the bright spot changes to the foUr cusp'ed evolute of. 
the ellipse.^ Such caustic sections were photographed by Coulson and BeckneU in 1922, 

The associated scale factor : In the above we discussed rays reflected from surfaces 
and rays diffracted by edges. The edge ray? as in Figure 2-17(b) are drawn radially 
outward from a point on the line representing the edge, bufthe scattered rays of Fig- 
ure 2-15 for th^ cylinder are not radial. If we visualize the cylinder becoming nar- 
rower and narrower we might expect on the basis of our remarks for^dge rays that 
in the limit the ray system of Figure 2-15 could be represented as a set of radial 
lines as in Figure 2-21, . \ 




' , ' ' FIGURE 2-21 

The situations pf .both Figures 2-15 and 2-21 are covered by [H^ » and both cor- 
respond to scatteringMby a Circular cylinder,' In order to distinguish them we must 
associate a scale factor for length with a ray. To do so, we could'assume that in addi- 
tion to the g^metrical property assigned to a ray by [HM i a ray of light (of a single 
color) has an^^sociated length X that is independent of the length of the ray 
path, *Wa^couia then distinguish the two different scattering situations for the cylin- 
drical obstacle of Figures 2-15 and 2-21 as follows: the ray system, of^^igure 2-15 
corresponds to a very large cylinder a » X , and thd ray system of Figure 2-21 
(Jdrresp^nds^to a Wry small cylinder V^'< ^ ^ . , ' | . ' 

The existence of a&i associated length might ha^^been guessed (from Grimaldi's 
experiments on light diffracted into shadow regions) but was not. We show subse- 
quently tiiat the ,required scale factor emerges naturally as part of a more general 
model for such phenomena. We mei^tion the matter now partly in anticipation, but 
primarily to stress that the present model is incomplete. 



a. Sneirs Law, Fermatjs Principle. ^ 

In the preceding sectibna-we-considei^pd a set of rays incident on^ reflecting sur- 
faces, and used [HJ to det^miine the reflected set of rays. We now extend tl?eMevel- 
opmenttb partially transparent surfaces and consider in addition a set of transmitted 
rays, A transmitted ray does mjt lie in general aloi^ the extension of the correspond- 
ing incident ray, but makes ica appdJpriate angle with the ray extension (e,g,, as in 
Figure 3-1); this kind of 'J^reak" in the ray path is called refraction . 




J 



FIGURE 3-1 



Observ^ions and studies of the broken' appearance of a rod partially immersed 
in water, and of a beam of light traveling partly in air and partly in water, go back to 
Euclid and Ptdlemy (second century of this era), but the complete description of such 
effects was first given by Snell (1591-1626), As the appropriate analog of [EJ for 
reflection, we have Snell 's law of refraction : \. /^.^i— 

]: A ray incident on the smooth plane interface between two trans- 
parent media gives rise (in addition to the reflected ray) to a ^ 
yefracted ray on the othel: side of the interface. The incident ray, 
the refracted ray, and the surface normal lie in the same plane, 

' *"-^nd the two r'^ys are on opposite sides of the normalT/^Tfi? slne^ 

* • * 
of the^kngle (/3) that the refracted ray makes with the normal i^ 

proportional to the sine of the angle (a) of inciden^. 
From (S], we specify the direction of the refracted ray by 



(1) 

♦ 

or equivalently by 
(2) • ^\ 



/i2Sin)S = /ij sin Of', 



\i sin^ = sin Of 



The constants /ii , and are called the indices of refraction , and = ^^r- is called 



Ml 



the relative index of refraction ,* The situation is shown In Figure ,3-1 for mi < (^^ 
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assifmed in all that follows): the ray travels from S to P via a point I on the inter- 
face. The M's ar^ Dhysical constants which specify the essential physical property .of 
the media for the topic at hand, they may be measured experimentally, and we'assume 
they are known. In particular for light (yellow ligjjt) passing from air to water we ' 

We may apply ( S ] to such problems as a point source above or.below an air- 
water interface. In particular the caustic for the system of refracted rays can be 

found by the method of Section 2. If we consider a point source under water (u = -) 

3 • * 5 \ 3/ 

and the rays for which sm /3 < - , we can show that the virtual cajkstic for the rays 

refracted into air is the evolute of an ellipse, and that4:he eikonals are parallels of an 

ellipse. {Axi object under water, viewed along different directions from above, appears ^ 

to lie along the corresponding rays tangent to this virtual caustic.) " 

♦ * 
, ^ Fermat assumed that in a given medium^light travels with a velocity v inversely 
' proportional to the index of refraction (^v - ^ , where c is the velocity of light in 
vacuum^ and rewrote (1) as 

n\ sinoc _ sinfi c c 

lichen derived (3) from the following minimum principle called Fermat' s Principle : ^ ^ 

^ [FJ : A ray'takes the least time to travel between ,two points. 

If the total ray path consists of two stra\eht lines Li and L2 in two media with 
velocities equal to vj aifd Vg respectively, thfen the correspondmg travel times are 
ti = with i = 1,2 ; from [F] we see that must be a^minimum, Fermat's^ 

princ^le [F] not only replaces the clumsy [SJ (the >yay [H] replaced it also 

-incluSbs [H] as^CKe^special caseHvhere v^ th^ taints S and p' are on the ' 

same side of the interface, . V 

We now use [F] and the geometry of Figure 3-2 to derive [S] , essentially* as we 

used [HI to derive The time taken to go the distance Li -from S to I at a 

^ - L Lo 
velocity Vj is tj = and, similarly, t2 = ^ is the travel-time between I and P . 

at velocity V2^ in medium' 2, Thus [F],ieji|ur^that , ^ 
(4) ^ hi ^^hl ^ y^t^ ^ ^ + yh| T'(d 1 

v/ V2 ^ Vi V2 

be a minimum. Differentiating (4) with respect to x ^tnd equating the'^result to , 
zero i.e., . • • . 



^ X ^ ' d - X ^ sing _ sing _ ^ 

7iy/h[T^ V2 yh2 + (d - X)2 • ^1 " ^2 • " ' 

■ • • ' r . ' • 
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* 

we obtain (SJ in the form (1), (2) or (3): ■ . 

i * ■ * 

(5> ' sina - — sin^ = sin*^ = /Isin^ . '^'^ ~ , 

V2 Ml ^ 

If = 1, and S and P are both in medium 1, then (5) reduces to [E] • 

sit is^Clear from our di^ssion of the replacement 'of [H] by [UM in Sectioif 2, 
thatwg should also generalize [F] by replacing leSst time by stationary time, Equiy- 
alently, if we define the optical path length ^to be piL, then aslh^ analog of [H'] jff^ 
•take * ' ' L ' 

. :* a ray follo\Vs the stationary optical path between points. 

Unlike' [S] , we niay u^e [F'] and (5) for refra^lij^n at curved interfaces. Thus 
we could now consider the refraotion analogs of the reflection i^roblems we considered 
^^eviously, For eXaniple, we could determine the caustic for^the two-dime^jsi^nal 
problem of a se|^f parallel rays inci^nt hose-on along the axis of a convex^bi^icir- 
cle of. radius a "capping" , a strii^ (whe^e^the semi^cijcle-^stri^^ ^. 
by M), or the three-dimensional analog of incidence along the axis of a hemispherically 
capped rod. For this case, the caustic of the rays refracted at the semicircular inter- 
face can be obtained by essentially the noethod <$. the previous section; th§ cusp of the 
caustic lies on the axis at a distance Ma/(fx-l) from the interface, so, that e.g,, for ^ I 
/i* = 4/3 the cusp is at 4a. Similarly we can obtain the caustic for the rays that . « 
undergo. t^yo refractions for incidence on a circle (or sphere); for this case the cusp 
is* on the axis at if distance (3M-2)a/2(/i-l) from the first interface, e,g,, for n = 4/3, ^ 
the cusp^is at 3a, ^ (Would there be a shadow ? Have yeu ever'illmninated a cylindrical 
gl^ss of water with a-flashlight?) 



Rainbow caustics, Newt6p (1719) shoWed that white light could be regarded as made- 
up of light of different colors, eaclx specified by a different value of some physical 
parameter (say w), and that in general the relative index .o'f refraction between two ♦ 
media (iepended on color, m =; m(w) . Thus a ray of white light incident at an angle a 
on an interface may be freated as a set of coiiy^identVays M different colors (u;) each 
being Refracted a dlfferent^gle i3(a)).as det^mined by ^ed^^rre^ponding 'index of < 
refracfion ^{u) , Consequently, a teingle ray of incident white li^ht becomes a fan of 
^coloreii rays {the spectrum) on i-efraction, the different colors appearing at angles <^ 
determined by . . . ' 




^ ' FI(?URE 3-3 . * y • 



^ - 9 ^ 

Relation (6) is strikingly e;diibited in thq rainbow formed by^sunlight incident on 
spherical water drops. In the following ^e use thg methods of calculus to determine . 
the' Angles of the primary rainfaovr ana secondarj^ainbow for circular cylinders and" 
spheres, ^ , 

A ray incident on a transparent circle (such as a cylinder of wate^ in air) gives 
vise to,an jnfinite number of rays. Some of «ipse are shown in Figure 3-4; initially 
we consider the ray p. If a system of^rallel* rays is incident on the cylinder, then 
we deal with incident rays making all angles a (fromvO to 90^ with the cylinder's 
normal, aAd to each cprresponds a different p(a) , We want to show that in the vicinity 

62 68 ••• 



^"FIGURE 3-4 



of some particular value of the angle a (say cx^) the rays p(Qfs) iwill have a caustic . 
(willlS^lfocused*'), or equivalehtly that the angle of emergence of p has\ stationary 
value corresponding to a,. > ' v 

' : ' i • ' ' 

" The. primary rainbow corresponds tQ rays that have undergone wvo refractions 

L ^ • ^ y ^.^ * . 

and one internal reflection as shown in Figure 3-5^ We now show that the angle <P 
(the angle between the emergent ray, and the incident' ray) has a st^ti?onary value (ps 
and ekpress in terms of the relative index ^ • ' ^ ' * 




FIGURE 3-5 



Vrom the figure, we Have: 



(7) - - 

Equating ^ to zero \ve obtain 



(8) 



I = 2^ -a! 



2f =vl 
da 



I ■ 



V 



-Iii addition, by differentiating the law ofTfifraction ^tsin/S =.sinQ!, we have 



so that (8) and (9) yield 



McoSiS ^ = COS& , 



(IP) ... ' ' X ^cosp = 2 COS a . 

Thug froiir'(10):and (5),, we obtain 

(11) , ^ 3cos2a3 ^= ^^^2.-1 



' 5 



whiph^etermines the stationary value of themngle^ incidence ,;^and consequently 



.the corresponding values of and (p^. In particular; 



(12) 



Fpr yellow light, ^(w) « | , and consequently (p^ « 42 ^ for the colors red through 
blue> ^e corresponding values of (p^ decrease through 42* . * > 

This result for a cylinder also bolds for a sphere, ajid is therefore' basic to the 
rainbow formed whfen sunligl^^^y||^inates a region of air containing many water drops, 
For one water sphere, if the sun is in back of you and youcan sec the ray through P 
of/igure 3-^5, uig colored rays will be at about 42*^ with respect to the direction of 
^tdence (in the plane of the sun, the drop, and your head), Iflhe sunlight illuminates 
a large. numl:>er of drops over a very large volume of space, tl^n you will see the 
familiar rainboW arc. ' < . ' 

secondary rainbow, corresponding to two internal reflections, we deal^ 
with the geometry of Figure 3 -'S.^ We now have 



(13) 




FIGUBE 3-6 
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Differentiating (13), and using (9). and (5), we obtain 



(14) • - 

and consequently 



8/xl 



(15) sin 2 

For ju = |,weh«e Vs« 51\ More "generally,, Jor n internal reflections, we.have 
(16) 



^« n(n^+2) • 



Using the fact that /x(u)) increases as the colors go from red to blue, one can describ 
the appearance of the primary and secondary arcs in space and account for the differ 
ent orders of the colors in the two cases. 

Stratified Medium : Jf'we apply the law of refraction to a ray traveling throifgh a set 
of parallel slabs as in Figure 3-7, such that e*ach slab has a different index of refrac- 
tion, we obtain • '-^.^^ 
(17). Mo sin ^0 = Misin ^1 = /i2^n 02 = • • • - constant = c ,* ^ 
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FIGURE 3-7 ^ . • * 

Similarly for the limiting' case of a continuum whose index of refraction is solely a 
function of x, we havf ^ ^ . / " 

(18) • M(0) sin 0(0) = M(x) sin0(x) = c , 



U^lng g'-= tane = , ^^^^ , we obtain 
^ y 1 - sin2 0 • 



(19) 



from which 



dx 



c 

Jf_ 



(20). 



dx 



Integrating (20) vb^fl^een 0 and x, we obtain 



y - yo = */ 



^^^^ 



w- 



.J 



Thus in terms of /i(x) we have derived an equation to -specify the set of.r^?s that 
"fstk&i a^O , yQ and arrive at x, y. ' ' ^ 



As an illustration, we assume 



(22) 




1 +bx * 



where b is assigned parameter. To ev^uate the integral (21) in terms of (22), we 
change the variable to (p , such that ^ . 



(23) . 

and rewrite (21) as 
(24) 



c(l + b^) = sin^ , 



sin"' 



bx) 



sm c 



sin^d^, 



Thus on integration, >ve obtain 



(25) 



. • ^ bJ "^u ^b— J : (S^ ' 



.i.e., the equation of a circle *of radius -r- whose center is-locateci'at 

Rotating the coordinate frame of Figure 3-7 (for convenience in the following 
application to rays in the atmosphere), we show ray paths ^n Figure 3-,8 for (25) with " 
b < 0 andj^id^ b > 0 . ' ' - , .-.T ' \ ' . ^' . 
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b>0 



FIGURE 3-8' 



The a^ve results serve to accoimt for mirages. Normally the density*of the 
atmdspiiere decreases gradually with increasing altitude; the index /a , which depends 
prin^arily on the density, also decreases.gradually. However, over a cold extended 
surface the density and ^ may decrease yapidly with height. An object on the surface 
m&y then be. seen at l^irge .distances by nleaiis of down-curving rays as' in Figure 3-9 



/ii \ ' 

T"^ X 




FIGURE 3-9 



(in which the cuiyature is greatly exaggerated). The eye^^ghts along the stifle of the 
ray^s arrival, and one imagines that the ship lies along the linB-<5f-sight. On a much 
larger scale and with normal decreaS'e ot >i "With altitude. Figure 8-9 accounts for 
our seeing the sun by refraction after it has^p^ssed b0low the horizon. 

r . . , \ , , ' ' ^ ' . 

A more common mir age occurs ove r a hot extended surface when ^he density and 
^' first, increase and thendecre^se with increasing height. For .such cases^he eye 
.may see the object by an upcurving ray as well as 'by a straight ray as in Figure 3-10 



4. Kepler-Lambert Principle * 

In Section 2, we assumed Hero*6%principle [H'] that the ray path be stationary, 
and used the calculus to reveal some of the implicit physics. Except for the discus- 
sjlonjof sl^adows, we did not associate a magnitude with the rays. We now do s®, and^ 
then supplement [H'l with an energy principle or flux principle . We introduce a flux 
density F = 1 Fl as a measure of the energy flow per second through unit £|rea nor- 
mal to a ray; indicating the direction of a ray by a unit vector R , we call FR the flux 
vector, , * ' . . ' 

Kepler in 1604 (by a mixture of mysticism, insight, and soite obseirations of 
light sources) proposed the inverse square law for the flux density associated with a 
source of light. He argued essentially as follows; If a steady source (one not varying 
with time) is emitting rays uniformly in all directions, then the total associated flux 
(total energy per second) passing through any spherical surface centered on the soUrce 
(as in Figibre 4il) is a constant, then, since the surface of a sphere increases as^the 




> 



FIGURE 4-1 



? 1 ' ' 

square of its rac^ius R , the flux density F must be proportional to ^ , ^Visualize 

the source as sometlwhg like a steady 'omnidirectional water faucet/) Equivalently, 



smce 



(1) ^ . ^' /F(R)dS = F(R)/dS = F47rR2 = c,. a constant , ' 

it follows that * . - • . 

Lambert (1760) generalized (1) by taking the jsomponent of the flux vector /r 
normal to a surface gis the measure 9C the energy flow. Thus if S is any surface en- 
closing any steady source, and if N is the outward unit normal on S, then from the 
work of Kepler and Lambert it follows that: 

*[KL1: '-/g FR ^^,NdS ^ /s F cos ^dS = constant = C , ^, *^ 

where 0 is the angle, between the ray. direction and the surfade notihal as in 
Figure 4-2, - 



FIGURE 4-2 ' 



Equation (1) is the special case of [KLJ corresponding to a uniform point' source 
at the center of a sphere; for this case- F depends only on'n , and R is parallel to ^. 
ylf we take the constant in'(l) to equal unity, then the corresp8nding''form of (2) is the 



flux density for a unit point source 
(3) 



Equation (3) corresponds to uniform radiation in three -dimensions. 



We may also apply [KL] to determine flux density for-^a unit source radiating 



; uniformly in or|y t\yo dimensions, i.e., to obtain F(R) for a unit line sourc'e . Thus we 
V considef an extended source along the z-axis emitting rays uniformj^ i^erpendicular 
' ^ xy-plaftes as iif Figure 4-3. We apply [KLJ fpr C^= 1, and S emilfto a coaxial 




y , 



FIGURE 4-a 



right circular cylinder having unit length^aloj^^ Zr and radius R as in Fi^'4-4 




FIGURE 4-4 

The [KLJ integral vanishes over the flat caps of the cylinders at 2 = ±1/2 ; fQ^ these 
pieces, we see that^l^is perpendicular to N ^ iz (wh^6^ is the unit vector and 
-consequently R • N = 'iR • z = 0 . We are thus left with the integre^l^over th6 circu- 
lar wall (oi hei^t unity and radius R), for which v^= R-> R =^ t : , * 

(4) . ^ 'FiR) • jdS = F, • 1 • «27rR =^ 1 , ^ ^ ^ 

where we took C = 1 to correspond jto a unit source. Thus from (4), 

\^ ✓ ' i >i ' * * 

(5) , , " = 2^ ' 
is tne flux density for unit length of unit lin^ouroe. 

vSimilaply a planaV source is defined as an infinite i plane (say zy) emitting 
rays.perpendicularly along = ix = ±^ -^^^^ ^ ^ ^^S^t cylinder 

as in .Figure 4-5, with faces of unit area parallel to the source (and "enclosing** it). 
Since R • N vanishes except^ovtfr these unit fades ^ [KL] for C = 1 gives 

(6) . \ • F/dS = F • 2 = 1 , ' , y • . • 

» * 

and conseqyiently " * 

(7) . ^ ■ • ■• F = i , ; • ^ ■ . 

is the flux denlsity for Unit area of source. ^ ■ 

It should be-kept in minBTthat all the -above Equations are very special cases of 
[KL] . In general FR is a function of all coordinates, and [KTi] applies for any • 
closed surface enclosing any number of steady sources. 
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. ' ' ^ FIGURE 4-5 

From [K-L] It also follows'.that.the integral over a surface S. that does not 
enclose any sources'must vanish:^ . - 

(8) ' 'L FR-NdS = 0 , 

, . J ttfe constant in [KL] is zere for a source-freje r^egion, (The source is outside the 
closed surface, sq, that whatever flows in through part of So flows out through anpther 
^ *' part,) We use ,this to define a pencil of rays (a narrow cone of rays) analytically. 

Consider the capped tubular surface Sq of Figure 4-6 which encloses a set of 




'2 5 



FIGURE 4-6 



rays. The ^curved surface is generated by the rays paaeiiig through the boundary\,','^^ ' 
curve of Si , and the entrance and exit faces Sj and S2 are taken perpendicuiarjta tiie.' ^ 
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rays, i.e., the faces are piepes of the. eikonal surfaces discussed in Section 2. Thus 
if f^c, Ni, and N2 are the normals to S^, Sj , and Sg, respectively, then R '"ti^ == 0 , 
R ^ fJJ = -1% and R . N2 = 1 . Applying (8) to So = S<, + Sj + S2 , we^see that the 
' integral over Sc vanishes and we are left with 

(9) /Sj FdS = FdS , . . 

where the integrals are over the entrance and exit faces of the tube.^ In general F _ 
varies from poiht to point on each face. However, except for special situations, we 
can take the fates small enough so that the variation of F over eacb is ^negligible, 
and approximate (9) by ' *' • 

; (10) ' ; FjSi = F2S2. , . , * 

The set of rays for which (10) holds is defined as a pencil of rays; the set^s ^enclosed 
by a tube whose faces are portions of eikonals-. In cleriving (10), the "$peciar^situa- 
tions which we excluded are those where a face coincides with a focu^ or caustic. As 
discussed in Sect^u 2, a focus corresponds to the intersection of many rays, so that 
' a closely fitting ^tube enclosing such a set would narrow down to ^2 ^ 0; for such cases 

(10) is not a valid relation for F. However, such cases are still covered by (8) pro.- 
vided So does not intersect the caustic. , 

Let us apply (3D) to the essentially two-dimensional problems of reflection from a 
.cylinder with generator along z discussed previously in Section 2. Dropping the un- 
essential z-coordinate (i.e., taking'all pencils as having height along z), we treat 
S of (10) as a small arc length equal to the local radius of curvature (p) times ]the 
small angle (ip) subtended by S at th^ center of curvature (tfie origfn of p), i.§., 

:* (11) . , s = pj^ . 

Since the two faces in (10) are chosen as portions of eikonals (surfaces normal to the 
vayky, and the centers of curvature are the limiting intersection of the common nor7 
— trials to th.e faces, we see that St and So have the same center of curvature,. Thus ^ 

* Si/i)V = S2/P2 = ^> and from (10) we obtain , . • 

* * ' * ' «. 

; = F. I F, ^ . ■. , : ., 

^j^>whifeh, specifies the variation of the flux density ^ith distance along Tays. 

^ We now consider the perfect (co^iplete) reflection <Sl a parallel pepcil of rays of 
"* . . . * • , „ 

' Wjdth'So* and flux density Fq from a^gonv§x cuxvpinealJ-'poftion Cj of a reflector as 

m Figure 4-7. Th6. length of the fei^cqnal of the CQrrespo^nding reflected pencil is Sj at 

Ci, and S2 at a distance 'R from "-CV. Perfeeft,rieflectloA ifjeans that'nq rays peAetrat6 

< ^he reflector, i.e.,^e tS^l incident flux is 690361^^%' tKe'prdbess andtpas^e^ 

* ' through the^erminal cap S2..Thus (tp>hold6; Fq Sq = FiJjk^ = 1^2 Sa' Appro^mat-W, 

^*\x^ the curye^by their tangent lines, we have S^ = Sq (anc^fcli equals^ cp^ a, 



< FIGURE 4-7 

where a is the angle of ^the rays with the normal at Ci); consequently Fi = Fq.. Since 
Ci is convex, S2'> Sj and it follows that F2 < Fj = Fq . To express F2 explicitly in 
terms of Fq , we now use (12); we write the radii of curvature of the ^ikonals as Pi = p 
' and P2 + R, and obtain . \ 

(13) ^,.F, ' ' / ^ 

where R is the distance along> the reflectejd pencil, and p + R is the total distance 
fyom the caustic (the locus of origins of the p's). Thus as discussed in Section 2, the 
^ rays ^Haf pass through S2 appeal* to originate at their virtual intersection point (on the 
caustic)^inside the reflecto^^ 

For the semicircular mirror of radius a . .(see Figure 2,-13), we, found previously 
that P 2 ^ ' where a is the angle of incidence with the surface normal; this 
klso holds for reflection from^ convex portion of a more general surface in terms of 
the radius of curvature a at the^point of incidence*. Thus for a convex point (i,e., a 
point on a convex portion of the^mirror), the reflected flux density equals 

a ^ ' ' 

cos a 

(14) F =.^-^- Fo , 

- cos a + R 



2 



V ** ' Fa 

where we dropped the subscript 2 We may alsdrewrite,the above as P = .l,r\\ 



where Q = - is the curvaturgi.ctf^the^^fional at tl^e 5ceflection point; ..^ , "r 



On the^other hand, for Reflection fi:om a concave point*, the<Jdti^tiQ is fd'aU.anci 
B and p are on the same side of the reflectoif'as in Pigure 4^4-^ For thi 
r.eplace p by -p in (13) and (14), .9J\S obt^ ^ M ^ ^2'^ 



iS'Case-we 
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(15) 



p -R 



^0 



|-cos a 
fcosd^ 



Fo , R ?^ P = 2 ^^^^ " 



where the absolute value is used Wcaiase"we-defined F as a positive quantity^ Ecpa- 
tion (15) specifies F except on the caustic R = p , the special situation (S = 0) excluded 
from the start when we introduced (10). ^ 

Equations (12) to (15) aj^ply for two-dimensional problems. The discussion of the 
corresponding three-dimensional forms can also serve as a vehicle for some addi- 
tional terminology *on properties of surfaces. If the neighborhood of a point on a sur- 
face can be represented by functions alLof whose derivatives exist at the poyit 
{regular at a point ), then there are two, orthogonal directions (the prVncip al directions^ 
on the.surface for which the radii of curvature have maximum and minimupi values ; 
these are the principal radii of curvature pa and Pb- In terms of pa and pb of the 
eikonal at the point of refleb'Son, the analog of (13) is ^ 

PaPb 



(16) 



F = 



(R + Pa) (R + Pt^ 



V 



so that Y is essentially a product of two terms of the form in (13). Equation (16) 
holds if the scatterer is convex at the reflection point; the reflected tube diverges in 
both principal planes (the planes through the principal directions and through^ R). Th^ 
extension of each ray within the reflector lies in general on two caustics; both caustics 
are vi^al, and require negative values of R for their specification, i^e., th^ey^ofre- 

'spond to the vanishing of -R + Pa And of -R + Pb- On the other hand, if the scatterer 
is concave at the reflection point, then pa and pb are negative and the caustics are 

^eal; the reflected tube is initially convergent in both planes. If the refTecting surface 
is convex iii one plane and concave.in the o^jhier (i,e.,^if it has a saddle point), then one 
caustic is real and one virtual; the Reflected, tube is then divergent in one plane, and 

^initially convergent in the other. The first twoS^ses (convex and concave) corre-' 
spond to ellirrtic points of the surface, and the third to a »hvperbolic point s Between 

■ / . ... 



JC 



thes^two classes of points li^s the transition case of a parabolic point : at' a paraboUc 
point, one radius of curvature is infinite, and (16) reduces to (13) as obtained previ- 
ously for two dimensions. If bot^p^ and of (16) |pr if p of (13)) become infinite, 
then we obtain the result, for a plane reflector/^ F Fq . 

In addition to the above points, there are special ))oints at \vhjich the principal 
^;directions are indeterminate: at an umbilical point the radii of curyature of any two 
>ormaI sections {the curves cut from the surface by planes contaiiung the normal) .are 
equaL In^general, the umbilical points of a surface are isolated points. However, we 
Have already implicitly considered one surface all of whose points are umbilics , i.e., 
the plan^ scatterer for which F - Pq , The only other surface having this property is 
that of the sphere.' ' . . 



We can discus^/^^ction of parall^ rays from a* sphere of radius a by exploitmg 
our analogous results for the circular cylinder. One caustic for the sphere is gener- 
ated by rotating that for the circle around the axis of symmetry, i.e., R = - | cos^. 
In addition, therfe is a line caustic R = - ^^^g^ 'which arises from the rotational . 
symmetry of the sphere: all rays, incident on the sphere at an angle a (the rays in a 
^circular tube) giva rise^to reflected rays\vhose extensions iriters,ect the axis of sym- 
metry at the same point (after grazing the epicyolOid) as in Figure 4-9. For"" 





FIGURE .4-9 



the circle in two dimensions we had only two such rays; for the sphere* we ha^rea full 

ring, '^hus in (16^1 the radii of Wrvatui^e p^ and Pb of the eikonal at the reflection 

^^i ... a a . 



point' eqifal ^ cos or and 



2 cos a 



(17) ' . F 



, and we obtain 
,2 



acos 
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^ Jlux and Path Length , A remarkable property of F as in (14) (remarkable only at the 
present stage of our development of a mathenaatical model for Scattering) is that it 
can be rewritten in the form 

Fo a^ cos^ Of • ^ 

,(18) RL" ' 

where L'^ is the second derivative with respect to,^a of the general path length 
I + R(0) = L introduced in Equation {2:3)\ and where the subscript H Indicates that we 
use the condition L' = 0, or 0 =,2Qf, as follows from Hero's principle* (We can re- 
.place Lyj/aZ* by the second derivative of L with respect to arc length along the reflec- 
tor J W6 may assure ourselves that (18) holds by retracing our derivation of the caustic 
in Section 2.3. Our equation g{ot) = 0 for a reflected ray corresponds to L* = 0, and 
our equation g*(Qf) = 0 for the caustic of the reflected rays corresponds to L*j^ = 0 . 

We mention this ^w to make more explicit that, F becomes singular oa the 
caustic LJJ = g^a) = 0 , which indicates a limitation of our present esentially geo- 
metrical model for the propagation of light, and as a preview of a deeper relation 
between flux and path length that must hold for a more co^nplet^ model. 

Partially transparent surface : We can extend the^esent flux considerations for the 
..^ case of a perfect reflector to the case of partia^' transparent media corisidereci in 
Section 3, and'obtain the corresponding reflected and transmitted fluxes when a gencil 
of rays is incident on the curved interface of two different optical media. At the pres- 
ent primitive stage of our model, we simply introduce a reflection factor 0 < ^^^^f 1 
as a multiplier for the values of the reflected flux: (e.g.,.^ in (14)) for the corre^nd- 
ing perfe^ctiy^reflecting surface: Applying (8) to a penci^ of rays incident on a plan^, 
intj^^ce'(witlj^.So "enclosing" the internee as in Figure '4-6), we then find that for 
the incident flux to equal the sum of that reflected and that transmitted we require ^ 
that the geometrically transmitted flux be multiplied by the transmission factor * \ " 

ijp(«): \ ( . ' j I ' If 'j 

I Scattering Applicationg . In the above we applied ( KL] iq obtain tlfe flux density f6r 
elem^tary sources (in one j, two-, a'nd thlpBe -dimensions^, arid to determinCi^the 
ch^e of flux density along a ray in a pencil of varying cjross section. We now extend 
oi^onsiderations of th^ elementary sources to the'analogous scattering problem^, 
. We define^he corresponding '^elementary scatterers** by Jhe previous stipulation that 
the total, radiated flux equal unity^and that it be distributed uniformly over the available 
directions; then we indicate generalizations. We do not solve any scattering problems 
explicitly, but exploit the preVidbs development to introduce terms and general fo^iris 
for subl^equent use. \^ , * ^ / 

• Thus if we have a set oi paraUel rays normally incld^t on a perfectly retledting ^ 
plan^ scatterer at 3? = 0, as in Figu2;e ^rTP^i ^^^n from. S^tidn 2, the incident set o| ^ 
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rays gives rise to a. refliected set of rays and to a shadow forming set of rays. We 
may say that the incideilWkys have "excited'^ the plane and converted it to a source 
of radiation; we call the incident set the primary radiation and the'scalFered set (re- 
fleeted plus shadow forming) the secondary.radiition. and say that the plane has 
^become a secondary soihrce. We define an eleme^tary planar scatterer as a secon- 
dary source fully analogous to the simple ^fehar source considered in Figure 4-5 and- 
Equations (6) arid (7), (In a later section we consider analyticafly the specific jproblem 



of (7) is that the flu^ does not depend^ 



that this corresponds to,) The essential featum 

on distance. Similarly for a more general plamr scatterer we write the scattered 
flux corresponding to the direction of inctdence x as * 

(19) ' . ' M(R,x) , = 

where the direction ol scattering R corresponds eith^er to geometrical reflection. 
R = -X, or to forward scattering R = x, (For parallel rays incident on a perfect 
reflector, it turns out that M is the absolute square of E discussed in Section 2,4; if 
the inciSint flux -density is unity, tlxen M = 1,) 

limilarly if we visualize rays incident pe^endicilarly on a fine cylinder as in 
e 4-ll>, and apply ( H ] essen^alty as for the discussion of edge diffracted rays • 
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|;in Section 2, we ^ee thatlhe scattered set of rays trave} radially outward from the- 
^catterer. We d&iine an elementary line scatterer as a secoiidary source fully 
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analogous to the line source of Figtfres 4-3 and 4-4 and Equations (4) and'(5), i.e., the 
total outgoing flux per unit length of scatterer is unity, aq.d the scattered flux density 
per unit length i? given by (5). Similarly for a more general Une-like or cylindrical 
obstruction, F Js inversely proportional to Rfbut the flux density is no longer the 
same in all directions: * * * * * 

(2or \. . •F= *^ ' 

where'^the direction of observation R may range over all* values in the xy-plane. ^ 

• : • ^ " * ' 

Filially in three dimensions, we visualize a point scatterer excited by rays, and 
> f* • 

define a secondary point source analogous to that of Figure 4-2 and Equations (1) to 

(3). More generally^ for an arbitrary scatterer in three dimensionsj the analog of 

(20)13 . ' . . 

where again M depends only on directions and not on distance. The three functions 
M depend on various parameters, and their determination requires a more complete 
mathematical model than the present one. However, the forms (l§J^(20)r, and (21) < 
give the appropriate,dependence of on R . 

^ — <_Z ^ * 

We ate now in a position to fufth&r-euridiscj/ssion of the relative n^agnitudesr of 
the different rays ot Section 2. Thus parallel rays incident pn a broad finite strip 
'excite essentially two kinds of secondary sources^ the. body of the strip becomes a 
secondary planar source with reflected flux density equal to that incident, and the 
edges become secondary line sources with flux density specified in general by (20).. 
The flux density of the rays geometrically reflected from a plane are irfdependent of 
distance,. |)ut the flux density of the rays diffracted from the*edges decreases in general 

lold on or near a caustic , of edge 

^ __ _ , in order id discuss magnitudes 

near caustics, not only of diffracted rays but of reflected and transmitted rays as 

weiL) \ . M ' - * / r 



t as ^ with increasing R. (The form|(20) does notjjio 
rays; we requi^a morejaomplete scatt#ing m'odel i 



5, Huygens^ Principle > 

In preceding sections we considered the reflection and refraction of a parallel set 
of rays. We started with Eijclidtte restrictive laws [EJ of propagafion ahd reflecticjp, 
and theu replaced [E] by the more.gener^ principle [H*] of Hero that the ray path (L) 
be stationary; using (HM and the calculus we determined the caustics and foci of, the 
rays reflected from'^unred surfaces. Similarly, to consider the set of ray^rans- 
mitted tlji*ou^ an interface df two transparent njedia (media^ specified by different 
indices of refraction /i) we started with the restrictive Snell's law fS] of refraction; 
and then replaced [S] (and [H*]) by Fermat's more general^^^iple [F»J that the jf 
^avel time or optical path (pL) be stgitionary. Thus all our results ontay paths and 
their envelopes (caustics) are c<5Vered by the one princijple [F*]- 

In the discussion' of system3.of rays in Sectioh 2, we alsb introduced a system of 
eikonal curves (eikonal surfaces in three dimensions) that were perpendicular to the 
rays; in Figure 2-16, we sketched ^ome of the ejkonals for reflection of a set of par- 
allel^ rays from a convex cylinder. From the remarks at the end of Section 2.4, we 
see that we can constru^J^ eikonal curve of Figure 2-16 geometrically as the, curve 
traced by the end of a taut string (taut against the epicycloid caustic curve of Fig-^. 
ure 2-12) ^hose oth^ end is fastened at the cusp. Thus if we measure length klong . 
'the string, then each point of an eikonal is at the same distance from the cusp of the 
caustic, r 

In addition to [F'J , we also used the Kepler-Lambert ^ux principle [KLl^to 
associate a magnitude with the rays (the energy per second crossing unit area normal 
to a ray). We used*[KL] td determine'the flux for unit symnietrical sources (three-, 
two-, and one-dimensicmal)^ and to derive the change-in flux density*for a pencil of , 
rays reflected from a culled surface* 

' ' Tlhis 'all our precgding discussion is coveref^/the two "laws of nature"^ [ F| J 
andnKL] plus some of the implicit physiis relevanj to^geometrical optics phenoiien^* 
The baste physics was contained in the two laws, the rest was mathematical manipu- 
lation based on a geometr}^ of rays and some procedu:*es of the^calculi/s. As a pre- 
liminary to the introduction of a:dditional structure inro our mathematical rnQdelj for 
the propagation of light, we now .supplement our previous geometrical construction of 
the eikoSals by an alternative construction called Huygens' principle. This principle < 
by itself does not give us any ;iew results but ^and'this is often much more significant) 
it gives us a new way of thinking about the results we h'ave already obtained, " 

1 , 6 * 

In Section l^we mentioned the two familiar forms in which 'energj; props^ates: 
packaged around particles, or associated with waves. At the present stage of the 



development the fliSc involved in. ["KLJ is iit some sense guided along the geometrical 
rays. It is easy to visualize the rays a& guide lines for very fine'' particles (a view 
held by the ancients, and refined by' Newton*- 1705), hut we may also regard the rays 
as the normals of a system of wave surfaces (the eikonals), / 

Many individuals (Hooke, Euler, ^d others) regard^ light ^s a waVe motion in/^a 
special medium, but it was Huygens (16^90) who introduced the subject as an analytical 
one. His intuitiqn was based on the analogous t\yo -dimensional problem of h9w dis- 
turbances travel on the surface of water, (Drops of water dripping off your fingers 
► above the surface of still water create disturbances at their point of impact that then 
'travel outward in circump ripples along the water surface,) % 




.Huygens used the fact th<t light has a finite velocity of propagation v (as estab- 
lished experimentally by ^omer, 167^) for the development of a wave thebry of light.^ 
He^ssumed that .in a given medium, light starting from an elementary Source at time 
to would Spread as a spherical, surface whose radius r(t) increased in time as v(t K tj). 

If we start a Kght source at time t© and leUve it turned on» the corresponding 
Huygen^s' wave surface is an outgoing spherical front —a discontinuous disturbance 
whose one -dimensional analog.is shown in Figure 5-1, In this figure we plot a magni- 
tude associated with the disturbance (say the flux density F introduced in Section 4, 
or a related function) as a function of time; at tiirie\i > to , the wave front has moved 



a distance v(ti - to) , and it keeps advancing wifh increasing t . (The discontinuous 
Unction drawn in Figure 5-2 is called a. Heaviside pulse') \ 




Starting with an^ advancing wave front (hereafter, the wave surface W) in th»e6 
ditnensions, Huygens regarded pach point 6n the wave surface W as a neW sourQe of 
an elemttta/y^ spl&erical wave (call it a wavelet, Iw) wHo^e^^radius also increased in 
time proportionally to v , Thus'^if the original wave surface W is a sphere of radius^ 



w(t^) 




'7 



, FIGURE 5-2 



r(t,) . the wavelet surface..w. spreads as a sphere of radius R(t) = v(t - t,) the two- 
dxmensiSK^log is shown in Figure 5-2. To obtain the wave surface' of the source's 
advancing wave froht^^gens prescribed ^ ' * 

IHuJ : to construct the wave surface W(t2) at timVtj > .t, . regard 
■ W(t,) at time tf .^as flie lofeus of the centers -. 

. of wavelets w of identical radius R = v(t2'- t,) , .arid tak^ 

W(t2) as the outer envelope of the set of w's. " } j ' 

Figure 5-3.- based on Figure 5.^2. illustrates [HuJ . The esseritial notion is thai; if ' 
we assign an appropriate magnitude function to' a wavelet; then only on the outward, en- 
velope of the set of Ws (i.e.. only o^n W(t,)) do the Magnitudes of the w-s add up 
(reinforce) to give a significant overall effect 



If a pla^iar portion of a wave surface is 'incident on a reflecting-surface, we can ' 
construct the reflected wave frprit by means of [Hu ] as indicated in Figure 5-4. The 




FIGURE 5-4, 



figu/e shows the intudent wave front at tj (plus the two rays or normals that bound it), 
and a dashed front at \ to indicate where the front would 'have reached #time 
t2 > tj in the absence of the refleictor. The actual^ reflected front at time t (the image, 
, shown imbroken> of the dashed front at the envelope of the wavelets generated 

* by tiie incident front as it encountered the reflecting .sur|ace. The dashed fi*ont is also 
the wave front of the shadow forming rays discussed in Section 2, Figure 5-5 shows / 




. Jiow Suj^ens' construction^pr scattering by a strip yields the closed scattered wave , 
su'rfaqe corresponding fo the reilected^lus shadow-forming plus diffract^ ra3!B.of . . 
Figui^e|2-17(d5; the result is pf .coursiJjsimply the cjosed elkonatSf Figure 2-17 (e). ^ ♦ 



.4 i . 



'\ ^ Similamr if th^ scattering surfacei|ls theSnteffac^betw■e^a two differ|nt. optical ^ 
media specined.by velocities v^d'v^we c&istinict#ie tramimi^ed portion's the 



yERic 



V. 



wavelets to taJ^itito account that these t)ortions are traveling at velocity v^ instead 
of V , and then construct their'errvelope lo obtain the refracted wave front as in 
^gui^B 5-6. * * ' . . - . 




^HGURE 5-6 



We have indicated th^ the Huygens* wave' surfaces are^slmply the eikonal sur- 
faces discussed in Section 2, We now apply [HuJ to reflection of a plane wave front 
, (parall4l ray system) by a perfectly reflecting convex semicircle and make this ^ 
k identification explicit. Since all waves in this problem niove with tfie same velocity, 
all distances (L) traveled are proportional'to time (ti, so that we may work with 
either L or t ; in order to exploit our previous figures and results, we work with • 
, distance L. The center of the circular sc^tterer in Figure 5-7 is at x =^ 0, y = 0, 
_ The xqrresponding incident wave is a plane wave front Whose position at any time 



fx-. 




-1 



t = to may be indicated by x =?.xo ; our reference time is t = 0 , and our reference 
position is x = 0 . ^-v. * 

^ We treat the^ Huygens' construction for Figure 5-7 essentially as we did that qf 

, Figure 5-4. We construct wavelets of different i:adii at different points *on the circu- 

lar sgatterqr, the radius, at a point being proportional to the time it wouldTEave taken 

the incident wave front to tray^el: from that point to th§ plane x = Xq . Using'Hu^^ns* 

, principle in this manner we may say that a point tj) = a(a) of the scatterer, where 

5 = -a cos a < Xq, under excitation by the wave front x =<-acosa (se^^Figure 5-;7) 

radiates a ci:jjcular wavelet of radius I ^ - x^ I =^ I a cos a + Xq I ; here ^ = Xq is the 

present position^ the incident wave front. The resultant waye front is the envelope 

■^s^ " - * « 

of all such elementary wavelets. To construct a wave front geometrically, on* draws 
-'■^ . ^ , , 

enough such wavelets to enable their envelope to be sketched. Figure^5-8 shows the 

. case X0 = 0 (i.e., for the time when the incident, front is at the origin) a^nd Figure 2-16 

* ' - ^ a 

of Section 2 shows additional curves for^di^rent values Xg > - g ^ each case, the 

Straight portion of the curve corresponding to the shadow wave front is dlso.tke posi- 
tion that" would have been reached by the incident f r©nt in the absence of the scattei^r. 
The curves of Figure 2-16 can be constructed^either by using the4)resent prorcedufe 
(circles centered on the scatterer) for different constants Xq, or by using the wave 
surface of Fi^re 5-8 as the-locus of circles of identical radii and then drawing their 
outward envelope. , ♦ * • , * ' 




Analytically, we find the envelope of the famil^^ of circles by the same procedur^ > 
,we used in Section % to obtain the ^nvelop^of a set of straight lines. Thus if we take ! 
Xj) = Q , we have the equation of. a Huygens* circlet . , r-' * 

(1) (x + acosa)2.+ (y - asina)2 = a^cos^a , ^ 

The derivative with respect to a gives y = asin^ - xtano;, and eubstttuting this^ 

expression for y into^ the^eguation of the circle (1) gives x lx(l +. tan^a) + 2a cos = 0 , 

'Thus either ' 
> 

(2) . = -2acos3a 



or 



(3) 



y 



X = 0 



.and 
and 



y = asina(l + 2^cos^) 



TT 

'.2 



a sm a 



< (i < 5. 

" 2 




parametric equations (2) describe the envelope 



. so that the corresponding curve (Wj^ of" Figure; 5-8) is half of a two^cusped eificycloid 
(twice the^size and rotated through 90 degrees, as coitipared with that for-Ue ravs 
shown in Figure, 2-12)^ This .portioV^f .the wave' front is generated by a point on ^ ^'-^ 
circle of radius j rolling on. tKe circle of radius*^. The equation^' (3) specify;the 4 
JVg portion of the envelope of FiguVe 5-8, which consists of a line segment of widtb 2a 
nprmal to the direction of incidence; this corresponds, to the shadow forming wave," ^ 

IJavlng one wave surface analytically as in (4) and (3), or graphically as in Pig- ^ 

*ure^-8, we can construct its normals (the rays of Section 2), jLd then obtain any ^ . 

other wave Sfronlj oi^flaying off ^ constant distance al|)ng the -normals and joining the 

points, Wq can construct the evolute.of the wave fronts (the caustic of the iraysj and 

determine that R + ^ fcos a is the radius of curvature, where R is distance alone 

^ \ « - \. * i ^ • 7 ' 1 - , 

the ray from the mirroi:; and, of course, we can ''discovery the law of geometrical re- 

flection by noting that at a given poiht, the reflected a^d^cident wave normals make 

equal and opposite angles with the scatterer*s normal; 

From a "pure"- v^ave view, in order to determine the scattered wave front when 
the/.incidteiit front is aA any position > -a , we use the wave surface of (4) knd (3) 
derived for x© = 0 (W = Wj. + Wg of Figure 5-8) as -thd locus of the centers of 
cir^es of radius l^l ,^and again determine fhe envelope, mechanically or analytically, 
te,, we need not refer back to the surface of the scatterer. If xq > - fj^we obtain the 
wave fronts s^o\vn in Figure 2-16 of Section 2.^ The point' on a wave surface (corre- 
sponding Jto the^cident front at xq) at a distance R = -^a&osc^ *from the mirror along 
a ray, may also -fed designated'liy the cylindrical coordihates r and ^ as intFigure 5-9, 



86- 



^92 



FIGURE 5-9 ■ 



ForAery large v|^lues^of ^, ,\ve see that R and, r are practically parallel 



I • ^1 




0 



that ^ -TT T 2a; we have r « R - acos a. ~ Xq - 2a cos a' f« Xo -2a j s 
which corresponds to a wave from a source at ,x = (the cusp lof the virtual 
'caustic). If we then neglect a « r , we obtain r :Xo and th^^ave fronts approach 
circles centered on the o;rigin of tne mirror,. * ^ ^ . • • • . 

For incidence on the convex semicylinder, these waye fronts are . real in thevs^e 

sense as wfe spdce of real interseptions foi"^he rays of Section 2; for inc^fen^^n tHe* * . . 

'.concave, semicylinder, the wave^ fronts ^f ^Figure 2-16 are virtual. 'I'he vittual ^yave^ 

• / - /' ^ ' 1 

fronts for incidence on the convex cylinder (the real ones for the..concave case/ are 

- ' ; i . ' ' v.; 

obtained for Xq < -a; these are the catves of Figure 5-JL0(b) 6>. (f) plus thqir ln>kgeS . 

in the xAxis.^ If -a < Xq < -f as in:Pigure 5-10(a), then fh^ wave system i&part ieal _ 

(the part near the axis) and part virtiJal. fot the s^i^ere we , obtain, th^virtuarwave 



4 



"fronts by rotating the curves of Fij 



. _^ffi40 arouncj the x-axis. 



In Figure 5-10(b) to (e), the ^traighttiines correspond to the snadow forming wave, 
t^e^9urvBS intersected by* the reflect )r suggest an edge wav^, and^thje remaining^curves 



suggest a 
curv^ c 
ure.5-10( 
figures ii 
wavelets) 
ufe 2-12 



«vav^ outgoing from an origin at x (the geometric set of ^ 

/ i ^ ( * '". "^i^ . a *' ' *^ * 

0 npletely inside the reflector alscJ includes the axial point ^^=?4-o to Fig^-^. ^ 

3a ' ' ' r> -y^'^ ^.Y^y 

) whicti corresponds to xa =<--:r • The most significant ieti^v^ of tne aet of ^ 



from -T- 




that the locus of the cusps (which correspond to. strong r^jppf^feh '€^t of the 
is the virtual caustic of the geometrical rays derived previdusly^i: (Jigr. , ^ 
The figures in reverse order, (f) to (a), illustrate th^t as Xq inls^^^ . 



to -I, the inner 'cusps of the virtual wave system trace the epitjrS^^^^ 
Virtual cpustic of the reflected ray system; f he outer cusj^s correspond to tbe.t^rtual | [ 



.u|^c" (x*= -oo)^of uie shadow iormpg rays. 



means 



we construct the;waVes of Figure 5-10 by the {>rocedure of Figure 6f:7|lE^/^b^ 
of Wavelets of different radii cente'iJwton the scatterer), ^en we fin&'tW;tiie\ 
"edge v^ye'fis generated by wavelets originating relativeljj close to the edge of the 




! ^ 

. . 9,3- : 



\ 1 



scatterer, and that the remii^ing 'curve (the "crater" generated by rotating the curve 
around the x-axis) is the envelope qf the'wavelets originating close to the axis. (All 
wavelets origijiatin^on the arc from the edge to the point (P) on the scatteifer cut by 
the very edge wavelet contribute to, the Inner part of the edge wa\'e, and the "crater 
wave" go r responds to wavelets on the arc between P and the axis.) Similarly i£ we 
construct the curves of Figure' 5- 10 by using wavelets of constant\adii centered oij^e 
wave front of Figure 5-8, then the edge wave arises primarily from wavelets originat-> 
^ing on the, curve near the edge of the scatterer. If we inverf'fhe construction and 
generate the real wave fronts from a virtual front (e.g., Figure/6- 10(d).),, then the* 
wavelets originating inside the rWlector generate all but the straight ]^:^ and ^ 
"rounded corner" parts of Figure 2-16. 




1 im 



\ 



6, Periodic Waves. 



In the preceding section we saw that |Hu] applied to a wav^ front gives directly 
the eikonals of ra^ theoity. fifuygens Represented light essentially as an irre^lar , * 
sequence of isolated disturbances or puls^sT^he essential feature of the mathemati- 
/ cal description of a wave pulse is shown irr Figure 6-1. which represents some 

• ^ X 




■ vt . 



FIGtlR^ 6-1- 



arbitrary disturbance j^ropagating with velocity v alpng fhe x*-axls. The significant 
uspect of Figure G^i. is that the shape of the*p(ilse d^eJ not change iii time. 

If we specify the pulse form^at t = 0 by y = f (x) , then since the pulse form at 
dny^time t is obtained by the translation x to ^ + vt , the pulse form, at time t is 



• . given by 

■ k ■ 



y = f (X vt) 



I 



ults we obtained 



Huygenfe* idea of 



Physically, I we see that thd function f(x - vt) represent^ the unchanged disturbance 
movifig along th^ x-axis (directly x) with constant v^Jpcity v . Similarly a'dist^- 

IbanC4j moving iti the direction -x wouli^ be j'eprelsentedjby f(x +^vt) . ^ " 

% r • * . • . 

By itself JHu] is merely' another method for rederjiving the 
geometrically. , However if we associate th6 idea of periodicity 
waves, then we will have progressed quite far towards .the full m^ematidal mo*el 
we arQ developing, . ♦ • * ' 

Periodicity : Newton (1642-1727), by refractingra pencil of wl)^ li^t through a prism 
of glass, showed tjhat a ray of white lighf could be regarded as m|kde up of rays each ^ 
, having a single color (an idealization c*aHed monochromatic light ), and that the relatwe\ 
index Of refraction ft depended on color. We touphed on (ta^^reviously in .our diWpus- 
sion of the, rainbow >vhin we worked with fi(uj)*, with' u; as^^"c(yior parameter." liis?^' 
studies't)n the colors obtained by illuminating thin transparent plates, essentially 
established that ^ ' ^ ' 




[N] : monochromatic light is periodic with period dependent on a; . 

Newton' s picture of light as a stream of fine particles* subj'ect to periodic "fits" that , 
follo>ved each other at regular mtervals is not appropriate for the visible phenomena 
he was familiar with, but the'idea of periodicity re][ated to color is as significant a6 
Hiiygens' idea of waves. - • 

Young (1801) combined Newton's idea of periodicity with Huygens' idea of waves, 
and regarded monochromatic light as made up of periodic waves. * ^-^ 

' * 

If we rewrite (1) in the form f (p) with - 

. (2r) P = lc(x- vt) , * k =M<^) , ' 

-where k(w) (the "propa*gation constant>')^depends on color, and where p/is called the 
phase of the wave, then Young's principle states 

• [Y] : monochromatic light can be represented by a wave that is a ^ 
eriodic function of the phase p = k(x - vt) . 

Analytically, we express [Y]' as _^ * - ' 

* • " ! 

(3) f(p) = f(27r + p) =: f(27m + p) ; ^ n = 0 , ±1 , i2 , - • - , . : 

. - r. ^ } 

where the period of»f is fixed at 2?: , i.e., f has the same value each time its argument 

* \ \ 

changes by 2n , the excess of p over an iritegral, multiple of 2tt gives position within ^ 

the cy<5le, the, base interval of length 2?:. , • . » ^ ' * I! 

If we add the constraint / * ' 

(4) ^ f(0) = A , • . . ': ; 

where A , the amplitude, is the maximum value of If1 , then the' simplest wave^unction 
satisfying (3) and (4) is the circular function 

(5) f(p) - Acosp = A^cos(k[x - vt] ) = •u(x,t) . 
We may write * ' ' \ - 

•(6) . -•\k.= .^ . 

where X is the wavelength associated with lightof a single color. If we increase x by 

Ax , then \\i increase p' by Ap = 2?: ^ ; each time x changes by th6 length X? we, 

Ax ' • A ' ^ ' * ■ 

^ have =^ = 1 and Ap , 2n , so that f (p) of (5) goes through a maximum and minimum 

in the process. The fac^pl l^x ^ -j- is a convenient dimensionless .measure of d;is- 

tance for a iponochromatic wav^; it gives directly the phase change in units of 27r ' 

.'"^ " ' t • ' 

corresponding to traversing a distance x. Similarly, we may^i^^itef kvt ^ 27r =• with 

(7) kv = f — 
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[N] : monochromatic light is periodic with period dependent on w . 

Newton* s picture of light as a stream of fine particles* subfeot tb periodic "fits" that , 
followed each other at regular intervals is not appropriate for the visible phenomena 
he was familiar with> but the'idea of periodicity related to color is as significant ae 
Hiiygens* idea of waves* - • 

Young (1801) combined Newton^s idea of periodicity with Huygens' idea of waves, 
and regarded monochromatic light as made up of periodic waves. 

If we rewrite (1) in the form f(p) with - 

. (2r) ^ p = k>- vt) , ' k =,'k(<u) , ' 

-where k(^) (the "propa'gation constant^' )iiepends on color, and where p/is called the 
phase of the wave, then Young^s principle states 

. . [YJ : monochromatic, light can be represented by a wave that is a ^ 
riodic function of the phase p = k(x - vt) . 

Analytically, we express [Y]' as • ' 

(3) f(p) = f (27r + p) =: f(27m + p) ; ^ ti = 0 , ±1 , i2 , • - - , 

where the period of'f is fixed at 2n y i.e., f has the same value each time its argument 

# i " ' . X j 

changes by In , the excess of p over an iritegral. multiple of 2tt gives position within 

the cycle, the, base interval of length 27r. , • . * ^ ' 

If we add the constraint . / ' 

9 

(4) ^ f(0) = A , • . , : - 

where A , the amplitude, is the maximum value of If1 , then the' simplest wave^unction 
satisfying (3) and (4) is the ciroul^r function 

(5) f(p) - Acosp = A^cos(k[x - vt]) =.u(x,t) . 

We may write ••. ' ' *. . , 

•(6) ■ . •\k.= .f . . . .. 

where X is the wavelength associated with light.of a single color. If we increase x by 

Ax , then \\i increase p' by Ap = 27r ^ ; each time x changes by th^^ length A? we 

Ay ' . A ' ^ ' • . 

^ have ^ - 1 and Ap , 2t: » so that f(p) of (5) goes through a maximum and minimum 

* 27rx ' . " 'i 

in the process. The fac^p! I^x ^ is a convenient dimensionless .measure of dis- 
tance for a iponochromatic wave; it gives directly the phase change in units of 

" ' t ' * 

corresponding to traversing a distance x* Similarly, we niay^^^^'ite kvt = 27r =• with 
(7) ' . , kv = f — 
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r 



as a dimensionless measute of time cbrrespohding fo the pHase cjiange in units of 2ir 



for a time interval, t . From (6) and 
(8) • . - 



i 




) we write v/^e phase velocity, as 



We are now in a position to interpret ;;he paraniefiers \ and T ^ X/v introduced in the 

^ i ! « 

above as ^^ell as the corresponding ''dolor parameter" we have mentioned previously. 

In Figure 6-2, w^ plot u of (5) versus x for fixed t = to\ and ip Figure 6-3, we 




FIGURE 6-3 



plot u of (5) versus t for fixed x = . We see that at a fixed point in/ time, U is 
periodic in x; the wave form, is repeated at intervals of length. X, which is why \ is 
called the Wavelength or the space period. 'Similarly, at a fixed posit/on on the x-axis, 
we observe that u is periodic in t: the wave form is repeatjed at tune intervaj^s T 
^ called {he time period or simply the perio^i. 'I ' 

At 2k given time to, we obtain the wave form of Figure 8-2. An instant later, 
to At , we <^btain the same system of crests and valleys ^ith ea^h point of the wave 



' • /x + AV to + At\ X tft ' <^ • 

shifted Ao x t Ax, such that ( — ) = ^ ' j * Thus we may visualize 

the wave as traveling in^he direction x = r- . 
' ^ -» A 

' The reciprocal of T is called the frequently (i ) of the source producing- the wave 

H f * 27r 

and it is convenient to measure this frequency in units. of 2n , i.e., to use w = -^jr 
where ^ is called the angular frequency^ Thus v/h rewrite (5) as 

(9) * . ' u = Acos(kx - wt) , ♦ 

and we identify the color paramete'r as the angular frequency of the wave associated 
with light of a single color. » 

Th6 angular frequency is a fixed characteristic of the source of the waves, and 
does not depend on the optical properties of the 'different media (characterized by dif- 
ferent v) through^hich a wave p^sses^, however, the wavelength % ^ 2;rv/^ do^s * 
depend on the medium. In ge;.neral the phase velocity v is a function of so that 
wave.s of different frequencies travel with different velocities v(u.) in the same mate- 
rial. Equivalently, since the indejCof r.efraction is defined as ^inversely, proportional 
to V , w6 may rephrase the above in terms of m(^) • T^kifhg the development to equa- 
tion (9) as applying to a medium with indgx pf^ refraction m ^ 1 (free spacje or „ 
'vacuum), we replace lo?*^acJ;he more general case by . * 



(10)' / • »kMx ^'I^E^ii ="|^x',;, \, = \ 



where* \a. function of ^ and the material) is. the wavelength in the optical medium 
defined by ^(o;) . - * * • * ' 



The wavelength X is the scale factor we anticipated \vhten we spught to distin- 
* '** * ' J- * ' -* V ^ 

'guish bfitweeji geometrical reflection ^hd edge diffraction at the^end of Section 2. ^We 

could have introduced much of the above structure into the ray picture by associating 

the idea of phase (periodicu'ty) with a ray. However, the wave jiicture is in general 

moi^e fruitful forj^he usual, visible phenomena. Fox, convenience jin the following, we 

may tise a mixed terminology^ with the rays^^ndejrstood as the corresponding wave 

normals. ' ' ' , < 

» ' ' ' . 

Thus we say that if ^ight of a single color travels a distahce L in free spaced its 

phase has changed by^.lcL. Cprresponding to the unit sources of Sectiop 4, the phase . 

at a distance R along a ray fr9m tlte source "differs from thfephase at the source by^ 

kR . Similarly for the*vefle.ction pi*Qblems of Section 2, the .phase of fhe reflected ray . 

at P* in Fig3ure,2-4a differs Jby k(Li + I^) from the phase at S , and the phase of the 

ray at in Figure 2-7 relative to its phase^t x = 0 4s giveiiby' ^ ' : . 

^ V . ' ^ p'= m + R)* =^ k(|, +• yfi^ ^X)2 ^ (r]- y)2) . / \ . ^ ^ ' * 

From the geometrical methods of constructing an eikqn^l (wave surfarce), we.see* that 
it i^ a curve (or surface) of Constat phase (i.e., there is no phase difference between 



any t\voi)omts on an eikonal), and may label a particular eikonal by a particular 
value of p. Similarly for the elementary sources of Section 4, e.g.^ we may su^^j^nd 
the point source of Figure 4-1 with a set of spherical surfaces of particular pkdii R„ 
corresponding to the particular phase differences kR^ 

Instead of working direQt^ with cosp i{ is mgre convenient to carry out manip- 
ulations with 



(11) 



cosp + i sinp , 



and take the real part Re(e*P) = cosp when we want to exhibit the periodic behavior 
explicitly. If we represent ^li) in the comple'x plane we obtain the vector diagram 
(Argand 'diagram) of Figure '6-4. As we progress along a ray (increase L)| p - 




FIGURE 6-4 

increases and the tip of the vector of unit length describes a circle of urvit radius. ' 
The projection of the tip. on the x-axis (the real axis) is the oscillatory function 
V cosp; each time p increases'by 2-n (each time the tip describes a full cijcle), the 
x-projection goes thaough its maximum (+1) and minimum (-l)<^alues. (The function^ 
e*P is oft^n called a phasor .) More generally, we work with . . 

(12) ' ' . ^ . . ^.Ur Ae*P , ' 
where the* amplitude^ A is posi^e. * ^ » 

In subsequent application^ we use the exponential form )^ 

(13) • ■ . . K u = Ae*(^'^-^) ' 

such that u of (9) corresponds to ReU^ J^^S'^^SlSB)^^^^^?^^^^ \vave traveling' 
in the x 'direction. * ' * ^ 

Tp tie in the present discussion with energy flux consideration§«i|-Se*ction 4, we" 
,4ioie that (in general) at distances from the source large compar^^Pravelength we "* 
may approximate A by a constant times /F , where F is the f^^^nsity introduced' 
fot the Kepler-Lambert principle. We write 
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(14) ^ ' « U « C>/F;eMkL-a,t) ^ 

where L equals x or r , and where F in general depends outdistance. . 

For the point source (or point seatterer) at the origin in three dimensions, we 
showed m Section 4 that F . We therefore write the corresponding wave as. ' 

(15) ; U^C3^-j^, r'= /x2 +y2 + . 

/Where we used the dim^sionless kr (instead of r) in the denominator. Similarly for a 
hne source (or line seatterer) along the z-axis, the wave corresponding to the flux 
density F = ^ is 



Qi[kr-L.t) J— 

(16)0 U ^ C2 ^ , r = Vx2 



+ 



-2 



In the same sense that we interpret (13) as a wave traveling along x , we speak of (15) 
and (16) as waves traveling outwardly along r - , or as outgoing wave§. For the 
planar sourte (or planar seatterer) at x ^ 0 (one-dimensional case), It is independent* 
of ^istapce; the analog of (15)^ and (16) is 



U= iCe'(^-<.) for « > 0 - •• , ^ - 

t ■ ) Ce-Mkx-^u-t) for X < 0 , 

w)iich we rewrite compactly as 

(J8j . U = CeMk'^'T'^O .A ' 

explicit dependence of (t5) and (16) on'f and t facilitates qualifying [Y], 
which holds rigorously only for a plane wave, more generally, the waves are periodic 
m t , and only approximately periodic i» r (or in p) because the denominators in (15) 

and (16) are not periodic. Thus an increase in r corresponds not only to an increase 

~ - cos^ki***wt^ ** 
in phase buTajlso to a "decreaseTin magnitude, e.g., so that although — ^Tias ^ 

maxima at^ual space inten^als X, the magnitudes of these maxima decrease with in-\ 

, CBeasihg,^ However, such magnitude effectfe for k >> 1 are insignificant for the 

problems at hand. - *. * , <>• 

^ ■ '-^ ■ ' . • ^ , , " \ . 

We should arso qualify the preceding .discussion of phase for the reflection problem 
by explicitly restricting it to perfect reflection. If the sca^ef er is partially transpar- 
ent^ then the raj' reflected at an' angle a. with the surface, normal undergoes in general 
an additional phase chahge ^(o^), which we add tothe,abo%e p^ . ' ^ 

interference : Xhe-ioncept of interference was introduced into .wave physics by 
Yoiing. mJ^'c cHscuss interference subsequently in detail but tnenti<)n it now to stress ^ » 
the mdst significant feature arising from associating a wave (or mor^ specifically a 
phase) with light. Tl^e essentials are indicated in Figure 6-5 for scattering of a 




monochromatic plane w^ve by a screen containing two very narrow, slits separated by 
d » \ . the waves (or rays with phase) thdt arrive at x,y (x » d) from the two 
sliis have traveled different paths Lj and L2 , and thetefore differ in phas^ by ' ' . 



(19) 



<P = k(L2.7.Li) « kd^in^' = ~^sin^ 



Their magnitudes differ little,' and .we may write the resultant wave at ^,y. in the form 

(20) , • • ^ = Uj + U2 W(l + e»^", ' * 

where we have absorbed e*^^i"**^ and other factors into W.. The corresponding eflergy 
flux density is pro|)or^onal to ' ' * . . 

(21) 'F = IU|2*= lw|2|i -he*^|2 = lw|22(l + cos^)'.'^ ' * 



'^Sfetting IWl^^equal to unity, ^ve show the^ssentiMs correspbrtding t6' 
(S^) ♦ i f'= ll + e**^'= 2(1 +-cos^) ' \ 

vect^rially in Figur6*6-6. 




FIGURE. 6-6' 



" ' t° ' that' if ^£>'/{^ (i.e., ,6^0 along the x-'axis) then F =4. (This corresponds 

■ ■.' esseritiali^to a caustic of edge rj^ys as discussed in Section 2; however, we now have 
, much^ore structure for the descriptio1ii*of light -in th'e^hadow region*) 'As we vary 



' and the angle of observation ^ ^ ^ , the intensity F goes through a maximum val^ 



of j4 .when 

(23) , = 2n7r, ;i = 0, ±1,'±2 

: ' 

and a minimum of zero when 

(24) ' . <^ = .(2n +i)7r . 



This behavior is clear from (22), and more graphically from Figure 6-6: iijp = 2n7r , 
theivthe vectors of unit lengjjji point in-the same direction along a straight ^ne and , . 
their resultant is 2 ; if </? ="(2ny+ 1) tt, then they point in'bpposite directions and ^ 
^cancel each other. The results F>witK variation of <P ate shown in Figufe 




FIGURE 6-7 



Thus- for a monochromatic wave, (fitted X), a parallel screen on the shadow side 
■ of the slit-screen (flie dashfed line at x in Figure 6-5) show -bright arid dark bands : 
the bright ^'fringes" corresponding to ^ ^ 2n7r are located on a' screen at ^istaijpe 
X from the strip by , ^ . ' • V • 



(25) > 



dsina « ^ ='0, X, 2X; . = nX , 



nXx 



i.e., 'when the 'path diffbrence i^ an integral number of wav^engths. Similarly the 
dark fringe's corresponding Xh <p ^ (2n + 1») tt are located by ' ^ 

We call (25) "construct!^ interference," and (26) "destructive interference." 

If we use white light (a mixture of tvav^^ of different X*s), theri along^the axis 
y = 0 , we obtairf^a whiie central fringe; ^however, from (25), the side fringes are dis- 
placed along, y directly 'in proportion to a pidwe tlfftefore see'bands of different 
colors* (Analogous phenomena in the shadow region of a wide strip were first noticed^ 
by Grimaldi.) jgomparing (25) with measurements we find,irom the displacement of ^ 
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the banUs of light of different colors, that the wavelength for red light is about twice 
that of blue light, i.e., • ' / • 

and that the cglors orange through yellow threugh greeiwh'ave wavelengths :\ of length 
intermediate to that of red and bfu^> • • 

' In the following sections we consider several elementary applications to scatter - 
ing'phenomena of ihe Huygens-^^Tewton- Young periodic wave theory of light. Thesej?^ 
applications are^kssociated with Fraunhofer (1787-1826, an experimentalist), Fresnel 
:^j[i78&-.1^27, a theoretician), and Rayleigh (1842-1919, bothT. . . ' 

K F3:aunJxofer Diffraction bv a Slit. >e'now apply the wave model to Fraunhofer diffrac- 

♦ ^ ' , * ' — 

tion of a plaijf wave by a slit of width, 2a in a perfectly reflecting plane ^s in Figure 

6-B: We»take the origin kt the cehter of the slit. 




-a 



FIGURE 6-8 
1^- 



1 We write the incident wave as * ' 

• . ■ ■ • • 'V - ■• 

(28)^ _ . ^„ . U, => e'('=^-'^«). ■ . ; 

and interpret (28) as a wave o^ unit flux dfensit^f traveling in the x direction. .Using 
(^Hu J ^implicitly, we regard U, as exciting wavelets in the jVf^e of the aperture, aYid" 
specify^a^vavelet originating at x = 0, 0 by the elementary <5utgoijig wave form 
^ 'yfkr ^^^^^ Similarly for a wavelet origiiaating at 0, tj as in Figure 6-8 

we use " - • " . . 



(29) 



eJ(kR-a;t) . . [ 

u ^;C^ ^ . , - R= ^'^2 -. 2rsin-e + 7j2_ 



ERIC 



98 



104 



■ . . . ^ ' 7 

Every point of the line y = 7],x = 0,for-a:2 7]^a (every line element of the slit)* 
corresponds to a wavelet of the form (29), We represent the net effect of all such 
wavelets at a distant point/ f by the integral 

/.a 

(30) 



U = L U(7])d7] 



Restricting'consideration to r » a ^ we 'apprpximate R in the exponent by 
(31) - / R » r - T]siTj.e\ 

In the denominator we use simply R ^ r , because lul Is much less sensitive to 
changes in the denominator than to changes of the phase, (Fropa Figure 6-6 we see 
that a slight change qf the magnitudes of the two vectors has little effect compared to 
a comparable changbsin the phase difference (p ,) Thus (30) reduces to 



(32) 



pl(l<r-a;t) - 



(33) • 



d "a 



i.e., Uis an^lementary cylindrical wave (source at the origin) as in (25), times 
function of angles G(6). (the scattering amplitude ). Thus for 0 ^ 0 or tt, we have 

^a 

(3^4) 



where 2a .is the 
'and obtain 



, . G(0) = G(7r) = jf^ dm, = 2a \ 
^idth of 



the strip. For the other angles, we^integrat^ the exponential 



(35) G{e) 



gika sl^d "Q^^^ sin 6 

iksin^ 



where S is the width of the strip,^nd T |s an oscillatory fuijction witli zeros at 
kasina*= n?:, n = ±1, ±2,. . . . The limit of 1(0) for ^ ^ (J is 1 , . - ^ 

Rayleigh-Born scattering by a/sphere . As another illustration ;ye consider Raylefjgh-^ 
Born scattering by a sphere of radius a \yhose optical properties, differ only very 
slightly from the ftee space in .which it is imbedded (a **tenuous" scatterer), W^use 
ihe geometry of Figure 6-9 with the sphere at the origin of the coordinates, and take 




FIGURE 6-9 
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the plane wave 

(36) ' , . Ui = e»(»^'-^0 

• as the incident field. ^ , 

^ > ... 
We regard the sphere as made up of elementary sources of spherical wavelets, 

, , ~ . ^ p i (kr - CJt) • 

such that the source at ttje origin produces a wavelet C ^— r as in (15). llie 

^ elementary sourxje at the position P (|, t], ^) excited by the incident field, 
produces an effect at r = (x, y, z) described*by ^ 



(37) 



R - >/(x - ^)2 + (y - 7?)2 + (z-02 = V (r2 -f p2) -2r • iS , 



where the phase is chosen to agree with that of .U^ when R = 0 . The net effect Is 
represented by the volume integral 6f u over the sphere of radius u: 

(38) - * ' U =• /u(p)dV(p) ^ ' • ' 

where dV is the V9lume element: ^ v , 
\ » 

Essenfially as for\the slit problem, we restrict consideration to r » a , and 
* approximate R in the exponent of u *15y ) ' 

(39) . V R « r - ? .-jo , ^= i- . 

where f is a unit vector ;\in {he denominator, we use simply R « r .• thus (39) ■ 
becomes 

; ^ \ - • * \; 

(40) ' \ ' TT J ^ e»(J^r-U) 



U-C^— -G(.), 

™. ^^^J^' ^ \ 

- ^ \GeO) = je»»^(S-p- r)dv , 



i.e., U is an^ elementary spherical wave outgoing from the origin as-in-(lb), tirpes a * 
scattering, amplitude Q{0) that is independent of r butMepends on the angle of obser- 
vation e (and, because of^ymmetry, on no other angle). 

Since t • z , where z = - is a.unit vector, we may rewrite the scattering 
amplitude as ' , ' v * V 

#2), ' * . G{0) = ie<->-?)dV . , , * r . ^ 

If r = z , i.e. , in the forward scattered direction 6 = 0 , we h^ve G(0) = JdV = V , whejre 
V ns the volume of the scatterer. For ^Itrary the easiest way to integrate over 
^ the spher^e is to introduce a'new polar axi^ in the direction z -,f, 



From Figure 6-10,we see that since ^ and r have unit lengt}i, z - r is a 




/ 
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vector of length 2 sin ^ . Thus we may write 



m 



(43) 



2 - f = 2 sin y Zo 



\vher6 Zq is ar'unit vector along the, axis Zq oi a new coordinate^ system. Consequently 



(44) , ' ' - ; '(z - r) p = p . Zo2sin| = Zo2.sin| 
\jhere Zq^^=^p * Zq i& the projection of p along>the new pol£|^r axis zj . 



We now have 



(45) 



q((0) = Je ^ = Je ®dV , K - 2ksin| . 



i2Kz 



We introduce the polar coordinate system of Figure 6-»ll, so that we may replace 



— I 
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' zo by p cos^ and use dV - dap^dp sin/3cl^ with p ranging f?om 0 to* a, a from 

. - r^TT .a i .TT V ' ' 

(46) G(e)-= J daf e%i<«'8sin^d^. ' 

The integral over a yields 2rr airectly. Th^ integral over p maybe rewrit^ as^ 



" ' . . 0 , ' -iKp ' , Kp ■ 

■Thus, we have reduced (46) to - ' ' ' : 

" • ' ' a ' '-a ' • ■' ■ 

♦(48)' G(0;=2j^V^%^^d^^^'psinKpdJ. 

\ , " • ■ ^ . 

, The^emaihing integral is simply the deriy^^ive of /pos Kpdp with respect" to' K : 

i: 

siu^l 



V 



which we rewrite" as w . - • ' * ' 

- (50) G(0) = V3[^ - ^ ^ VJC9) =.:G(0)jp-.'' V = x'h 2kasin^^ 

Where V is the voluml^^of the sphere. The value ^0) = V obtained from'(42) for 
• ' ^ r = z , e = 0 , is also the Umit of ^50) for X — 0.:^ Wore^generally, the result 

' ''(51). . , ■ ;■' / G(o) = /dv = y', 



• which follows from (42) for f = z (i.e!, for forwdfd scattering), holds for a scattei^l^^' 
of arbitrary shape subject to\the present restriction tha^it be "tenuous"^^-^ - , >^2^^ 



From (40)rand (50) we hive 



(52) , , ^ -u =,c;^-^,vj(el) 



^ ^ SO that the c'orresponding^flux density is.proportiQn^^k) V^. If a «.\ (ite., for a 
"small scatterer"), we have J(0) 1 , and consequently 

• (53) / U ^ C ^ V \/,>. - ^ . ' r 

In^the next^section, xve^ohpider ail applicatioii of scattering which hinges essentially 
* on the fact that U oc v/r for a variety of different soatterers. p 



7. ^ Rayleigh Spattering. 

.In the pre^nt section, we digress from the development of the mathematical , 

i ' ^ 

model for scattering in order to discuss a beautiful application to nature, this mate- 
rial serves as a supplement to Sections 4 and 6, and provides a more complete discus- 
sion of one of the topics touched on in Chapter 2. The development is based essentially 
on the flux principle [KL] , an^d on the fact (determined by interference experiments, 
a^ discussed previoirsly) that sunlight mayHDC decomposed into light of different 
colors fVom red to blue wUth associated wavelengths \ to Xb such that (approximately) 

(1) ; ; „ -X, = 2Xb . . - ^ •* * 

The wavelengths X associ^t^ with the intermediate colors of the vi'silJle spec£rum 

(prange, yellow, green) satisfy X^ > X > Xj, . 

* • • u • 

In 1871, Rayleigh developed a mathematical model to account for the blue color 

of the sky, and^for the red color of clouds>near^usk. The essential feature of 

Rayleigh's model is that when rays of different colors (different wavelength X) are ^ 

scattered by the molecules of the earth's atmosphere (mainly nitrogen and oxygen), 

th^ scatt^i^ers may be re^^rded as seconidary sourcies in the sense^of Equation (4:21,)'^ " 

with flujF^S^nsLty Inversely proportional to X^ , i.e., ^ 



where r is distance from th^^catterer, and C is independent of r and X. Rayleigh 

** i 
used a more complete form of Equation (6.38) (one we consider m the last section of 

this^lbtpter) to show l^at U — - iot various scatterers^ (e.g. , broacf range $ of rela- , 

> " rX^ * ^ ■ * 

^tive index of refraction m) provided their dime*hsions were small compared to X, and 

alj^ c6nstructed a simplified intuitive derivation of (2). » ^ ^ 4f 

The .simplified derivatioa^is l>ase.d on dimensional analykis. Thus, if w6 divide 
the scattered fjji^densitj^ ^ 1 Ul? = incident value Fq , then the result is 

ind^enderitxjLDiaj^ltiLin which w^ the ratio ^ » 

(3) -^i^-.. ^« ^ ^. — r - — X ■ 

is simply a nunoher andjdoes-aot depend on units or dim^sions, (since F and F^ are 
simply different values of the same physical quantity). For a scatterer of volume V 
whose dimensions are very small com'pared to X, Rayleigh assumed, 

^ [Rl: . ■ G=.BV. ' f = ^ 

where B is independent of the length dimensions of the scatterjer. (In Equation (6:53) 
we considered a specif case of [R], in Section 9. we consider an exception.) He could 
then obtain (2) from (R] ^by using dinaensional analysis. 



We incficate the length dimeirsion involved in [R] as {l} and consider the powers 
^of* L that enter the various terms. We have * > 

,(4) • _r2*={L2} V2=:-{L^^ 'f ' ^ ' 

whi^h state merely that since r is a length, r^ must be the square of a length; since 

a volume is a ifength cubed, Tnvolves the sixth power ©f a, l^engtb/vThe fact that the 

f7 isr^ii^ensipnless (a "pure*' nuntber) means that B of [ilj' must satisfy. . 

I ♦ 



~<5) 



B2 {h-*} . 



In view of the irestfiction on B mentioned for [g] , ^he only length pakmeter avail- 
able for (5) is \ . Thus, follows that ' ; , ' / 



(6) 

... ..... , ^ / 

where g may depend on optical properties (through a relative index bf refraction u) 
andjdirections. ' - * > . 



From (6) and [R] we have Rayleigh^s inverse "fourth-power law" 



Substituting (1),, we see that ^ > • . 

Thus, if a beam of~whH^;;i^t is incident on a small scatterer as in Figure.^7-1, the" 



0<' 




I 



• . • FIGURE 7-l» • 

% ■ . ■ " . ; ■ » _ ' / 

blue component Qf jvhite light is scattered sixteen timjes as strongly as the red, i^, 

the flux ratio of the blue and red components olJserved at an angle Q from the direc- 

^ • * ^ * 

tion of incidence is given ^yj(|^. 



Equation (7) specifies the scattered flux density at a point v{Q) as in Figure 7-1.' 
The flux scatter^dlthrough a spherical cap on a cone of half -angle a around the for- 
ward direction 0 = 0 , as m Figure 7-2, is obtained by integration, ^he area -of the 
cap, the portion of the sphere of radius r , is given by r^ tiipes the solid angle 
the cap subtends at the origin (i.el, r^Q is the surface analog of arc length vO 
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encounteifed in earlier^^^diroensionaL^oroble^ * . 

measured from the' direction of incidencye z , anji in terms of ihe angle tp (the' azimuthal 
angle)^measured in a plane normal to z , the area of the cap is gvven by 

r^J dQ ^ r2J d^l sini^d'0 = 2Trv^ (1 - cos a) . 
The corresponding flux through the cap may^be wriften^ 

M 

(9) 




Substituting (7). into (9) we 'obtain - ... 

(10) a(a)='^r ^.dfl . - * 

' so that <r is ^independent of distance If a = tt, theri th,e c'ap becomes a complete 
^ sphere; the value a(Tr) (called the "total scattering cross -section^') is the total flux 
scattered in response td^ incident plane wave of unit f>ux density. If the incident 
flux density is (say) Ij , then the total scattered flux is Ii^Tt) . ^ ' . ^ ^ 

• I^t us now visualize' an incident beam of rays flowing through a tube of cross- 
.sSctiorfal area S as in Figure 7-3; and' apply tlie 'Kepler- Lambert fliix principle [KL}. 
The flux through S at Zi is Ii S, anfi from [KL] :this equals the fliix IgS through the 





. S ' ^ 






• 
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surlace S at Z2 plus the flux through the incomplete spherical surface consisting of 
the sphere minus the cap os.^ i.e., 

where the integral is over the incomplete spherical surface. Thus, from [KL], 
J12) . I,S,= I,Q(a) + IjS. 

Regrouping terms ^we have 

(13) ' ' I2S = I,S - I,Q(aV ^ ^ J 

which states simply that the flux in the beam L S at Z2 equals the initial value 

IjS minus that diverted to" other directions- by scattering out of, the beam. , 

* Jrbm {91 (10), and UiXjWe have 



(14) 



which we rewrite as 
(15)/^ 



Q(a) .= m . '^("l- S-L.' ^ dfi . 



in order to emphasize the dependence on A, If we neglect ,o"(a) dnd, approximate Q. 
^by the total scattering cross section cr(7r) (i.e., if we ignore the "hole" in the spne^pal 
surface)^ then the difference 'between the initial and final values of the flux afong a ^ 
parallel ^eam intercepted by^ scatter er anywhere along the beam is 

a thelfS We^ N ^lich i^catterer^ in fiie bearn^* theiuindjer appropriate restrictions we * 
n^avusp (16) with tlie i:igbt-h^' side multi£>lied by' N.j09appro5cimate4he net ^ff-ectr? 



(17) 



'If there are n scatt'er^rs in unit volume in the geometry oi Figure 7-4,.cthen we have 



» * ' : i_: ±12 ij-t. 



"t^. Is.-..; 



Zs 
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N'= nS(Z2^ Zj) , and^(17) reduces to/- 



^<18) 



I2 - Ii = rnQIi(Z2 -Zj)\ 



In the limit as zj j^pproadnes Z2 , we obtain 



(19) 



and consequently 



/ 

/ 



(20) 



with Iq as the vaWe at say i^-p . ^ 



From (36), ^e have .Q'= /so that 



-n(AA*)z 



^21) 
4 

which is^the form considereji in Chapter 2. Consequ^tly I(X) decr,eases rapidly as X 
decrease_s,^,and white li^t becomes tprogre^sively with distance) relatively stronger 
in the red -yellow^ corrjpoiients bepause it is losing itsiblue components via scattering 
to other directions. Thus, as (7) accounts for the blue! of the sky in directions other, 
than towards the sun, (21) account^ for the. re^in^Sis. of the clouds. illuminated l^y sunlight 
near dusk. See the discussioyi of (6)^ in Section 4 of Chapter 2; 

' ' Equation (i) gives the scattered intensity for one scatterer excited by a wave' of. 
unit intensity. If one scatterer, of a collection is.a^ a diplance z |rom ihe entrance 
faqe of the region of sdatterers (as in JFigure 7^5), then we-multiply (7) by 



» • • • f • 

> > / ^ 



• 0 



) FIGURE '7-5 



> ^ 



*I(z) - Ioe""Q' of (20) to account for the intensity^ Idss^ of "the excitation that reaches it. 
Similarly, if we observe a scattered ray from this scatterer at a distance r from its' 



•s ^ ^ y ^ < 
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center, we incorporate an a,dkitional factor e*"Q''Uo account for the additional loss. 
Thus, the scattered intensilij^ for one scatterer as in, Figure 7-5, becQmes ' 



• (22) 



r^X* 



where z + r is the total ray path within the region of scatterers..' 
. Let us rewrite (22) for z + r = 1 as 



(23) 



e- 17/ A. 



"a form that shows that F vanishes for both \ ~ 0 ani X ~ " , and has a maximum at 
a definite value of \, say X = A. Differentiating (23) with respect to \ , we obtain 



(24), 



dF ^ _d|l dv;*-^. Be-D^-^a - x-<nv i' ■±\ 



wMch vanishes for the wavelength 
(25)/ 



corresponding to a maximum scattered intensity 



(261 



(27 



F •= ^e-i 



B 
A*e 



Dividing (23) by, .(26) we write the scattered intensity as 



1 / fA^-f ,V \ 



SO that F is expressedffin terms-' t>f the maximum valie. F^^ and the corresponding 
wavelength A , This simple model ap|ilied to skylight gives a maximum wavelength 
in the'Mue-gre^regi^n <x|^he'«unt^^ — ^....-^ ^^""-^-^f* 



I 



8. Method .ol Stationary Phase; ^ 

. * >Ve5jare tiQwJA a pQsitiQR t€k:bridge the gap between the geometrical optics ray 
procedures of the early sections and the elementary wave procedures of Section 6. 



/ 



We do so initially by considering diffraction by a strip. We work with the wave forms 



of Equations (28), (29), and (30) of, Section 6, and the geometry of Kigure 8-1. Thus 




x,y ^ 



' i. * 



f;gure 8-1 



we take thte "Incident wave proportional to 



(1) 



Ui 



the wavelet from a secpndary line source oj/the strip as 



' (2) 



ce^ 



R = \/x2 + (y - 7))2 , 



a^d represent the net effect of the wavelets at r as the integral 



(3) 



.u(7?)d7? = c J 77= dr} 



-a /kR 



' We have suppressed the time factor ^e"*^ for brevity. The actual wave forms are 
obtained by multiplying ^he above by e"*'^'* and then taking the real part of the result. 

, .^Xhe situation m Figure 8-1 is an^ogous to that of Figure 6-8, and if we 
assume R » a ^.1 I max obtain ^the same forms as (32) to (35) of Section 6. Thus 
"^if we expand — to first. order in (a$, previbusly) , we obtain R « r - 97 sin0, and 



consequently the j^revious procedure«yields 



J ce 



ikx' 



y/kr - 



..a 



ce* 



\/kr 



ce 



ik sin 0 



y/kr 



2a 



sin (ka sing) 
ka sin 0 



ce 



ikr 



2ar(0), 



Thus'except that the„constant c may^'differ iathe two .cases, the present result (the* > 
Fraunhofer approximation fof scatte'ring or diffraction by a strip) is the%ame as tfiat ' . 
for the slit. (The relation of the ;result for the strip to that of the slit shown by'the 
abo\^ is a special case of what is called Babinet^s principle.) 

We, introduce y = kasm0 to represent the phase difference between the rays 

from the cente!: and from a^i^edge to the observatioji point, and write F = (the 

■^Tragnhofer, pattern factof»'). The principal maxima of rxorrespond^to y = 0 (i.e., 

to the fonvard and6ack diFectionsT^^rrdnd- ^ respectively). The secondary max-- 

ima of r occur ;4rlhe^efo&ibf Ian ^'y, which are given^approximately by y ^ 'l.437r, 

2.467r, 3.47^ and, for larger vaiues,-by y « (n + ^j'i] for^-the first three of these 

ze^os of ^ , we have £j^ -^^^^ respeotiveijj. The zeros of F corre- \ 

spend to y = n;:. The angular half -width of the principal maximum (obtained from 

the position of the first zero, y = Rasing tt) is'" sift ^? ^ o =^JL = Jl 

, : • ka 2a • " 

The form (4) is i/estricte^^s^^- » a. \p. order to consider situations where r 

and a are comparable iii^magnitud^, or r < a, Ve us^ ^ different approximatie)n for 

R in (2). Thus^we now relax thp requirement r » a ^d assume instead; that we 

.restrict tjie direction of'obser^aticjp.Jo the neighborhoods of the forw;rf^scattered ' 
(the direction of incidence) and back^ scattered directions. Vojj>flg^llcity, we assume 

' (y - 7?)^ « , and use 




(5) ^ ' * R - 



We use (5) in the exponent of (3), but iri^he denominator (ess^tially as for Equation 
(32) of Section 6) we use only the leadutelerm R » x 'jthus 

.The* present integral describes Fresnel diffraction by a strip. 

After we have analyzed t,he behavior pf (6), we treat the analogous plroblem of 
' scattering by a circular cylinder^ A limiting case of the result We obtain will corre- 
, • spond to our geom^trigal optics results pf Sections 2 and 4 in the form 

■ • ; ] \ \ ^ ; ' 

where Ljj is the ray path obtained previously by using Herds, principle ifi^l^ LJ^Ii-' - 
its second derivative with respect to a parameter along the circle ,^arid F is the'flux 
density obtained by the Kepler-Lambert principle [KL] . , -^iZ^^:::;:^r:zr:ir:zzrz 



Our deprivation of (7) from (3) will obviate the previous special assumptions. nL- 
particular, we will not have to assume Heroes "principre of the extrailTm pltfiT 7""" " 



•i 
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Statements such ^'nature prefers an extSemum/y'nature abhors a vacuum/' etc., " * 

* r 1 • 

may be .useful aids to memory, but "natureT neithe^ ^'prefers" nor '^abhors," and 

• 11' *' ^ 

such statenients make our m}5sticism too ejipltpit. We will show that -Hero's principle 

[H'] is merely a corisequence of an approxiinate evlluattion of integral. 

The approximation procedure is known jiis the mtethod of stationary phase . Jt was / 
introduced^ by Kelvin as^ a mathematical niet!l :od for approii^atihg\^Qia^s <S mte^als* 
that come up fij^quently in wave problems^ '^he integrals we are ^ncerned with are 
of the form U =^ Juds, essentially as in (5^.^ * 

Intuitively, ihe method is based on Young's concept of interference, and the 
essentials \veYe discussed foj Figure 6-6, A cottiplex mmberjAe*^ may be repre- 
sented as a vector of lengtt3^*Avand1^irfe^on angle (phase) on an Argand diagram. 



and the resultant of a set of numbers A^e " is simply the vector sum Te*^ = ZA^^e 
shown in Figure? 8-2 or 8-3. If the phase anglbs are alt the same, then the elementary 



J 




i: ; 



/T^fGlJRE 8-3 



er|c --^ 



I'll 



vectors all point along.a straight line and T " SAn : we say that the vectars reinforce^ 
each other, or that t^ie elementary^ wave^''lnterfer.e construc^vely." On the'^er" 
hand,^if the angles are sfiiofl that. the nose of the last vector^ends up at the tail^Khe 



fin^fo^Vca^a^closed^polygon, then T 0 ; we say that the vectors cancel, or that, > 
th6 elementS^fWayes^"interiere destructively,"/ In thjg sY&alibn shown in Figure 8-2, 
the pliase (f^ .changes sli^tly (varies slowly) with increasing Ifi and the resultant 
magnitude T is qlose to its maxiraum vdlue. In th,e situation shown in ^Figure 8-3, 
the variation of (p^^ with n is large and T is small. The and may depend ori 
a parameter 0, and the magnitude T(0) of the resultant, 

(8) T(0)e'y(^) = > A„(e)e*''"^^^" 

n=l 

* • 

may assume any yalue between 0 and SAn with variation. of 

Sim^ilarly we may use the same idea for ^ integral of the form 

(9) : \ . ^ /^^A(r];0)e*^(^'^)d7] . - . - - 

The integral of (4) is of the above form, and the series of maxima and minima shown 
by the riesctttant can be interpreted graphically by means of a veptbr diagram such as 
Figure 8-2. Similarly.for the integral in (6). . 

In we can express (p in terms of a parameter 0, and then find that 4? = 0 
(and ^ 0) for some set of values of 0 (say 0^) , then we say the phase is station- 
ary at 6^ ; for these values' the situations are ;analOgous to Figure 8-2 '(nq_t .Figure 8-3), 
and we expect the resultants T(^s) to be large. Befpre discussing the general method 
mathematically, we first consider the strip jnroblem without using the method ex^c- 
itly. The results We obtain initially by relatively familiar procedures will provide a 
basis for introducing lab61i^ and concepts for our subsequent more^ggneral discussion- 

Fresnel Diffraction : We may Express (6) in terms of the tabulate d^Fresnel^ntegral 

whose path in the complex plane (i.e., the trace of the pointy ^7^:^). at§ f function t)^*^ 
i7i) generates Cornu^s spiral of Figure 8-4. The magnitude I5(77j)l - T^tji) is an' 
.oscillatory function ot 'jy\ . Figure 8-4 for (10) is a continuous analog.of such cases 
of discrete v^ectors shown in Figures 8-2 and 8-3. Although the "elementary vectors*^ 
of Figure S-4 form a curve of infinite length, the curve spirals ih around the point 
+ 4i the resultant T = t5 1 approaches the limit -p= as tj • 





. ' . ^ FIGUflE 8-4 . 

We no\V consider i7(T)) analytically. For large values of tji , we use 



CO CO 

{ -i 



Integrating the Second term by parts, we develop (^11)' as the series 



lirj)i/2 



(lay 



Thus the error 'in using the leading term ^(<«) ia proportional ro — . The leading 
term itself, i,e J, " - ^ " -1 / * ' ^ ^ 



(13) 



caabe det ermi ned as follows. 
Let us consider fii;st 



(14) • •' ^ -''^ " -^i ^"''^^ 

The squ^are of 'J may be written as' • 



(15) 



, 0 V 



Introducing polar coordinates' r^ = + , tan </? = ^ , we rewrite 
(16) . ^ ' K ' K e''^rdrd.?> =,|*.J^ e-'^rdf 

Since die"*^) = -2re~^dr , we obtain 



(17) ' ' . ' ■ 
Thus |;i4) equals ' 

(18) • 

Similarly, we have 

(19) 



J e'^^dx 



which is known as Laplace's integral . We treat the integral (iS) heuristically by 



ITT 



using B ^ - to obtain \ - , ^ 



ir/4 



which corresponds to the vectolr from th^ origin to tl>e terminal point of the up^per 
spiral in Figure 8-4. (The result may be verified by Cauchy s theoreni. If we use 
7) = ^^(1 + i)/ in (13)j we obtain J e7^ dS(i where the path is ,along tKfe; ' 

45'' line in the complex f pl^uie. We replace this path by the positive real axis plus 
the arc from 0 to ^S*" at infinity] the integral over the arc vanishes^ and the integral 
over tM real axis leads essentially to (18).) , . - 

On the other hand for small values of 77, we expand the exponential of (10) as a 
^ series and integrate term by term ^ ^ 

(21) ^(r„)=/''*[l^^^..".] 



Ih terms of (10) we rewrite (6) as^** 



•'(22) ' 
(23) 



/ 



; U - /F cl , 



Before applying thj^^ result to the strip problem,, let us first i^pply it to an "infinitely 
^wide slit,'* and determine c fpr the Huygens* tree-spacfe wavelets that, simply serve 
to regenerate the incident wave. For the limiting case ka «^ , we have t)^ i'* . 
Using the limiting values of 7 obtained from (10) and (20) we obtain 

(24) • • " I - 5(«>) - 5(-«) = 2 J(«>) = >/T e*^/^, 

and consequently > - ' ^ 

(25) *'"""* . ^ ~ c -^e?/^e»»« ^ 

But an "infinitely wide slit" njeaSs no obstruction, so that U of (25). must equal tlie # 
incident wave e*^ . Consequently, the unknown constant c of (2) and (22) for the 
elementary Huygens* sources is • ' 



4 



(26) 



/2? 

More generally, if we areldealing with the secondary sources on a scatterer excited 

by the incident wave, we. may write . - *> t * 

* ♦ • * * 

N ' f ' 

where g may d^end on the material of the scatterer, and on directions. Thus we 
may rewrite (22) as * ^ ^ { . * 

We now apply (28) to scattering by a strip as in Figure 8*1, There are essentially 
three different ranges of y that we consider. • ^ 

^lu Eigure S-.SDve specify tl^ee different ranges of y at a fixed value of x , ^ay^ 
three different portions of a screen placed parallel to the strip. We will.use (28),,(12), 
and (21) to obtain explicit approximations of U for the .three ranges, of y correspond* 
ing to th&- braces shown in the figure. The range y^ is centered on the geometrical 
proj.ection of thje edge of the strip (equivalently the neighborhood of the shadow boundary 



y = ^), and the range Vi includes much of the geometrical projection.i(the shadow) of 
the strip on the screen. 

The range of yj corresponds to t]^ » 1 and -t]^ >> 1, and include? y y 0 as the 
Special case t)^ 1 ; the range of y2 corresponds to t]^ « 0 and -t)^ » 1 ; 

the range^of y3 corresponds to -t]^ » 1 and -7)_^ » 1, If we replace x,y*by^ 
Ixl.,. , then the results will apply not only for the three sets of points (x,yi) in 
the first qttandrant shown in Figure 8-5, but also to the sets obtained in the other 
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quadrants by.^eflqjjUng the three given sets in the x-axis and in the y-axis. If, Say, 
*^^^'n\^^^> we approximate the integral hy the first two terms of (12) obtained by 
usitig (20), i.e., , ' 'v.v' • V 



(29). 



^(T?i) 



i7rni/2 



/2" I 



, • T?i » 1; 



. however, if 771 is very^mall then^we use only the leading term of (21), 

'(30) ; . i(r?i) «'tj,', tj/« 0. ' . 

,^ . At a point ,x,yi (or more generally a^thelfour points Ixl ,']yjy,^we have 
= ^\f^\(^ ± y) » 1 ; substitutii^ (29) into (28), we obtain , 



(31) 



U-g^p-*'^/4lJ(T/^) +5(.r,J] . _ • ' 

where ^g is nqtnecessarily the same for the forward scattered direction Ixl ^ x as 
for the back scattered direction Ixl = -x . Essentially as for (5), we require . 
I (a*± y)/xl << 1 ; subject to this v^e see that 



(32). 



U ge^^J^' as lx/k(a^±y)2J ^ 0 



3fe 



116 



122 



Thus in thisL]irnit;fli^strip is a ofie -dimensional seconaSry source. Taking into 
account. that U i§ in general a function of direction, we write 



(33) 



ft>r 



X ^ ^ 



In the forward direction, we require g^ = -1 imorder for a-geometrical shadow to 
exist in the sense of Section -2, i.e., ^ -e**»« '^Similarly for the cas^e of a perfectly, 
reflectii^ material, we require Lg^ I = 1 in order that the ratio of tl^e reflected to* 
incident flu^c density (^uyUi|2) equal the previous result unity. Thus we have 



(34) 



U r g.e 



•itot 



for 



X < 0 . ^ 



where 5 is a real number determined by the material of the strip (or, equivalently, 
by thQ boundary conditions ) . For pi^ent purposes we take 6 = 0 , so that 



-Ixl.y. 



(35) 



U ~ e'"°5 for X < 0 ; 



physically this corresponds for example to a water wave^or a sound wave m'cident on 

a rigid immovable strip, and also (subject to additional conditions) for the reflection ^ 

of an electromagnetic wave from a metal sjrip. 



We introduce the abbj;eviation 



(36) . . MY)- V27rk 
so that we may rewrite (31) as . ^ s 

(37) - . \ .U^ w ^ [e^^ + Ue(a - y) + u^ (a ^ y)] . 

The present results apply for the geometry of Figure '8-6.; From (5) we see that 



U^(a+y) 



Ue(^y) ^^^^^^ i- 

. . HGURE 8-6 




* R w X t (a ^ y)V2x so that the exponehts of Ue(a t y) are ^prpximations of kR^ , 
' ' ■ ' The factors kR are the phase changes of waves traveling from the edges 
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of the strip 0 , ±a to the observation point x, y So that we may interpret Ue (a ± y) as 
the edge waves corresponding to the edge rays discussed in Section 2. Thus the nor- 
mals of thQ different waves 0^(37) shown as directions of propagation in Figure 8-6 ^ 
are also the rays corresponding to such po^ints as Pg in Figuce 2-17 (c). 

In the region of the geometrical shadow the, total wave Uy is the sum of 
Ui = e^^and U ' 



+ 



(3S) Ut+ = Ui + = -Ue(^ - y) - Ue(a + y) , 

* ♦ ' ' • *■ 

ij.e., the shadow forming part of cancels the incident wave and we are left only with 

the edge waves or diffracted waves. Thus, corresponding to (38), since IUt+I^ cc 1 ^ A_ 

, * • ' k # 27r ' 

a "perfect shadow" does not exifet for nop-vknishing \ , the diffracted field is small for 

relatively small x (in the immediate vicinity of the obstacle), but itancreases in magni- 
tude as X increases and the shadow "disappears" with ir^rfeasing distance from the 
scatterer. (However, alternative forms are required for either x 0 or x «>.) The* 
field U-j-^ is oscillatory both in x and y. The flux density along the axis (y ,^ 0), cor- 
, responding to lUf+Pwith . ^ ' \ • ^ 



(39) 



Ut4- ^ , R = >/x2 + a2 ^ X' + aV2x., 



is a relative maximum, this is the analog of the Arago "bright spot" discussed for the 
•disc. In the back-scattered region enclosed by projections of the strip edges parallel 
to -x ; we have ^ 

(40) Ut- = Ui + U. = e*^ + e"*^ Ue (a - y) + Ue (a + y) 

where e""*^ is the geometricaliy*re^cted wave. 

The above Results (31) to (40) are subject to "two restrictions: the first, 

J (y - a) ^xl"^< 1 , enables us to use the approximations in (5), and restricts us to ob- 

jservation near the back and forward scattered direction^; the .second k(a ± y)Vlxl » 1 

is required in order to u"^*approximation (29) for ^. The first bounds 1x1 fro.m below 

(Ixl » a), and the se'cond^from above. Together, the two restrictions also state .that >^ 

ka -.27ra/\ »> 1 (i.e., the strip is very wide compared to the, incident wavelength, and 

^ that y cannot be near ia (i.e., y cannot approach the shadow* lines or their analogs in 

the "lit region"). The corresponding edge wave Ug'is the "near-cauj^tic" foran. ' 

♦ 

We cannot use (31) to (40) for a point y ^ a in the range of y2 Qf Figure 8-5. J[n 

that region although -7? - \/-?7 (a + y) « AK-^2a » 1 ,'we see that ' 
" . V ^ixl V ttIxI 

" V 7^ " y) ^ 'l\i\xs, in the general form (28), we still use (2?) for J(-r?j\ 

but we must use (30) for 9(r\^. Consequently at tlie four points Ixl , lyft-, we obtain 



t - 



124- 
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(41) U = qFe^^''^' 



V27rJxt^ . <^ 2 V27rk^ 2a' 



The second and thicd terms'aile of the form considered in (31), The first is a cylin- 
drical wavje, i.e.^ tKfe^A^a^^f a line source on the edge — a "true*' edge wave dfecreas- 
ing^as l/\ir with increasing distance. The third term becomes negligible for very * 
large ka , for whldruasnBtlife total fteld for x > 0 reduces to 

(42) Ut+ = U+ « f el^; H(kx)^^^i%:^ , H(kx) ^ \e'^J^ e'^"^^ , ' 



Xhus the field m th_g fteighbo.rhoocl of jHe shadow -line is half the incident wave plus a 
cylindrical wave corresponding to a line source with source strength proportional to 
k(a -y)\ As y a , we «ee that Uf hxTl linearly; this holds, whether y Approaches 
a fi^om above from*" above" or "below*' in Figure 8-5, Similarly for x ^ 0,ye have 

J (43) Ut« « + ^ ^ii^^^^-l-Q-iy^ •f'H(klxl) ^i^^^ . 

In the range of y3 in Figure 8-5, we have -t?^ » 1 , and again i^e (29) foi? 



both Fresnel integrals in (28). Howevel^^ in contrast to (31),, the scattered wave at 
the four points Ixl , lysl , is given by 



(44) . ^,:,em^i^—^im^^_^ — r^i^, Fft— J 

,, . ? ^-Cue(a -y) + Ue(a..+\;ifl , - 

SO that such points^ receive only the edge contributions of (37) , 

the above explicit approximations suffice for present purposes. A more com- 
plete discussion of the problem of the strip is given in introductory texts on optics in 
which the field at any point in space is usually computed graphically from" Cornu's 
spiral Figure 8-4: • *^ * , ' 

Thus up to moderately large Ixl the scattered wave is largely "confined" to the 
*^tr4p iy{n< a;.'ifl,lhe senSe.ihSt*ivithin this strip we obtaiathe geometrically, reflected 
and shadow forming waves. These two waves correspond to the waves scattered by* 
an infinite plane, however, superimposed on these are the additional waves, that we 
, interpreted as edge waves, jgecausfe of the additional waveis, the^ shadow is not * 
**perfect;" careful observations (subject t^he present restrictions on distance ^ 
parameters) on th.^ shadows of scatterers having very regular edges show a system 
of bright and dark bands parallel to the edges of the scatterer (a "fringe system"). 
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We now determine the.numbex' and separation of such extrema that may be observedL. 
in .the shadow region on a screen_^arallel to the strip, , ' 

-We use • ' . ' 

(45) ^ \. ^ F(7))dTj =-T(£y2£ , &=£(y), 

, . / 

to differentiate'^U. as' given by (22) and (23), and obtain ' y 

The^ extrema, oC-^O obtained from dU/dy = 0 , correspond to 
(47) " ^ e*'^'»5/2 = 

Consequently !the exponents must satisfy \ 

(48i ' . f (7)2. - 7)b ay = 2n^; f n = 0 , ±1 . •••-^ . . ' . 

The separation of extrema i-Spttius^'^'^ \ . |\ * £ 

(49) -\ |y„., Ay = 

and there are N extrema, with N given by , 

(5b) >^ ■ N = ^ = 2k^ = • 

^ . , Ay TTx Ax ' 

in the geometrical projection of the slit. Thus the number of extrema increases with 
increasing strip width, pr<.with decreasing .wavelength, or decreasing axial distance,. 

Before continuing the main line of this section, we consider the range of very ^ 
large I xl excluded in the discussion based qti (29), If Ixl becomes very large, so ' 
that 7^ i^i ' ^^ ^0, then both n+ ^ ^ and 7]_ ^ 0 in (28), and we approximate both 
Fresnelintegrals by means of i(ai)XThus * - ^ 



Thus for this casp,,the scattered; field is essentially that of a line source of strength 
ka We derived this resultivia (5) which means that it is restricted to angles near, ► 
.the forwaird, and back directiojisr^cotiaparing with (4)> we see that (51) is merely the, ; 
specisil case of corirespondingjto 6 r 0 , tt , and that for other values of 0, the 
appropfiate form of U^t^large distances is simply dbH(kr)kar(0), ^ 



Method of" Stationary Phase > We have discussed the prelimiiia'ries, and can'now^ turn 
lo the main, topic of this section. In general, we consider an^ntegral of the form 



a 



Where G is a slowly varying function of x compare^ to e^^^W in the sense that -CL, 



changes billy by a small fraction of itself when kL changes 27r . If there exist'Sne ^ 
oi; rpore values of x for which ^ vanishes, then the principal contributions to the 
value .of the mtegral arise from the neighborhoods of the extremf {or stationary 
vilues) pi L , elsewhere the contributions cancel through destructive interference as 
defined «)reviously. We reiterate that the intuitive basis of this idea is the recogm- 
tibn that on an Argand diagram (as in Figures 8-2, 8-3, and 8-4), ,1 is a sum of ele- 
mentary^ vectors whose direction (essentially the phase kL) is in general a rapidly 
cl anging function of x, so that the resultant I is consequently small. However, if 
there exists a value of x for which' ^ vanishes, then the phase is stationary at this 
.value and only slowly varying in its viciiuty; the elementary vectors near this value 
are almost in phase and add to give a large result. 

Proceeding analytically, the Xaylor's expansion of L(xj around some v^ue Xa is 
,($3) ' . L(x) = L(Xa) +• L»(Xa)(x-x^) +X"(Xa)(x -x»)V2 + , 
v here L*(Xa) meanj^differentiate -J.(x) with respect to x and then set x equal <t^ Xa*, 
Jiixiilarly. for the second „derivative L** , etc. If a value x, exists for which V (Xs) = 0 
men ; ' • " *^ . * ^ / > 



r . ' ' , (X - x^)2 ^ 
L(x) = L(x,) +;L"(xw) o + 



r Assuming that JUHx^) ^ 0,vand that the higher ^rder terms are negligible, we keep 
^ only up to quadratic terms in the exponent of (sa). We replace the slowly varying ^ 
, function G(x) by its value G(Xs') att^the statl9iiary poinjfcjthe point marking the center 
of the region in, which the integrand contributes significantly)T^and wotk with 



I «^ I, = G(Xs)e 



ikL"(x3)(x-xJ2^2 



dx 



(55) 



where^fc J is the Fresnel integral ate in (lO)ff. 



.4 



[C 



In' particular,?! ir)^ - «>, then 9 - j7(«) of (20); we have 5(«>) - 9 (''^) ^ yf2 e'^^/^ 
and consequently ' 

s 



(56) 



Gse 



Where the subscript s indicates that the function is evaluated at the stationary point 
* Xs. (If there is more th^n one stationary point, then Ij is a sum of such forms.) - 
' -Compare (56) with (7). v ^ . ^ 

Before applying (56), lit us show that it is actually much larger than (52)^ for the 
^ case where the integral haJ no stationary point. .We show that if there are no statioif- 
ary^ values of L , then the integral (52) is only of the order ^ as compared to I3 
which(from (56)) is proportional to";,^ ; since k > > l , is therefore much larger. 
< To see this we introduce. L as th^'new inte^ation variable: 



and integrate by parts in order to develop I in«{5owers of ^ : 



(58) : 



^ . JL(-a ) • k^dLl L- ) 



.ikL +' + 
- L(-a ) 




Thus as^ long as L'/G does not vanish in the range -a'-^ to a^, the. integral is 6hly of ., 
order j;^' and is therefore much smaller than the stationarj; case Ig of (56) * ' - 



-^ As a fii:syillustration, we^pply (55) to the original integral (3) for the strip with 
the constant ^iven by (27): - - 



(59) U = g-^/^e-"^/" I , 

Comparing witH (52), we'W ihat 0(77) =^ 
<?'..asln Flgiire 8-7, we differentiate to obtain 1^ 



1 



, and that L(r?i = R())) . Inti:oducing . 




, FIGURE 8-7 



i 



(61) R" =. , ' , -= 



•The stationary values correspond to R' = 0 : 
(62) 

Consecjuently 
« 

(63) 



R2 j" « • 



,2 



R. R, 



Rj = ±x = Ix 

Substituting iijto the integral of (59) via (55), r/e obtain 
(64) 



1 

Ixl 



,g ii>xo me integral ui \o^} via vi>ii;, uutcuu 



•y) 



Introducing (64) into (59) yields the earlier form (28), The limiting fOrm based on 
(56) gives U = ge*^'''' as previously, with g = -1 from our "shadow condition," 

We may generalize tl^e above direqtly to an arbitr^y angle of incidence as in 
Figure 8-8« The'ihcident plane wave may be written 



a 


-^-^ 




_j — ■ — — 


0 


« 


! -a 





FIQURE 8-8 

(65) Ui = e"«=~°^"^*'"". , . . 

which is merely the form of (1) obtained by routing the xy .coordinate frame through^ 
an angle -a , Since the phase of (65) is zerq at the origin (tifie center of the strip), 
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the wavelet originating at the origin has the same phase as obtained from \% iSdr-''^' 
ever, the wavelet originating at tj is excited by U, (0, t;) = e""> T,so"that its phase'^ ' 
contaiHS the additional term ,k7f sin a. Thus instead of "(2) we have ' J* ^ 

- U(7?) =»C = , 

and we replace (59) by ♦ ' 

where the appropriate value of c will be aeterrhined from a limitirig case, 

CorrespcHidingto (52), we take G = c/./fcR to be slowly varymg, and differenti- 
ate the phase ' • 

* - • ' * '* 

(6^) L = 1] sin a + v/x^ + (t? - 'y)^ , * ' 

with respect to tj We now have • 



(69) ■ ( L' = ,sina + simp , L" = .co£iP' 

R 

♦ • * 

The stationary values correspond to 

(70) L» = sina + sin<^ = 0 ; sin<^ = -sin a ; » </> = -'a, tt + a , 

which contains Euclid's principle of^ef lection and the principle of stiadow formation. 
Since sin (/> = J)^ = ^cj4<P , we see from (70) that y - r?3- = Ixl tan 
Rs = xsec(/? = Ixl sec cv, a^ . ^ 

(71>L3 - nsSina +xsec|^3 = (y - Ixl tana) sin a + 1x1 sec a =^ ysina + I xl cos a 
Similarly^, ^ ^ ^ 

(72) ( LI! ^ cbs^ g 

Substituting Wo (67) in terms of the limiting form (56) we obtain 

^TT +fkL +iff/4 27r 

U e 



= g±ll«cos a+llcy slna \/2t[ ^^^/^ ~ TT 

kcosa ^ _ ± • 



Comparing with Ui of (65), we'determine Cg from the shadow condition = -Ui: 
(74) / ^ ' • C3- = l*£Osa , 
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which differs from our earlier result by the additionarfactor cos a . Sirfiilarly thecor- 
responding value of g for a perfect reflector as for (31) ff, is now replaced by gcos a 

Using the more general form corresponding to (55), we now tiave 



xl sec a 



>/2 

■cosa(±a +.1x1 tana- y) ,y^^ 



Circular Cylinder > Let U3 now apply the same method'; to consider scattering of the ^ 
plane wave (1) by the convex cylinder as in Figure 8-^. The point a(<^) on the 




FIGURE 8-9 



cylinder has the coordinates a cos (p , a sin t^* , and the excitation at the point is • 
gika cosc>^ write the corresponding scattered field _as the integral of u over lh6 
arc ad</? : ^» - , , . 



(75) 



U = / ,^plk(R^aco3<y>)ad(P . 



\e do ilol know .c - completely, but the result ;(74)''^iid .the corresponding form gcos a 
suggest the generalization ^^"whjch-.f^^ (thS^,^^^of with respect to the sur 

fac^ normal) is replaced by the^prejs^ent analbg^^r^- V ^ Thits for convenience (and as 



can heJustified with a mor^ coniRl^te model):we use 



cos<^d<^. 



ancTwe shainrfee that cos'^, Vthe integrand is appropriate fo^^botji^geometrical re- 
flection airf^hadow fonpation. We consider only the range 0 < '6 < t explicitly; 
however, the resutts maybe extended to all'e by introducing absolute values of the ; 
* trigonometric functions (as required ta preserve symmetry). 



The phase of I proportional to ' ' ' 

(77) - , L = R V acos.^ .= s/x^ + - 22ivTos(<p - 0) + acos"^ , 

' and its derivative ' / ' 

*- f ♦ 



vanishes for the two values <P ^ <Pi^ or <f>^ , such that 
: (79Lr 



\ ^ 
sin^^L - ^) sin^L 



(79D). 



sin{(p^ -6^) ' ^ sin(7r - ^j^) 



: P9r a given value of 6, the phase has only one stationary value: the value may cor- 
respond either to geometrical reflection. as in Fi'gure 8-10, or to forward scattering 
as in.Figure 8-il. The first value applies for y in the "lif^egion L in Figu're 8-10; 
the seco^id value , which yields the shadow forming j;^j^rc5h:esponds to y in the ' 

^ "dark" region D as in Figure 8-11.. Equation (79L) corresponds to Euclid's princ'iiDle 

^oCrreflection, and (79D) to the prinoiple" of shadow formation. 




. ^ , FIGURE 8-1^ 



.the second ^ierivative of L equals 



- X „ - a^r^sirt^fcp - fl) , / ■ 

^ + p[-cos(^P - e) - apos^ , 
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and substituting (79) leads to the Special values of the stationary points. Thus for 
either'case 



Ly =-J- [-asin2 9>, 



(81) Ly =':^ [-asinVs + rcos(^s- 0) -Rs'cos^s], 

where Rg and t/?s are the special values shown in efther Figure ^-10 or 8-11. At a 
forward point, we 'see from Figure 8%11 that rcos(7r - + 0) + a = R^^ cos(7r - (p^); 
consequently, rcos(<;^p - 0) - Rqcos^^ = a, ^ (69) reduces to 



(82p) 



On the other hand, at the reflection point, we see from Figure 8-10 that 
rcos(^L - ^) = a + Rl cos^tt - <P^) = a - R^cos <p^; consequently 



. R 



-' [a(l r sin^ <p^) - 2Rj^ cos ^^J 



(82L) , 



' V . 2acos(7r <p^) p 



R, 



Similarly the Stationary value of the phase function at a forward "poiiit- equals 
= Rp + Aicos acos(7r - ^y^) = x . ' 

Although we could eliminate and (p^ in the reflected value, it is simpler to leaVi^ 



(85D) 



^ f ' 



4»|-.JiiJhe^orig^ form v 

(83L) , — ~ " acos(7r -^L) * 

' . ' ' ' , - ^* 

.here (for a given value of r and P) Rj and (P^ are determined as in Figure 8-10, 
' • " ' ' * • \ 

i.e.,, by noting which ray (or the corresponding value of yj of the incident wave front 
*^ ^ ' - • * *>*• • / ' 

^. can reach v{Q) via reflection in the cylindrical surface. " ^ -> 



^ 1 



1 -*-v 



lis 3 



Finally, the required slowly vary^g Rar^^thfe^^^^gtlg^^ at the , 



Stationary points are given-by 
(84) 




(85) 



Substituting (56)^into (76), aoH then entering (84), we obtain ^ 
Us = -gaG, #,irgaco^'^s; - 



Thus for the perfect reflectorN^J^ 1 ,Ve enter the above D-values in (85) and obtain 
the shadow forming wave • ' , ^ _ < 

(86) ' Up = -e*"« = -Ui . ' ' • 

Similarly, the L-values give the geometrical reflected wave 



(87) 



Ul = 



a COS a 



~. lk(R, - acos^a)' 
^ ^ *- 

2(R^ + |cpsa)-y'- 



Of 



) 



which we considered in Sections 2 an^ 4 ^rom a geometrical basis. In particular, the 
flux density ratio F/go of (4:14) is simply the present IUlI^ , and, similarly thp 
previous (4:18) corresponds to lUgl^ of (85); the present forms are richer in that 
they make the ray path (Ls)^s well as the caustic (L'J) explicit, ' 

The corresponding-Fresnel approximations are obtained by using (55) instead of 
(56).. We chaise variables to t? ' - a i^n<^ , D L = aciDS <^ D 'l , and obtain 

; (88)^0 W>= Uj, ^ - 9(^j\\ r,^ = ^(±a - i) , 

wl\ere is given in (86); thus (88) is^ simply (28) for the range x > 0, Similarly, 
in terms bi^^ of (^7^" we have* 



(89) 



2k(RL + 2 cos a) 



ttRj^ a cos a 



(±a - a sin a) . 



.a = TT - 



where'i7(T7) is the Fresnel-integral^e in (10) ff. 



The field on caustics : Although supplemented by phase considerations, equations (86) ^ 
and (87) are still res^ults of geometrical optics which we could ctJnstruct piece-by- 
piece from the spegial; "laws'* of the earlier sections: [HM «^ves the directions, IKLJ 
the. magnitudes, and the phases may be' obtained from Newton's idea pf periodicity, 

^ > V * • • A 

However we have now obtained these results ^seatially from the single idea\ of 
> ^ * • 4 ; 
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periodic waves that evolved through the work o| Huygens, Young and Fresnel into 
' (YF]: U = /uds , ^^^^^^ ^ . ^ 

"which represents the scattered wave'as an integral of elementary cylindrical waves ^ 
over a line or an arp (as in (3), (e7), and^TS)), or of spherical* waves over a surface 
or volume (as in (6:38)). Stai;ting with [YF], we used mathen(iatical procedures to 
approximate the integral, in particular, we obtained the general short -wavelength 
approximation (56), two of whose i^ecial cases are given by (86) and (87). Aside from 
. . superseding the earlier .special laws^ [YF] i6 applicable to many other phenomena than 

' covered by [H'j and [KL], or by (56). We have already applied [YF] to obtain the 
Fraunhofer, 'Fresnel, and Rayleigh-Born scattering approxiiiiaj;ioiis , more of which 
are covered by '(56). We now apply [YFJ to supplement the short -wavelength form (56^ 
by determining the magnitude of the field on a caustic, the case L'J ^ Uj'j ^ 0 excluded 
in (54) and in our discussion* of [KL]. Thus, for example, the anal^)g of (87) for reflec- 
tion from the concave semicircle has a replaced by -a (see 4:15)), and does not hold 
on the caustic ^* 0, Rl = -(a/2)cosa. However, we do ndt require a special "law" 
to obtain a non-singular form, but merely a more appropriate approximation of [i^F] 
than given by (56). ^ , ' - , 

On a caustic, both L* and L" vanish; in addition, on a cusp of a caustic (where 
the derivative of the equation of the caustic vanishes, since the curve changes direc- 
tion) the third derivative must also vanish. Thus for such cases we c^n^no 
longer approximate L(x) by jneans of (54), i.e., we riiust keej), Additional terms in the 
Taylor series in order* to obtain th^ first correction to L(Xs) . 

Let us assume that the first n.- 1 derivatives of L at vanish, and approxi-- 
mate L(x) by the ?aylor polynomial of n-th order: ? 

* * n / \ L ^ 

' (90y ■ * - L(x) « fi(x,r + ni^'^" lT^ * ' ^' " "^"^^ ^ ^ ^ 

For this case, for infinite limits ^corresponding essentially to (56)^ , we have 

If we introduce a new variable y through y" - k(x -'^iJ^L^^Vnl , then w^ may re- 
write (91) as . 

where {he remaining integral (the ^gamma functioni* integral) depends only on n . 
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' Using (92) instead of (56) in (76), we obtain t . ' 

_\'rm)} : : ' u ~ c„e"'^»(k)^"n , 

* ^'^ - — > ; . ^ ' ■ 

where we have suppressed practicfdly everything but the dejSendence on k = 2n/X . Thus 
' away from a caustic, we have n = 2 and luf is independent of wavelength (e.g., as in 
i^r^^ f"^ correspond to "true^eometric,aI optics." For the line caustic of • 

circular cylinder, we h^n = 3, an|fsince ' R(?0 « a(?) on the cau&tic R = |qoso, '"^ 

^7^-^fi^-ftSPj93)-«ftd from dimen^iojial considerations tha^ U (ka)6; simiIarly,for a 
'^."^^^^^^'^ circle), we have n = 4, aid U oc (ka)"?. More generaUy>since 
j - • -^Cx^-aTwrsse-that- k2 " n increases with increasing n and approaches ki (at which 
limit the- phase is "completely" stationary). Thus U increases on caustics as the 
\yavelengtR A decreases; the field is always finite since the case of zero wavelength 
(thg implicit assufnption pf craventional geometrical optics) is not physically reaUzable. 



9, Mathematical Model for Scattering s 

/T5y>revious sections we considered certain aspects of wave theory but based the 
development'on several supplementary "la^vs of nature/^ In the pi'esent section we tie 
..these special postulates together into a mathematical model for scattering. 

W^e Equation : We considered the plane*waves ^os(±kx - u;t) , and, for convenience, 
worked with the real part of the Qorresponding expolientials: 

(1) ' e^*^'*^^ ^e***". e"*^^ = f(x)g(t) , \ . 

Equation (1) is a product of a function of x times a^'function t, each of Which satisfies 
- a second order differential equation: , 

(2) . • k = f ; ■ ' ■ [_ ' 

(3) = -a;2g(t) = .k2v2g(t) , . a; ='kv ^ ^ ^ 2irp ^ 

• As discussed previously, \"Ys the phase velocity of the wave generated by a source 
vibrating at a frequency p, and X is the wavel^n^.— the distance between crests^ 

fhe general form of (2) and (3) ia** ' m # . ' * / 

(4) ^-TT^^^ ^2F(y) = 0 , 

whose, general' solution equals ^ 

(5) Cj cospp^ C2 sin/3y ; or,'equivalently, e*^ + D2 e"^^ , 

where cthe first con^ants (Cj, C2) are linear^combinations of the second (Dj, D2). Thus 

%in choosing the particular combinations that led to (1), we used some selection rujes. 

•J 

* We discuss these rules subsequently. - 

Now let us use the*||)bove to construct more general equations.^ Our attitudeaa.^^ 
the following. We know of phenomena that can be described bj^3^lif^^tihctienB^uch . 
as (1). -Let US seek a general wave equation that yields (1) as well as more ge^neral 
wave forms. The more general waves may well correspond to phenomena not covered 
by(l).' ; »^ 

From (2) and -{S), we have ^. , 

^ j^'d^m = _i_'di£(ti ^ 2 ' 

" - . dx2 ^ \ V2g(t) dt2 • 



fiyf". Subtracting one from the other, we obtain 




i_ = 0 



(X) ('2g(t) dt2 , 



I 



\ - \ 



or equivalently 
(7) - 



diiM-.lMdi^.o. 



where ^^perates only^on f(x) , and ~ on g(t) . The present notation is awkward. 
We would like to combine f (x)g(t) in a single form E(x, t) . To do so and4>j<'eserve the 
idea that 'the differentiations with respect to x and t are independent, we introduc^ 
the notation — to represent differentiation with respect to x while t is fixed — 
partial differentiation ; similarly for t. Thus, we rewrite (7) as 



(8) 



^ /iL 
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v^at^ 



E(x,t) = 0 



This is called the wave equation. The wave functions of (1) are special V)lutions of 
(8) corresponding to periodic waves. 

We generalize (8) to two spatial dimensions x,y by introducing an additional 
operation — r into (8): 



(9) 



-2 ay2 



32 
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E(^,y,t) = 0 . 



The plane waves e**^ cos a* iky sin a - ia.'t^^at we considered in Section 8 ere solutions of 

(9). Similarly for three spatial dimensions we introduce an additional operation — 

X . 8z2 
in (9), We write the general foian as ^ 



(10) 

Where 
(11) 




E(r,t> -V2E(r,t) 



v2 at2 

1^ 



v2 = iL^ 

8x2 



ay2 az2 ' 



(which is also frequently written as A) is called Laplace*s«pperator. 

We ma^ also recast (11) in polar coordinates r, 0, <p. In particular, Jhe elemen-^ 
^tary spherical wave that we considered previously. is the special solution of ^10) that 
depends only on the magnitude r but is independent of 0 and ^ . Thg^simpier equation 
for the elementary spherical wave ^ 



(12) 



E(r,t) 



c^be obtained by comparison with (1) and (8). Thus if wes.replace E(x,t) by rE(r,t) 

a2 a \ ' * * i 

and — by — , We obtain the corres{fcnding equatioti for (12): * 



ar2 



(13) 
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We may rewrite this directly in the form of (10): ' , . ' 

The general solution of (8) may be writteil • . . 
4^15) • . E(x,t) = F(vt.-<x) + G(vt + X) "^^^ ^ :r '"^^^f 

where F and G are arbitrary twice.^amdyentiable fimotiong.. Similarly, the general 
solution of (13) is ^ • . ' ' ^ 

•(16) rE(r,t) = F(vt - r) + G(vt^t r). . ' ' • 

The corresponding solution for the equation^f the line source in two dimensions and of 
the general equation (10) caraiot be expressed so simply. W« mention the general 
solutions onl> tp stress that the solutions correspondii^ to periodic waves are special 
case6. 

Let us now ignore practically e^verything that led us to the wave equation (10). We 

^ accepjt (10) as fundamental and seek its periodic solutions, For completeness, we 

repeat the definition^ of the fundamental paraineters given in previous sections^. 
\ • * 

The time-periodic waves we jsonsidered correspond to solutions having the product 

form * • ^ ' . 

(17) E(r,'t) ^f(r)g(t) . 

If we substitute (17) into (10), the "variables separate" in the sense that we obtain 

f(?) v2g(t)- dt2 

essentially as in (6). Sin^re the left hand side of (18) is a function only of r , and the 
rigljt hand side only of t , each side must equal the same constant; call this constant 
-k? . Thus (18) reduces to ' • • * ' 

(19) r . ^ ' . . 

(20) • ' : — v2f<?->4M^ff5'h=-o^^ • ; 

where (20) is known as HelmhoM^S^uation.'^ ^ ^*v* ' : 

Equation (19) is the form (4) we considerSd^reviouslyJ Its solutionsja^ the 

, periodic functions in.. (5) /;Wffl < " 

'A 

(21) ' 'i^^^^^F^^ ^" ^ ' 
to work with, fn jequaiiott (XOr^^sigtisnSO!^^^^^^ the distance an element of 
the wave coders in unit time. From (21I,jve see that g(t) is periodic in t , i.e., if the 



time t changes by multiples of the constant T= -l!^, then g is unaltered: 
•(22) g(t) = g(t+my) ; T = ^ = .1 / n = 1,2.3 

■This J.is the periodicity of the wave in time, and i> (the frequency) is the number. of 
^timd'tHat g has the same value in unit interval of time. The analog in space, the" 
wavelength \ ,= vT = is the distance covered by" an' element moving with velocity 
V for a time T . But from (21J . we have J = i . Consequent^ k = f .is the rela- 
tion between the "separation constant" (the propagationJactor.o.r''wave number) and"^ 
wavelength. ^ ; , ~ • 

The space equation (20) is known as the reduce^ wave equation or»Helmholtz's ; 
equation. The one dimensional case . = 0 is given in (2), and the spe- 

cml case of the spherically symmetrical \Kayfi is implicit in (14), i.e., 
■^(rf) + k^rf = 0 , f = f(r) . 

The above equations specify propagation of waves in a medium whose pro^rties 
are determined solely by v . For the periodic cases, once we fix the frequency factor 
CO. the corresponding wavelength in the medium is determined. If we are dealing wfth 
several such media specified by different velocities v„ ; m = 0, 1. • . . . then we ' 
obtain the same wave equations with v replaced by v„ ; the corresponding reduced 
wave equations for frequency factor co involve k„ = ^ = 21 ^^^^^ ^^^^ 
earlier equations in v . we take Vo = v = constant as aVference. and write 

where is the relative index of fefraction; consequently k. = M.k . The oQjrespond- 
ing space equation for a mediuip specified by /i . is * 

(24) - ■ (V2 +/i2k2)f(r) = 0 . ' 

e.g., 



(25) 



for the one-dimensional case,. If "is ihdependent of* x . then the solutions of (25^ are 
fhe forms (5) with /? = /ik . ' , - • ^ 

Conditions on the ^lutiWn: All the problems we considered are described by functior^s 
E(r)g(t) = E(?)e-^^S where E(r) is a particular solution of the reduced wave equation 
[I]:. (V2 + ii^k^)E{r) = 0 . 

The particular solution is determined by constraints that have been implicit in our 
development. The constraints are of two kinds:^ ' * 

in]: . restrictions on the solution at the Scatterer's surface \ 
^ my.. restrictions on the solution at large distances from the scatterer. 

. , ,140 ' ^' 



The adcJiticm^rconstraints are necessary h^caus^ the wave equation nferely 
describes the local properties of the inedimn and how^a w^ve travels froni' point to 
::,point, But what if th^ medium is discontinuous? e.g';, ^ ' - ' 

^ {l}: • suppose we have a glass of water and consider waves on the surface of 
. - ' the waterl bounded by the unyielding rim of the glass; 

{2}: ^ suppose we are in a^b6at on a very large lake and the boat is an obstacle 
for an* incoming wave? 

. Cases {1} and {2} illustrate two essentially different kinds of wave problems we may 
be concerned with. J ^ -^^^ * > • ' • ' ' , 

In {4'^e,'aeaT^th a bounded medi^nl.• we are givfen v^ the shape of the boundary 
and constraints on the soliiti<Jii at the boundary, and then mav,seek to determine the 
^-forms and periods of the Wves that caabe hiaintained*^ sim enclosed media. These 
are free vibration problems: the waves on a taut clotheslixj^, ,^e waves on the sur- 
face of a glass of water, ^ sound waves in a closed i^m, th?eleoti:omagnetic waves 
in a metal cavity, etc., are illustrations, and analogous probl^s exis! in the quantum 
theory of atomic gtates. ^ \ * 

If the bouiidihg surface is on\ that y^ds, then waves on the inside create waves 
on the bounding surface, arid they kay propagate in a region external to the surface. 
We may also set up vibrations on a surface and use the surface as a source of waves 
for the external medium, e.g., a vibrating- drum h^ad as a source of sound. All musi- 
cal instruments, strings, drums, pipes are examines of "vibrator-radiator" systems 
for sound. (We have switched from talking about light to talking about ^ound and water 
waves; this is*partly for conyemence^ but also to stress the fact that as far as wave 

physics goes there are analogous phenomena in all branches of science.) 

<k 

In {2} we deaHvith Abounded obj^bfthat represents an obstacle to a wave travel- 
irife iti an essentially unboSided medium. We require conditions that tell us the shape 
^ and size the obstacle, whether its surface is penetrable l>y waves, and, if so, then- 
; what is the medium inside its surface. Such boundary conditions or transition condi- 
tions specify the kind of discontinuity the obstacle represents in the imbedding 
medium. Depending on the phenomena we seek to model, we may require boundary' 

conditions such as , * • 

it \ 
ina] E = 0 on surface • ^ • ' . 

or ^ • , . " 



inb] -r-- = 0 on surf ace < *\ - " ^^ ' ^ 
on " . ^ 

9E * 1 - ^ .*••'< 

where, ^ is thei rate of change of E along the riorihal at a point on the surface. 

^ These conaitiotts correspond to surfaces impenetrable fo waves. If the surface is 



penetrable (partially transparent), then many phenomena correspond to the following, 
the waves outside the scatterer's surface travel in medium- 1 and the ^vave function^ ^ 
satisfy (V^ + ^^^^)Ei = 0; within the scatterer they travel in medi5)im-2 and satisfy 
(V2'+ ^^k^)E2 = 0 ; at the surface and are related by the transition conditions 

8E- ^ 8Eo. 

where A is a supplementary ph/sical constant. Thus in general, the wave problems 
we coasider are specified by two physical constants (or "physical parameters*') M2/M1 
and A whose values must be assigned at th'e start. 

Having [I] and [11], we complete the mathematical statement of the scattering ^ 
problem by conditions at large distances from, the scatterer [HI] . These specify that 
,\ve seek a solution consisting of essentially two terms: one teym corresponds to the 
incident field, e.g., 'a plane wave ' ' 

inia) . . > Ei = e^^ , . 

which i3 the space part of e*^'*^'*; the other term, say Eg, corresponds to the outgoing 

w^e radiated by the obstacle in response to Ej , In (6:40) for scattering by a 

tenuous sphere, we saw that. the wave at large distances from the scatterer was the 

product of a function of directions and the elementary outgoing wave of a point source, 

- — . Similarly, in (6:32) for scattering by a strip, the wave at large distances 

was proportional to that of a4fn^sojirce - — p= — . The Corresponding wave surfaces 

V r " . . 

are symmetrical, and, in addition, \^e s^Wlh Sectibn$ 2 alid 5 that the eikonals cor^ -vr.. 

responding to geometrical reflection from a semicircle and hemisphere, although com- 
plicated in shape in the vicinity of the obstacle, became more and more symmetrical 
with increasing distance. Similarly for ^planar s.catterer (6:17), the wave e*(^'*'*^*) 
is symmetrical in 1x1. Suppressing the time dependence e"*^* , we summarize all 
such c ase s by the statement 

[nib] ' ' Es„j - grn -J—— as r «; m = r, 2, 3 

where r is measured from some point in the scatterer, and where g (called the 
scattering amplitude), is independent of r . Thus at la^ge distances from the scatterer 
^r^'^ 00), Eg reduces to the* elementary sjonmetrical wave times a function of angles. 
The condition [Illb] is a consequence of the weaker Sommerfeld tadiation condition 

Limr(""" ITg^ -ikE5m| = 0 , as .r «, .ahd of ^till weaker condition that th^, ^ 

integral over a sphere of radius r of the absolute square of the function in braces 
approach zero as x The condition for r «> insures that we deal with outward 
radiation, and that the scattef er correspond. to a source of waves (instead of sink of 



waves, and also rules out free vibrations as in {l}). Collectively we write the total 
fielcras' ' / » ' 

imj . E - E, + E3 ; E, = e^^ ; - gm i ^ = 1,2;3 , 

where i and s stand for incident and scattered respectively. 

P^uations [I], [11], and [III] 'constitute the mathematical model for scattering. They 
replace all the, special princij;. we considered previously; they cover all the cases, 
where the principles apply, and many additional ones as well. TJrfey incorporate the 
essential gross physics of the effects of an obstacle oi^^he radiation from a source, 
and the cumulative fruits of two thou§and yeajs^ e^mtion in providing a framework 
of well'posed problems with unique solutions. In suniitiary,. the wave equatibn [1] de- 
scr^be^ the local properties of the me^ia, the*sUrfaceit)nditions JII] take account of 
obstacles (interfaces, transition regions), and Bie. conditions at la^ge distances [III] 
specify that the field consists of a wave (Ei ) from a primary source perturbed by a 
wave (Eg ) outgoing from tlj^e, obstacle. We can now seek analytically the redistribution 
of the radiation of a source arising from the presence of an oHstacle. 1 

Point Scatterer : As an elementary illustration let us consider thje scattering of a 
plane ^ave e^^ by a small sphere of radius a for the boundary fcondition^flla], E = 0 
at .r = a. For the gfeneral case of a sphere of arbitrary radius a we^would work with 
the complete solution of [I] for m = 1^ subject to (1113], We would represent Ei ^and E, 
in terms of angle -dependent functions and initially unknown constants, and ^hen use 
[Ila] to determine the constants. However, the restriction a « \ , or equivalently, 

Ml 

(26) ka 0 Ir • 

simplifies the problem. From the geometry of Figure 9-1, the incident field^ e^^^ 




' ^ FIGURE 9-1 , 

' ^_ -1 , • . • 

equals at the surfaqe of the sphere; using the restriction (26), we have 

gika cos v> « 1 so that jwe may work with the approximation 

(27y , .1 : Ei(a).= . ^ ; 

Thus the exciting fiel^l aVthe surface is independent of angles, and the corresponding 

scattered wave must |)e similarly; indeppndent of angles: E^{r) is a solution^of 

^(rE,) +k^rE, . . 

dr^ 



0 , and the only one satisfying imb] at large distances is 



i lis 



(28) . ^» - " r 

. ^ e^kr 1 

where C is a constant. At the surface of th^ scatterer r = a . we have ^ « ^ ' 
and consequently.the total field at r = a is approximately . 



-(29) 



E(a) = Ei(a) + E,(a) = 1 +" 7 



Applying the boundary condition E(a) = 0 , we get 



(SO) 

The scattered wave for r > a is'thus 



C = -a , 



(31) 



This corresponds ph^^sically, for example; to scattering of under^vater sound by a 
^mall air-bubblfe (an exception to (R) of Section 7) . 

Slab Scatterer: As another example, let us consid^ scatt'ering of k plane wave by a 
- partially transparent slab as in Figu're'9-2. The conditions on the problem are; 



(32) ; 

(33) ' 

'J 

(34) " 

(35) , 



+ kA El f 0 , 



Ixl > a ; 



E2 



= 0 , K = ; Ixl <.a ; 



dEi 



• ^1 " 1^2 » dn 



= A 



Ixl a 



El = El + Ea 



El = e 



Hoc 



dE2 
dn ' 

Es ge as Ixl. 

0 » 



From (35), we write 



r 4 



(36) 



, ikx 
e 




X > a 



= g.e 



-ikx • X'< -a 




I * • 



> q e 



ikx 



1 ^ 



-a 0 a 
' FIGURE 9-2 
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From (33), we take the most general solution in the form 



(37) 



E2 = b eiKx+sij_e-iKx , 



^ 



We thus have four constants (g^, g , b^, bj to determine, and we do so by applying 
the surface conditions (34). 

At X = ^a , we get 

(38) * ' * • e-^^ + => b^^-*Ka + b_e*Ka ^ 

(39) < k(e-*»« -.g^e*^)'= AK(b^e-*Ka . ^yK^^ . 
similarly, at x = +a, ' - ^ ' ' 

(40) - b^e*Ka + b_e*'*^* = (1 + g+)e*^ . 

(41) ^ . KA(b^e*Ka . b_e-*Kaj = + g^)e*^ . 
fRfil we have four algebraic equations for the four uhknowns. ^ 



Solving these, and introducing the abbi^viatioft^ 

Z - 1 



(42) 
we obtain 

(43) . 

M) 



Q 



z -*^i ' 



^7 - KA'- A 



g = -Q 



1 -Q2ei4Ka 



^ 1 (1 -Q2)e^(K-k )2a 



1 - Q2 e 



14 Ka 



where R and T are called, the reflection and, transmission coefficients. The x5 or re 
sponcjing internal field is ^ ' * ^ ^ - ^ , 

• • . . • 1 - Q2e*4Ka . . 

Expanding tKe denominators in (43) and (44) enables us to interpret the solution in 
term^ of multiple reflections inside th^ slab. ' ^ i ! 

If *we are dealing with a s^lngte interface at x = 0 as in, Figure 9-3, theniWe 




FIGURE 9-3^ 
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ii . 



, ol^in simply 

•; • (46) . >^ 

■ (4t)- * 



E2 = (1 - Q)e*Kx ^ 



The results *and all the above may be generalized by inspection to an arbitrary , 
p angle of incidence a as in Figure 9-4. Thus in equatio^s (32) and (33) we may replace 




tor-. 



FIGURE 9-4 



by k^cos^a and by K^cos^ j3 'and.obtain the same final functions' in t^rms of 
the new constants kcos a and KcosjS , e.g., (46) be^o^mes 



.(48) 



gikxcosa . Qig-I 



lkx\osa 



Z' - 1 
Z» + 1 * 



21 - KAcosg^ 
kcos a ' 



* ,and (47) becomes 
(49) 



(1 - Q*)eJ'^'' , 



We itiay now multiply (48)* and (49) by. the same factor e^?^?. This converts 

(48) to ^ I 

(50y ; El - e**"'*^***"**^^*^"^- Qte-i>«c<»a<tifcysiiftx ^ Ei(ay r .Q^Ei(7r-a:) / , ^ 

• ^ - ' i, \ ' ' * ^* ' 

where ]^x(Oi) is a plane wave incident at an angle a, and Ei(7l^a).is its^mirror image* 

in the plane x = 0 . From (49), multiplication by e^^^*^"" gives 

- Qij>eiKxcU0+Ikysina ^ 



(51) 

If[we require that 



ksince = Ksin/? , ; i.e.,* 



K 



sin Of -r-sinfi.^ fi sin i3 



then (51) equals- 
(5^) 




^ E2 = .(1 - Q^)e*^*^^^**^y**"^ 



which is a plane wave traveling at an angle jS . Equation (52), which we recognize as 
"Sne.lUs La\v" [S], is thus an artifice for converting the solution of a one-dimensional 
problem to the corresponding two-dimensional solutioui it insures that corresponding, 
wave fronts match at the surf ape* , - * 

Integral representations ; In Sections 6 and 8 we used Huygens' principle in orde^to 
represent, the total scattered field as the integral of "wavelets" arising from elemen-- 
tary sources distributed over a surface or thrppghout a volume. To round out the 
previous intuitive discission we should indicate how such forms follow from^the 
present mathematical model [IJ^, [H], and.[m]. If \ye had ayailablei a theorem of 
Gauss (which relate^ certain surface and volume integral forms) we could prove that 
scattering functions satisfying [I] and {Illml can be represented' in terms of ele- 
mentary sources H as * 



(54) 



= /[H:„(klr-pl)a„E(p) - E(p)a„ 
s ^ 



H„(klr - pi)] dS(p) ./ 



where p Is 'a pgint on a surface S<p ) as in Figure 9-|5 that encloses the aoatterer but 
excludes the observation point f , and wher^n = ^ is the outward normal deriva- 
tive. The function ECp ) = Ej (p) + E,'(p ) , th^ total [field at p , and its normal deriv- 
ative, are weigljting factors for the surface distributi|)n of elementary sour^^ 
Hnj(lf -'pi) and^anHm which rdSdiate from p to r, i ^ 



r - P 




The elementary sources are essentia^l}^ those we worked with i^ earlier sections. 
Thus, in three-dimensions (i.e., if all three ^pace dimensions are significant in the 
p^blem), w^ require a pbiht source « i ' > ^ 

(55) H3(kR) = , R = \r-'p\ = /.(x - + (y - 7,)2 + (z - S)2 , 



and in one-dimension, a planar source 



(56) 



,UcR 



Hi(kR) = 2__^ Rflx.; 
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In hvo-dimen^ons, for argument kR» 1 , required line source approximates 
(57) H2(kR) ~ ^OTeiicR R- /r — 1,2 ^ r ;t> 

which is th« form ^e worked with previously; the fact that (57) holds only for large 
-values of the argument accounts for the restriction kR » l thai we mentioned ^pr the 
^strip andK^ylinder.problems. For sgiall values of J.R ,.the elementary line s^rcy 

behaves quite differently- ' ■ ' • . -^Z^ 

^ ^"^ ^ 

(58) 



1 . 
H2(kR) - 2^lrikR , Off, 



^' : Its exapt representation Is given by 

Where. Ho (which is knbwn as "HankeHs function of the first kind of order zero) 'is the 
special solution of (20) for two-dimensions corresponding to angle independent out- 
gomg waves, i.e" it plays the same fole in two-dimengjon^s ^ and e'""" piay in 

. the other cases. If we specialize S(p ) to the surfacex^f-th^catterer its'elf , then we 
can use such surface conditiois as [IIJ to obtain equ^ions (integral equations) for the 

. unknown values of E,(p) and a,E,(p) • for simple surfaces-, the procedure is analo- 
gous to that we followed for the slab* . - . " ' 

Although we will not 'prove (54), we wilf show how totbbtain the approximate forms 
, we worked with in earlier sections. Thus if we specialize S( p) to the scatterer's . 
■ .surface^and use th.e boundary condition [HbJ that a„E'(p) = 0 , we reduce <54) to 

- - .:(^^) • E,(r) - -/E(p)a„H„,'(kl?-?|)dS(^) ' 

In particular, in tWo-dimensions and kfr - pi » l , wetise (57) in (60) to^btain ' 

' which is of the required, form [nic J . • If we knew the field E(p ) on the seattWer's sur- ' 
face, we could obtain the scattered field E,(r ) by integration. If we do not know the 
• .field (and it is only for very siniple shapes that E(?) is knpwn exactly), then we may 
^ seek heuristic physically motivated approximations. 1 

;n particular if the scatterer is very big compared to iavelength , then it is • 
plausible to, approximate E(p) by elementary geometrical Utics considerations. , 
Following essentially Kirchhoff, one approximates the totaljfiey E(p) on the "lit' ' 
Side" of the s^rer by twice the incident value E. , and bV zero oi> the "dark ^ide." 

Thus if Qiihcfi'fiif/^ ^ ; 




Thus if we substitute 

'(62) ^ E(p)«B.(p) on lit side"; • E(p>« 0. oji dark side 

^" / ' 1'48 
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into (61) we obtain the general case of the integrals we considered jn Sections 6 and 8, 
i.e., for the stFip with dS(p) = dj] and the circular cylinder with dS(p) = ad(p . 

From (54) we could also construct the volume integral of elementary scatterers 
that we worked with previously for the case of a partially transparent sphere. First 
we specialize (54) to p on the scatterer and then use the transition conditions [lie] 
to replace the external surface fields E(p) = E(k,p) and ^nMp) by the correspond- 
ing internal fields E(K,p)' and A8nE(K,p) where K - k/ii is t^e^internal wave number. 
We then use^ the same theorem of Gauss to convert the resulting surface integral to an' 
integral over the volume of the scatterer. In particular for constant ^ and A = 1 , 
we would obtain ' 



(63) 



E, = (k2-K2)/H^(gf^^^p|)E(K,p)dV(p) 



where V is the volume of the scatterer. If we add Ej to both sides of (63) we obtain 
an integral equation for E which can be solved for simple shades. We will not prove 

(63) , but we Will show how this rigorous result leads to the previous appro^amation 
(6:38). - ; ' ^ 

For tenutus scattdrqrs in the sense = ^ k2 , Rayleigh replaced the 
unknown internal field E(p) in (6^ by the incident wave: ^ ' 

(64) ' m.py^ Ei(k,p)'- e**^^ . . ' 

If we substitute j(64) into (63) and specialize to three -dimensions by using (55), 



we 



obtain 



P ^. kUu^-i) r e"^"^"^' ;,c^,,,., 
Es. ""JT^J -|^T^e*»^^dV(p) 



which is the more complete version of the form wq worked with previously/{n (6:38)ff . 
Equation (65) gives. directly the result that the flux (IeJ^) is inversely proportipnal 
^ to \^ for various scatterers whose length dimensions are small cotnpareti to X , and 
was used by Rayleigh prior to [1r] of Section 7. 



i 



It should be stressed th^t such Approximations as (62) and (64) are adequate only , 
for limit^ ranges of the parameters. However,. wWn their limitations, they provide^ 
u^efuj and instilietive explicit results for probld!ns that c^ot be solved rigorously. 



ERIC 



, 1-' 



143 1 4 9^ 



14 



In concluding this chapter, we should reiterate that we have covered merely a 
Selected sequence of topics in wave physics. The wave equations that we*^t)acked into" 
are, in general, geherated in physics prog^:ams by operating on first order differentiat 
equations that relate the physical obseryables regarded as basic within a particular 
discipline, e.g., particle^ velocity and excess pressure in acoustics, electric and 
magnetic intensities in optics, radio, etc. The wave function E We have dealt with ^ 
represents one of theSe different physical obse^ables (or one of its Cartesian com- 
ponents), and may also stand for the probability density function of quantum mechanics. 
Similarly the media that we specified by an index of refraction m and interface, condi- 
tions represent quite different concepts in different disciplines, and involve different 
kinds of physical parameters, i^nplicit continuity rec^uirements on appropriate pljy^- 
cally observable fields, etp. We have discussed neither the physics implicit in the 
above nor in the much more complex question of sources ard the geaeration of fields. 
We began with let there be light, and followed a narrow thread of concepts. 

Although we made "light" the theme for much of the development, we have not 
* covered an essential aspect that di^stinguishes wave models for light from the models 
used for sound: light, and all electromagnetic waves,. must also be characterized by 
polarization , this requires in general that we deal with^vectbr wave functiohs with 
amplitudes perpendicular to the direction of propagation instead of the scalar functio^ 
we have considered. However, out discu^ssion of light >v as in no sense meant to be 
comprehensive, and as stressed in the introduction of this chapter, there are many 
phenomena involving* light that are not described by a wave model at all. In' illustrat- 
ing different applications of calculus, we have used light as a vehicle for an introduc- 
tion to wave physic §, not only because we have many visual experiences to draw on, 
but because the adequacy of ^a wave model for such phehomepa was far from oQifious 
to the early investigators (and not particularly obvious even to us. without some care- 
fill observations). , For water waves, the appropriateness of the ma^epjatical model . 
would have ^een cle^r frotn the start, and even for sound waves the intuition leads 
relatively directly from the visible waves on stringed instruments and on drum heads 
, to waves 4n air. Thus in discussii^ light, we could introduce key topics leading to the 
development of the wave model essentially in th^ir historical order, and thereby indi- 
cate the, greater generality of'the wave itiodeL as v^ell as the domains of applicabilit)^ 
^of the ^arlier i^pecial "laws of nature" that are now exhibited as consequences. How- 
,evet, initial reservation's that "light" is neither wave nor particle, and tfiat only 
.certaii^ classes of phenomena involving light are adequately ^described by a^yave 
model jsnould not l?e lost sights of. Light is one of the most complex characters in the 
J^BC of ^mathemajtical physics. * ' ' ( , ' . ' ^ 



The preceding two chapters illustrate attempts at systematic ^approa'^hes to appli- 
cations of mathematics in sci^ence. Chapter 2 considers simple -equations for growth 
and competition that arise, again and again in superficially unrelated, cdntexts of nature, 
it shows that phenomena and processes that occur in all the sciences are linked by one 
mathematical modeL We select a narrow thread of mathematical njethods and follow 
it through various sciences. Chapter 3 is qujte different; there we follow science as a 
thread.. We select a narrow sequence of physical cpncepts leading from geometrical 
optics through wave physics and exhibit various methods of the calculus that further 
the development.^ ^ ^ - , 

Thus our two chag^ll^^iwnatjiefhatics and science atre very different. They 
supplement each other i;j[^l^jjjicating the ways that mathematics and science interact. 
The first chapter follows £^^athematics thread, the second a science thread, and the 
two together may suggest the crossing threads of a fabric, two threads intersect 
at Rayleigh^s theory for the color of the sky: in Chapter 2, it is a special case of a 
general attenuation process; in Chapter 3, a spegial case of general scatt.ering 
process, 

* One thread suggests that mathematiqs }ntef sects every science. The other thread 
suggests that every science intersects all of the mathematics. Together they mj^y 
suggest that the interactions of mathematics and science are profound indeed^ These 
aLJljd the threads of the ^fabric of our universe, the structure of our perception of nature. 
Mathematics and science are the very warp and woof of the universe our intellect has 
created, ~ - - < - • - » < . - - . • 
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